Метод золотого сечения.

#include <iostream>
#include <math.h>
using namespace std;
 
double f(double i)
{
    return  -1*(i*i*i)+9*(i*i)-24*i-2;
}
int main()
{
int i, n;
cout<<"Kol-vo shagov n= ";
cin>>n;
    double a=-5, b=3, g, h, fu1, fu2;
    g = a+(b-a)*(3-sqrt(double(5)))/2;
    h = a+(b-a)*(sqrt(double(5))-1)/2;
    fu1 = f(g);
    fu2 = f(h);
    for(i = 1; i <= n; i++)
    {
        if( fu1<=fu2 )
        {
            b = h;
            h = g;
            fu2 = fu1;
            g = a+(b-a)*(3-sqrt(double(5)))/2;
            fu1 = f(g);
        }
        else
        {
            a = g;
            g = h;
            fu1 = fu2;
            h = a+(b-a)*(sqrt(double(5))-1)/2;
            fu2 = f(h);
        }
    }
    cout<<(fu1-fu2)/2+fu2<<endl;
system ("PAUSE");
return 0;
}

/*
 Метод Золотого сечения (отношения)
 Программу сделал Сеченов Слава
Разработана на Borland C++
Если вы прогите в другой програмулине, 
то удалите строки где в коментах написании // Borland C++
   Удачи.
 */
#include <vcl.h>// Borland C++
#pragma hdrstop// Borland C++
    #include <math>
   #include <iostream>
#define M_E 2.71828182845904523536 //в Borland и без этой строки работает
   using namespace std;
 
double f(double x)
{return  pow(M_E,x)*sin(x*x-2*x+11);} //Здесь ваша формула
 
//---------------------------------------------------------------------------
 double ab(double i)
 {if(i<0)return -i;else return i;}   //это модуль
 
//---------------------------------------------------------------------------
 
#pragma argsused// Borland C++
void main()
{
 double a, b, g, h, fu1, fu2,eps=0.000001;
 
 
        for(int i=-1;i<=1;i+=2)//Сначало ищем минимум, затем максимум
        {
 
        cin>>a>>b; //Вводим начальную координату и конечную
 
    g = a+(b-a)*(0.38);
    h = b-0.38*(b-a);
    fu1 = i*f(g);
    fu2 = i*f(h);
 
    do{
 
        if( fu1<=fu2 )
        {
            b = h;
            h = g;
            fu2 = fu1;
            g = a+(b-a)*(0.38);  //Основная формула
            fu1 = i*f(g);
        }
        else
        {
            a = g;
            g = h;
            fu1 = fu2;
            h = b-0.38*(b-a);  //Основная формула
            fu2 = i*f(h);
        }
    }
    while(ab(a-b)>eps); //Условие выхода
    i==-1? cout<<"Minimum ":cout<<"Maximum ";
        cout<<"OTBET "<<(fu1-fu2)/2+fu2<<" v tochke "<<h<<endl;
 
 }
 
system ("PAUSE"); //Консольная команда паузы проги
 
 
}