Question

想用Matlab求解方程组t = 1 - (1-p)^(1/4) ;t = 2*(1-2*p) /((1-2*p)*33

Answer 1

%main.m

x0=[0 1];

options=optimset('Display','iter');

[x,fval] = fsolve(@myfun,x0,options) ;

p = x(1)

t = x(2)

以下存为myfun.m

function f=myfun(x)

f = [1-(1-x(1))^(1/4)-x(2); 2*(1-2*x(1))./((1-2*x(1))*33 + x(1)*33*(1-(2*x(1))^5))-x(2)];

end

运行main.m

p =

0.1773

t =

0.0476