метод золотого сечения!не могу найти ошибку в коде(С++)

[prostdj]

#include <math.h>
#include <iostream>
using namespace std;

double f(double x)
{
return x*x*x+10*x*x+31*x+30;
} 
double max(double a, double b, double eps)
{

double g, h, fu1, fu2;
double fi = (1+sqrt(5.0))/2;
//cout<<"Input a for max"<<endl;
//cin>>a;
//cout<<"Input b for max"<<endl;
//cin>>b;
g = (a+(b-a))/fi;
h = (b-(b-a))/fi;
fu1 = f(g);
fu2 = f(h);

do{

if( fu1<=fu2 )
{
b = h;
h = g;
fu2 = fu1;
g = (a+(b-a))/fi;
fu1 = f(g);
}
else
{
a = g;
g = h;
fu1 = fu2;
h = (b-(b-a))/fi;
fu2 = f(h);
} 
}
while(fabs(b-a)<eps);
cout<<"Maximum "<<f((a+b)/2)<<" v tochke "<<h<<endl;
} 
double min(double a, double b, double eps)
{

double g, h, fu1, fu2;
double fi = (1+sqrt(5.0))/2;
//cout<<"Input a for min"<<endl;
//cin>>a;
//cout<<"Input b for min"<<endl;
//cin>>b;
g = (a+(b-a))/fi;
h = (b-(b-a))/fi;
fu1 = f(g);
fu2 = f(h);

do{

if( fu1>=fu2 )
{
b = h;
h = g;
fu2 = fu1;
g = (a+(b-a))/fi;
fu1 = f(g);
}
else
{
a = g;
g = h;
fu1 = fu2;
h = (b-(b-a))/fi;
fu2 = f(h);
} 
}
while(fabs(b-a)<eps);
cout<<"Minimum "<<f((a+b)/2)<<" v tochke "<<h<<endl;
} 
int main()
{
double a, b,eps;
cout<<"Input a"<<endl;
cin>>a;
cout<<"Input b"<<endl;
cin>>b;
cout<<"Input eps"<<endl;
cin>>eps;
max(a,b,eps);
min(a,b,eps);
system ("PAUSE");
return 0;
}

вот написал кое какой код...препод говорит что неправильно считает.... прога работает..писал на dev cpp
буду благодарен за любую помощь в коректировке кода...