C++. Метод золотого сечения для поиска екстремумов

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

#define epsilon 0.0001
//------------------------------------------------------------------------------

double func(double x){
return x*x*x + 10*x*x + 31*x +30;
}

//------------------------------------------------------------------------------

double double_find_zp_min(double a, double b, double eps){

if ( a > b )
{
double temp = a;
a = b;
b = temp;
}

while ( fabs( a - b ) > eps )
{
if (
func( a + fabs(b-a) * 0.38)

< // МІНІМУМ

func( b - fabs(b-a) * 0.38)
)
{
b -= fabs(a-b) * 0.38;
}
else
{
a += fabs(a-b) * 0.38;
}
}

return (a+b)/2;
}

//------------------------------------------------------------------------------

double double_find_zp_max(double a, double b, double eps){
if ( a > b )
{
double temp = a;
a = b;
b = temp;
}

while ( fabs( a - b ) > eps )
{
if (
func( a + fabs(b-a) * 0.38)

>

func( b - fabs(b-a) * 0.38)
)
{
b -= fabs(a-b) * 0.38;
}
else
{
a += fabs(a-b) * 0.38;

}
}

return (a+b)/2;
}

//------------------------------------------------------------------------------

int main(int argc, char* argv[])
{

double min_x=double_find_zp_min(-3, -2, epsilon);
double max_x=double_find_zp_max(-5, -3, epsilon);
printf(" Min of func ( %.4f ) = %.4f \n",min_x, func(min_x));
printf(" Max of func ( %.4f ) = %.4f \n",max_x, func(max_x));


system("pause");
return 0;
}
вот написал кое какой код...я пересчитовал на бумаге и оказалось что прога считает почему то неправильно.. и еще препод требует чтоб а и б вводили сами ... то есть чтобы поиск екстремумов проводился на любом промежутке заданом нами..буду благодарен за любую помощь в коректировке кода...