%采用fib序列法，單變量搜索
a=0;
b=1;
e=0.0005;
M=100;
%l=(a+b)/2;
xbest=0;
ybest=0;
tic
c1=1;
c2=2;

for i=1:M
    c3=c1+c2;
    c=c1/c3;
    
    c1=c2;
    c2=c3;
    
    l=b-a;
    if l<e
        break;
    end
    x1=a+c*l;
    
    x2=b-c*l;
    y1=f(x1);
    y2=f(x2);

    %for i=1:10
    if(y1>y2)
       a=x1;
       ybest=y2;
       xbest=x2;
    else
       b=x2;
       ybest=y1;
       xbest=x1;
    end
end

i
l
xbest
ybest
toc