level 12
黄金分割法（0.618法）：
function y=hjfg(a,b)
n=0;
while (norm(b-a)>1.0e-6)
x1=a+0.382*(b-a);
x2=a+0.618*(b-a);
y1=x1^3-2*x1^2+1;
y2=x2^3-2*x2^2+1;
if y1>y2
a=x1;
else
b=x2;
end
n=n+1;
end
x=(b+a)/2.0;
x
n
黄金分割法（0.618法）：
function y=hjfg(a,b)
n=0;
while (norm(b-a)>1.0e-6)
x1=a+0.382*(b-a);
x2=a+0.618*(b-a);
y1=x1^3-2*x1^2+1;
y2=x2^3-2*x2^2+1;
if y1>y2
a=x1;
else
b=x2;
end
n=n+1;
end
x=(b+a)/2.0;
x
n
level 12
level 12
level 12
level 12
level 12
level 12
level 12
level 12
level 12
level 12
level 12
level 12
level 12
level 12
level 12
level 12
level 12
level 12
level 12