Градиентные методы / zs2

int ZS_2(double a, double b, double *x, double e)
{
long k=1;
double n618=((-1)+pow(5,0.5))/2;
double x1, x2;

x1=a+n618*(b-a);

while((b-a)>e)
{
x2 = a+b-x1;

if(f(x1)<f(x2))
{
if(x1<x2)
b=x2;
else
a=x2;
}
else
{
if(x1<x2)
a=x1;
else
b=x1;
x1=x2;
}
k++;
};
*x=(a+b)/2;
return k;
}