Метод деления пополам (дихотомии)

[Килямов Сергей]

вот что я написал:
//---------------------------------------------------------------------------

#include <vcl.h>
#pragma hdrstop

#include "dihotomia.h"
#include "math.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
: TForm(Owner)
{
}
//---------------------------------------------------------------------------

void __fastcall TForm1::Button1Click(TObject *Sender)
{
double a, b, x1, x2, y1, y2, eps, e, d, En, xmin, xmax, ymax;
int i = 0 ;
eps = 0.001 ;
e = 2.7 ;
a = 0 ;
b = 1 ;
d = 0.001 ;
do
{x1 = (b + a - d) * 0.5 ;
x2 = (b + a + d) * 0.5 ;
y1 = 1 / pow(e, x1) + pow(e, (-2 * x1)) ;
y2 = 1 / pow(e, x2) + pow(e, (-2 * x2)) ;
if (x1 >= x2)
{b = x2;
y2 = 1 / pow(e, x2) + pow(e, (-2 * x2)) ;
}
else
{a = x1;
y1 = 1 / pow(e, x1) + pow(e, (-2 * x1)) ;
}
En = (b - a) * 0.5 ;
xmax = (a + b) * 0.5 ;
}
while(( b + a) * 0.5) ;
ymax = 1 / pow(e, xmax) + pow(e, (-2 * xmax)) ;

Edit1->Text = xmax ;

Edit2->Text = ymax ;
}

//---------------------------------------------------------------------------
Ничего не считает кто что посоветует? писал в C++Builder 6

[Vago]

О, боги! Вы что вообще ищете? Корень или экстремум?

[Килямов Сергей]

ищу максимум функции

[Vago]

#include <math.h>
#include <stdio.h>

#define DEBUG_PRINT 0

double fn( double x ) {
   return ( 1. / (exp( x ) + exp( -2. * x ) ) ) ;
}

int main(){

   double   a = 0. , 
            b = 1. , 
            x1, x2, xStar,
            y1, y2, yStar ,
            eps = 0.001, 
            delta = 0.0001 ;
 
   while ( (b-a) > eps ) {
      x1 = (b + a - delta ) / 2. ;
      x2 = (b + a + delta ) / 2. ;
      y1 = fn( x1 ) ;
      y2 = fn( x2 ) ;
      if ( y1 >= y2 )
         b = x2 ;
      else
         a = x1 ;
#if DEBUG_PRINT == 1
      printf( "%.4g %.4g\n", a, b ) ;
#endif
   }
   xStar = (a + b ) / 2 ;
   yStar = fn( xStar ) ;

   printf( "Max. value F(%.3g) = %.3g\n", xStar, yStar ) ;

   return 0 ;

}

Привет составителям вашей методички, написавшим в "Шаге 2" действия, предпринимаемые при поиске МИНИМУМА и ни слова об этом не сказавшим. (О том, что весь этот так называемый "алгоритм" попросту содержит формальные ошибки, ибо, в реальности, на каждом шаге a либо b должны меняться - деликатно умолчим).

Преподавателям программирования, не рассказавшим вам про старшинство операций, и не сообщившим, что для вычисления "e в степени" в Си существует ф-я exp() - отдельный привет!

Ввод вместо присваиваний и хотя бы элементарный контроль введённых данных - уж сами, пожалуйста...