function [xmin, fmin] = Untitled2( f,a,b,e,alfa )
l=a+(1-alfa)*(b-a);
u=a+alfa*(b-a);
k=0;
hold on
while abs(a-b)>e
    if feval(f,l)>feval(f,u)
        a=l;
        l=u;
        u=a+alfa*(b-a);
    else
        b=u;
        u=l;
        l=a+(1-alfa)*(b-a);
    end
    plot([k,k],[a,b]);
    k=k+1;
end
xmin=(a+b)/2;
fmin=feval(f,xmin);
end