[sharonyue]

Hi guys,

I m doing simulation of PBM. In QMOM, all the paper suggest to use PD algorithm. http://www.tandfonline.com/doi/abs/1...86829708965471
But I failed. So I wanna try another way. the problem Im facing is:



How can I get r1,r2,r3,w1,w2,w3? Thanks.

[sharonyue]

solved:

function co
x0=[1 2 3 4 5 6];
options=optimset('display','off');
[x,fval]=fsolve(@han,x0)


function y=han(x)
y(1)=x(1)+x(2)+x(3)-1;
y(2)=x(4)*x(1)+x(5)*x(2)+x(6)*x(3)-6.417994;
y(3)=x(4)^2*x(1)+x(5)^2*x(2)+x(6)^2*x(3)-48.9333;
y(4)=x(4)^3*x(1)+x(5)^3*x(2)+x(6)^3*x(3)-426.86446;
y(5)=x(4)^4*x(1)+x(5)^4*x(2)+x(6)^4*x(3)-4255.826960;
y(6)=x(4)^5*x(1)+x(5)^5*x(2)+x(6)^5*x(3)-48347.809106  ;
y=[y(1) y(2) y(3) y(4) y(5) y(6)];