function x=fsxgold(f,e,a,b)
% fsxgold函數 黃金分割算法 
% x=fsxgold(f,e,a,b)為輸入函數 f為函數 e為允許誤差 (a,b)為搜索空間;
L=(5^(1/2)-1)/2;
u=a+L*(b-a);f2=subs(f,u);t=a+b-u;f1=subs(f,t);
while (abs(u-t)>=e)
    if (f1<=f2)  b=u;u=t;f2=f1;t=a+b-u;f1=subs(f,t);
    else  a=t;t=u;f1=f2;u=a+L*(b-a);f2=subs(f,u); 
end;
end;
x=0.5*(t+u);