function [ a0 b0 ] = zlatni( funkcija,a0,b0 ,epsilon)
flag=0;
n=0;
c=(3-sqrt(5))/2;
while abs(a0-b0)>epsilon
    x1=a0+c*(b0-a0);
    x2=a0+b0-x1;
    fx1=feval(funkcija,x1);
    fx2=feval(funkcija,x2);
    if fx1<=fx2
       b0=x2;
    else
        a0=x1;
    end;
end;
end