function podzial()
alfa=(sqrt(5)-1)/2;
epsilon=0.001;
a=0;
b=5;
xL=b-alfa*(b-a);
xR=a+alfa.*(b-a);
ip=0;
ik=10;

while((b-a)>epsilon || ip==ik)
    if(f(xL)<f(xR))
        b = xR;
		xR = xL;
		xL = b - alfa.*( b - a );
    else
        a = xL;
		xL = xR;
		xR = a +  alfa.*( b - a );
    end
    ip=ip+1;
end
xMin = (a+b)/2;
disp(xMin);
disp(f(xMin));
end

function [y] = f(x)
    y = x.^3-6*x.^2+11*x-6;
end