одномерная оптимизация методом дихотомии

[DiGris]

имею код,// Optimiz dihoto.cpp: определяет точку входа для консольного приложения.
//
#include <stdafx.h>
#include <iostream>
#include <conio.h>
#include <cmath>
#include <stdio.h>
using namespace std;
double f(double x)//объявляем функцию f
{
double f;
f = pow(x, 2) + 8*exp(0.55*x);//pow вычисляет х возведенный в степень у, в данном случае х возводится в квадрат
return f;
}

void main()
{
const int m=1000;//задаем максимальную координату точки минимума
const double eps=0.001,d=eps/10;//задаем точночть
int i,j,l,N,n;
double a,b,x[m],y[m],yl,xl,X,Y;
N=0;//N-количество разбитых отрезков
cout<<"a=";cin>>a;
cout<<"b=";cin>>b;
while(b-a>2*eps)//задаем границы разности координат отрезка
{
X=(a+b)/2;//нахождение абсцисы точки минимума
x[1]=X-d; y[1]=f(x[1]);
x[2]=X+d; y[2]=f(x[2]);
N=N+2;
if(y[1]<y[2])
b=x[2];
else a=x[1];
//a=x-d;
//b=x+d
}
X=(a+b)/2;//искомая точка по Х
Y=f(X);///искомая точка по У
cout<<"x="<<X<<" y="<<Y<<" N="<<N;
getch();
}
работает нормально, но препод говорит изменить его под алгоритм из учебника... заранее спасибо)

[Serik2010]

Погляди вот такую реализацию дихотомии на С++ тут
А тут параллельно еще и золотого сечения метод тут