Метод Рунге-Кутта и многое другое

#include <vcl.h>
#include <iostream.h>
#include <conio.h>
#include <math.h>
#pragma hdrstop
 
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
 double f(double u, double a, double s, double x, double w, double d)
{
  return a*cos(s)+2*s*u+w*w*x;
}
 double f1(double u, double v)
{
  return u;
}
void __fastcall TForm1::Button1Click(TObject *Sender)
{
double a = Edit1->Text.ToDouble();
double u = Edit2->Text.ToDouble();
double s = Edit3->Text.ToDouble();
double x = Edit4->Text.ToDouble();
double w = Edit5->Text.ToDouble();
double d = Edit6->Text.ToDouble();
double v = Edit4->Text.ToDouble();
int N = Edit7->Text.ToDouble(), i;
double h = 0.1, k1, k2, k3, k4, m1, m2, m3, m4;
for(i = 0; i<N-1; i++)
   {
                   k1 = h*f(a, u, s, x, w, d);
                   k2 = h*f(u + h/2, x + (k1)/ 2, a, w, s, d);
                   k3 = h*f(u + h/2, x + (k2)/ 2, a, w, s, d);
                   k4 = h*f(u + h, x + k3, a, w, s, d);
                   x+= (k1 + 2*k2 + 2*k3 + k4)/6;
                   u+= h;
                   Series1->AddXY(u,x);
   }
      Edit8->Text = x;
for(i = 0; i<N-1; i++)
{
                   m1 = h*f1(u,v);
                   m2 = h*f1(u + h/2, v + (k1)/ 2);
                   m3 = h*f1(u + h/2, v + (k2)/ 2);
                   m4 = h*f1(u + h, v + k3);
                   v+= (k1 + 2*k2 + 2*k3 + k4)/6;
                   u+= h;
                   Series2->AddXY(u,v);
}
      Edit9->Text = v;
}

//---------------------------------------------------------------------------
 
#include <vcl.h>
#include <iostream.h>
#include <conio.h>
#include <math.h>
#pragma hdrstop
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
double f(double u, double a, double s, double x, double w, double d)
{
  return a*cos(s)+2*s*u+w*w*x;
}
 double f1(double u, double v)
{
  return u;
}
 
void __fastcall TForm1::Button1Click(TObject *Sender)
{
double a = Edit1->Text.ToDouble();
double u = Edit2->Text.ToDouble();
double s = Edit3->Text.ToDouble();
double x = Edit4->Text.ToDouble();
double w = Edit5->Text.ToDouble();
double d = Edit6->Text.ToDouble();
double v = Edit4->Text.ToDouble();
int N = Edit7->Text.ToDouble(), i, E=5;
double k1, k2, k3, k4, m1, m2, m3, m4;
double h = 0.1;
for(i = 0; i<N-1; i++)
   {
                   k1 = h*f(a, u, s, x, w, d);
                   k2 = h*f(u + h/2, x + (k1)/ 2, a, w, s, d);
                   k3 = h*f(u + h/2, x + (k2)/ 2, a, w, s, d);
                   k4 = h*f(u + h, x + k3, a, w, s, d);
                   x+= (k1 + 2*k2 + 2*k3 + k4)/6;
                   u+= h;
                   m1 = h*f1(u,v);
                   m2 = h*f1(u + h/2, v + (k1)/ 2);
                   m3 = h*f1(u + h/2, v + (k2)/ 2);
                   m4 = h*f1(u + h, v + k3);
                   v+= (m1 + 2*m2 + 2*m3 + m4)/6;
                   u+= h;
                   double T = f(u,a,s,x,w,d);
                   Series1->AddXY(T,x+v);
   }
      Edit8->Text=printf("%1.2f \n", x);
      Edit9->Text = v;
}