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

float pauel (double e,double a,double b)
{
int k=1;
double c,d;
c=b;
b=(a+b)/2;
d=((a+b)/2+0.5* (((func(a)-func(b))*(b-c)*(c-a))/ (func(a)*(b-c)+func(b)*(c-a)+func(c)*(a-b))));
while ( (fabs((d-b)/2)>=e) || (fabs((func(d)-func(b))/func(b))>=e) )
{
if (b<d)
{
if (func(b)<func(d))
c=d;
else
{
a=b;
b=d;
}
}
else
{
if (func(b)<func(d))
a=d;
else
{
c=b;
b=d;
}
}
k++;
d=((a+b)/2+0.5* (((func(a)-func(b))*(b-c)*(c-a))/ (func(a)*(b-c)+func(b)*(c-a)+func(c)*(a-b))));
}
cout<<k<<endl;
return ((b+d)/2);
};