метод рунге-кутты 3-го порядка в билдер

//---------------------------------------------------------------------------
 
#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)
{
}
//---------------------------------------------------------------------------
 
 
 
 float f1(float x, float p)
{
      return(p);
}
 
float f2(float x, float y, float p)
{
      return(pow(2.7,-x));
}
float x0=1;
float X=1.3;
float h=0.1;
float y0=0.77;
float p0=0.44;
float N=10;
void __fastcall TForm1::Button1Click(TObject *Sender)
{   for (int i=1;i<N;i++)
    {  float p,x,y;
    float k1=f1( x, p);
    float m1=f2(x,y,p);
   /*  StringGrid1->Cells[i][1]=k1;
     StringGrid1->Cells[i][2]=m1; */
       float k2=f1( x+h/2, p+m1*h/2);
    float m2=f2(x+h/2,y+k1*h/2,p+m1*h/2);
      float k3=f1( x+h, p+2*m2-m1);
    float m3=f2(x+h,y+2*k2-k1,p+2*m2-m1);
    x=x+h;
    y=y+(k1+2*k2+2*k3)*h/6;
    p=p+(m1+2*m2+2*m3)*h/6;
    StringGrid1->Cells[i][1]=x;
     StringGrid1->Cells[i][2]=y;
       StringGrid1->Cells[i][3]=p;
    }
 
}

//---------------------------------------------------------------------------
 
#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)
{     StringGrid1->Cells[0][0]="x";
     StringGrid1->Cells[1][0]="y";
       StringGrid1->Cells[2][0]="y'";
}
//---------------------------------------------------------------------------
 
 
 
 float f1(float x, float p)
{
      return(p);
}
 
float f2(float x, float y, float p)
{
      return(pow(exp(1),-x));
}
/*float x0=1;
float X=1.3;
float h=0.1;
float y=0.77;
float p=0.44;
float N=10; */
void __fastcall TForm1::Button1Click(TObject *Sender)
{float p[30],x[30],y[30],h=0.1;
    float N=12;
    x[0]=0;
    y[0]=0;
    p[0]=0;
 for (int i=1;i<N;i++)
    {
    float k1=f1( x[i-1], p[i-1]);
    float m1=f2(x[i-1],y[i-1],p[i-1]);
  /*   StringGrid1->Cells[i][1]=k1;
     StringGrid1->Cells[i][2]=m1;    */
       float k2=f1( x[i-1]+h/2, p[i-1]+m1*h/2);
    float m2=f2(x[i-1]+h/2,y[i-1]+k1*h/2,p[i-1]+m1*h/2);
      float k3=f1( x[i-1]+h, p[i-1]+2*m2-m1);
    float m3=f2(x[i-1]+h,y[i-1]+2*k2-k1,p[i-1]+2*m2-m1);
    x[i]=x[i-1]+h;
    y[i]=y[i-1]+(k1+2*k2+2*k3)*h/6;
    p[i]=p[i-1]+(m1+2*m2+2*m3)*h/6;
    StringGrid1->Cells[0][i]=x[i-1];
     StringGrid1->Cells[1][i]=y[i-1];
       StringGrid1->Cells[2][i]=p[i-1];
    }
 
}