【无标题】
f=@(x)x^2-3*x-4;
a=0;b=2;
eps=0.05;
r=a+0.382*(b-a);u=a+0.618*(b-a);
while 1
if(feval(f,r)>feval(f,u))
if((u-r)<eps)
f=feval(f,u);
disp(‘f’);
break;
else
a=r;
r=a+0.382*(b-a);
u=a+0.618*(b-a);
end
else
if((u-r)<eps);
f=feval(f,r);
disp(‘f’);
break
原创
2021-12-19 12:48:46 ·
51 阅读 ·
0 评论