Интегрирование системы уравнений Гамильтона для модели Хенона-Хейлеса методом Рунге-Кутта 4-го порядка

#include <math.h>
#include <iostream.h>
#include "runge_kutta.h"
//---------------------------------------------------------------------------
#pragma hdrstop
//---------------------------------------------------------------------------
#pragma argsused
double m = 1.0;
void rhs(double *y,double *dy,double t)
{
    dy[0] = y[2];
    dy[1] = y[3];
    dy[2] = -(y[0]+2*y[0]*y[1]);
    dy[3] = -(y[1]+y[0]*y[0]-y[1]*y[1]);
}
double EDS(double t)
{
return 0;
}
 
int main(int argc, char* argv[])
{
    int n = 4;
    double ht = 0.05;
    RungeKutta *rk = new RungeKutta(n,ht,rhs);
 
    double
        Etot=1/8,
        V=(1/12),
        w=(Etot-V)*2*m,
        *in = new double [n];
    in[0] = 0.1;
    in[1] = 0.2;
    in[2] = m*sqrt(w);
    in[3] = 0.3;
 
    rk->setInitialValues(in);
 
    double t = 0.1;
    for(int it=0;it<100;it++)
    {
        rk->doOneRKStep(t);
        t += ht;
 
        cout << "x=" << rk->y1[1] << "  y="<< rk->y1[3] << "\n";
    }
    return 1;
}
//---------------------------------------------------------------------------

//---------------------------------------------------------------------------
#include <vcl.h>
#pragma hdrstop
 
#include "Unit4.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm4 *Form4;
 
//---------------------------------------------------------------------------
__fastcall TForm4::TForm4(TComponent* Owner)
        : TForm(Owner)
{
    DrawingBoard=new Graphics::TBitmap;
    DrawingBoard->Width = this->PaintBox1->Width;
    DrawingBoard->Height = this->PaintBox1->Height;
}
//---------------------------------------------------------------------------
void __fastcall TForm4::PaintBox1MouseDown(TObject *Sender,
      TMouseButton Button, TShiftState Shift, int X, int Y)
{
    IsDrawing = True;
    StartX = X;
    StartY = Y;
}
//---------------------------------------------------------------------------
void __fastcall TForm4::PaintBox1MouseUp(TObject *Sender,
      TMouseButton Button, TShiftState Shift, int X, int Y)
{
    if( IsDrawing )
    {
        DrawingBoard->Canvas->Rectangle(StartX, StartY, X, Y);
        PaintBox1->Canvas->Draw(0, 0, DrawingBoard);
        IsDrawing = False;
    }
}
//---------------------------------------------------------------------------
void __fastcall TForm4::PaintBox1Paint(TObject *Sender)
{
    PaintBox1->Canvas->Draw(0, 0, DrawingBoard);
}
//---------------------------------------------------------------------------
 
 
void __fastcall TForm4::Button1Click(TObject *Sender)
{
}
 
//---------------------------------------------------------------------------