function [x,b,yx]=golden1(a,b,l) functionyx=fun(x)
yx=x^2-2*x+1;
end
a=0;
b=3;
l=0.15;
n=0;
x1=a+0.382*(b-a);
x2=a+0.618*(b-a);
y1=fun(x1);
y2=fun(x2);
while (b-a)/3>l
if y1>y2
a=x1;
x1=x2;
y1=y2;
x2=a+0.618*(b-a);
y2=fun(x2);
else
b=x2;
x2=x1;
y2=y1;
x1=a+0.382*(b-a);
y1=fun(x1); end
n=n+1;
end
n
x=(b+a)/2
yx=fun(x)
end
function [x,b,yx]=golden1(a,b,l) functionyx=fun(x)
yx=x^2-2*x+1;
end
a=0;
b=3;
l=0.15;
n=0;
x1=a+0.382*(b-a);
x2=a+0.618*(b-a);
y1=fun(x1);
y2=fun(x2);
while (b-a)/3>l
if y1>y2
a=x1;
x1=x2;
y1=y2;
x2=a+0.618*(b-a);
y2=fun(x2);
else
b=x2;
x2=x1;
y2=y1;
x1=a+0.382*(b-a);
y1=fun(x1); end
n=n+1;
end
n
x=(b+a)/2
yx=fun(x)
end