4Â¥: Originally posted by crazyiyy at 2014-07-14 16:39:24
ÄúºÃ£¡ ÄúµÄ´ð°¸µÄÈ·¸ú1stOptÒ»Ö£¬»¹ÓиöÎÊÌâÏ£ÍûÄúÄÜΪÎÒ½â»ó¡£Èç¹ûÎÒÐèÒª¶Ô·½³Ì×éÖеÄijЩϵÊý½øÐб仯£¬½ø¶øµÃµ½Ò»×é½â£¬ÎÒÓ¦¸ÃÈçºÎ²Ù×÷¡£
ÀýÈç ÎÒÐèÒª¶Ô·½³Ì1ÖеÄ48*2^3 ±äΪ48*2^x  ÆäÖÐx=1,2,3,4 ÕâÑùÎÒÓ¦ ...

function jie_fun
u0=[0.5 0.5 0.5];
u=fsolve(@myfun,u0)
f_val=myfun(u)



function f=myfun(u)
x=u(1);
y=u(2);
z=u(3);
f(1)=1/((4.5+(8.5)*(x+y-x*y)+(33/2)*(x+y-x*y)^2+(33/2)*(x+y-x*y)^3+(33/2)*...
    (x+y-x*y)^4)*(1-(x+y-x*y))/(1-(x+y-x*y)^5)+(1-(12/(48*2^3-4)))*(1-x)+6*...
    (12/(48*2^3-4))+(1-(12/(48*2^3-4)))*(1-x)*(12*(1-y)+1))-z;
f(2)= (13*(1-(1-z)^8)+9*z*(1-z)^8)*(1-x)*(1-y)-x;
f(3)=(1-(1-z)^8)/(2-(1-z)^9)-y;

function jie_fun
clear a;clc
format long
u0=[0.5 0.5 0.5];
u=fsolve(@myfun,u0)
f_val=myfun(u)



function f=myfun(u)
x=u(1);
y=u(2);
z=u(3);
f(1)=1/((4.5+(8.5)*(x+y-x*y)+(33/2)*(x+y-x*y)^2+(33/2)*(x+y-x*y)^3+(33/2)*...
    (x+y-x*y)^4)*(1-(x+y-x*y))/(1-(x+y-x*y)^5)+(1-(12/(48*2^3-4)))*(1-x)+6*...
    (12/(48*2^3-4))+(1-(12/(48*2^3-4)))*(1-x)*(12*(1-y)+1))-z;
f(2)= (13*(1-(1-z)^8)+9*z*(1-z)^8)*(1-x)*(1-y)-x;
f(3)=(1-(1-z)^8)/(2-(1-z)^9)-y;

3Â¥: Originally posted by ÔÂÖ»À¶ at 2014-07-14 16:14:10
ÒÔϳÌÐòµÚ¶Ô2Â¥ÖгÌÐòµÄ¸Ä½ø£¬Ìá¸ß½âµÄСÊýµãºóµÄλÊý¡£


function jie_fun
clear a;clc
format long
u0=;
u=fsolve(@myfun,u0)
f_val=myfun(u)



function f=myfun(u)
x=u(1);
y=u(2);
z=u(3);
...