黄金分割法matlab程序

自定义函数：
function y=koko(x)
y=x*(x+2);

主程序：
a=input('请输入初始区间下端点:\na=');
b=input('请输入初始区间上端点:\nb=');
e=input('请输入计算精度:\ne=');
t=b-a;
while t>e
a1=a+0.382*(b-a);
a2=a+0.618*(b-a);
f1=koko(a1);
f2=koko(a2);
if f1b=a2;
else
a=a1;
end
t=b-a;
end
X1=(b+a)/2;
F1=koko(X1);
fprintf('最优解为：\nX1=%8.6f,\nF1=%8.6f',X1,F1);

运行结果：
请输入初始区间下端点:
a=-3
请输入初始区间上端点:
b=5
请输入计算精度:
e=0.3
最优解为：
X1=-0.973876,
F1=-0.999318>>