- 积分
- 0
- 注册时间
- 2011-5-23
- 仿真币
-
- 最后登录
- 1970-1-1
|
请高手指点啊,下面的程序解不到想要的结果,请大家帮忙看看!
%不用m文件法;
R0=zeros(1,281);
n=zeros(1,281); ;
Hv=zeros(1,281);;
syms H;
for k=1:281,
H=k+199;
[email=f=@(x)([((0.08*log(exp(x(1))+0.002))/(x(1)*exp(x(1))-0.998*x(1)))^x(2)-H*(exp(x(1))+0.002)/(630000*x(1]f=@(x)([((0.08*log(exp(x(1))+0.002))/(x(1)*exp(x(1))-0.998*x(1)))^x(2)-H*(exp(x(1))+0.002)/(630000*x(1[/email]));
210000*x(1)*(x(2)/(x(1)*exp(1)))^x(2)-3*H-22.3]);
options=optimset('MaxFunEvals',1000);
x=fsolve(f,[0.0012;0.152],options); |
|