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

void ZS1()
{
int k=1;
double m,l,e=0.001;
m=a+0.618*(b-a);//a и b задаются глобально
l=a+0.382*(b-a);
do
{
if (fx(l)<fx(m))//ф-ция в кот ищется мин.
{
b=m;
m=l;
l=a+0.382*(b-a);
}
else
{
a=l;
l=m;
m=a+0.618*(b-a);
}
k=k+1;
}
while((b-a)>=e);
cout << endl << "Minimum: " << (a+b)/2 << " k=" << k << endl << endl;
}