Реализация метода хорд по блок-схеме

//---------------------------------------------------------------------------
 
#include <vcl.h>
#pragma hdrstop
 #include<math.h>
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
    : TForm(Owner)
{
}
//---------------------------------------------------------------------------
void __fastcall TForm1::Button1Click(TObject *Sender)
{
 
float  a=2,b=3,e,y1,y2,x,y3,y;
 e=StrToFloat(Edit1->Text);
 y1=3*sin(sqrt(a))+0.35*a-3.8;
 y2=3*sin(sqrt(b))+0.35*b-3.8;
 
 
 x=a-y1*((b-a)/(y2-y1));
 y3=3*sin(sqrt(x))+0.35*x-3.8;
 
 if ((abs(x-a))<=e)  Memo1->Lines->Add(FloatToStr(x));
 else
  {
    if (y3<=e)  Memo1->Lines->Add(FloatToStr(x));
    else
     {
       if (y1*y3<0) { b=x; y2=y3;
 
       }  else {a=x; y1=y3;}
     }
 
  }
 
 
 
}

class Tricky_form : public TForm
{
__published: // IDE-managed Components
 
    TTimer * tmr_1;
 
    void __fastcall FormPaint( TObject * Sender );
    void __fastcall tmr_1Timer( TObject * Sender );
 
public: // User declarations
 
    double a0, b0, a, b, eps;
    double x_scale, y_scale;
 
    int x_scr( double x ) // функции для вывода на экран - масштабируют координаты точек для формы и инвертируют ось y
    {
        return ClientWidth / 2 + x * x_scale;
    }
 
    int y_scr( double y )
    {
        return ClientHeight / 2 - y * y_scale;
    }
 
    double f( double x )
    {
        return
            // pow( x, 4 ) - 13 * pow( x, 2 ) + 36;
            3 * sin( sqrt( fabs( x ) ) ) + 0.35 * x - 3.8;
    }
 
    __fastcall Tricky_form( TComponent * Owner );
} ;

__fastcall Tricky_form::Tricky_form( TComponent * Owner )
    :
    TForm( Owner )
{
    a0 = a = -7;
    b0 = b = 10;
    eps = 0.001;
 
    x_scale = ( ClientWidth / 2 ) / std::max( fabs( a0 ), fabs( b0 ) );
 
    double y_max( 0 );
    for ( double i( a0 ), i_limit( b0 ), y; i < i_limit; i += 0.1 )
    {
        y = fabs( f( i ) );
        y_max < y ? y_max = y : 0;
    }
    y_scale = ( ClientHeight / 2 ) / y_max;
}
 
void __fastcall Tricky_form::FormPaint( TObject * Sender )
{
    Canvas->Pen->Width = 1;
 
    // оси координат
    Canvas->Pen->Color = clGreen;
 
    Canvas->MoveTo( 0, ClientHeight / 2 );
    Canvas->LineTo( ClientWidth, ClientHeight / 2 );
    Canvas->MoveTo( ClientWidth / 2, 0 );
    Canvas->LineTo( ClientWidth / 2, ClientHeight );
 
    // график функции
    Canvas->Pen->Color = clBlue;
 
    Canvas->MoveTo( x_scr( a0 ), y_scr( f( a0 ) ) );
    for ( double i( a0 ), i_limit( b0 ); i < i_limit; i += 0.1 )
    {
        Canvas->LineTo( x_scr( i ), y_scr( f( i ) ) );
    }
 
    // хорда по текущим концам отрезка
    Canvas->Pen->Color = clRed;
    Canvas->Pen->Width = 3;
    Canvas->MoveTo( x_scr( a ), y_scr( f( a ) ) );
    Canvas->LineTo( x_scr( b ), y_scr( f( b ) ) );
}
 
// ---------------------------------------------------------------------------
void __fastcall Tricky_form::tmr_1Timer( TObject * Sender )
{
    a = b - ( b - a ) * f( b ) / ( f( b ) - f( a ) );
    b = a + ( a - b ) * f( a ) / ( f( a ) - f( b ) );
    Refresh( );
 
    if ( fabs( b - a ) < eps )
    {
        tmr_1->Enabled = false;
        ShowMessage( "a = " + FloatToStr( a ) + ", b = " + FloatToStr( b ) +
            "\nf(a) = " + FloatToStr( f( a ) ) + ", f(b) = " + FloatToStr( f( b ) ) );
    }
}

    int count( 1 );
 
    do
    {
        a = b - ( b - a ) * f( b ) / ( f( b ) - f( a ) );
        b = a + ( a - b ) * f( a ) / ( f( a ) - f( b ) );
 
        if ( fabs( b - a ) < eps )
        {
            ShowMessage( "a = " + FloatToStr( a ) + ", b = " + FloatToStr( b ) +
                "\nf(a) = " + FloatToStr( f( a ) ) + ", f(b) = " + FloatToStr( f( b ) ) );
            break;
        }
    }
    while ( ++count <= 1000 );
 
    if ( count > 1000 )
    {
        ShowMessage( "За 1000 итераций корней не найдено." );
    }