【无标题】

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;
else
b=u;
r=a+0.382*(b-a);
u=a+0.618*(b-a);
end
end
end