function y=fungold(a)
y=a^2+2*a;
end
function goldmin
clear;
clc;
a=-3;
b=5;
tolerance=0.01;
n=0;
i=100;
for k=1:i
if(abs(b-a)<=tolerance)
break;
else
a1=b-0.618*(b-a);
a2=a+0.618*(b-a);
if (fungold(a1)>=fungold(a2))
a=a1;
else
b=a2;
y=fungold((b+a)/2);
end
n=n+1;
end
end
n,a0=(b+a)/2,y=fungold((b+a)/2)
end