Решение диф. уравнений 2-го порядка методом эйлера

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, AxCtrls, OleCtrls, VCF1, Math, TeEngine, Series,
  ExtCtrls, TeeProcs, Chart, XPMan, OleCtnrs;
 
type
  TForm1 = class(TForm)
    Button1: TButton;
    F1Book1: TF1Book;
    Edit1: TEdit;
    Edit3: TEdit;
    Button2: TButton;
    Chart1: TChart;
    Series1: TLineSeries;
    Series2: TLineSeries;
    XPManifest1: TXPManifest;
    Label2: TLabel;
    Label4: TLabel;
    Label6: TLabel;
    procedure Button1Click(Sender: TObject);
    procedure Button2Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
 
implementation
 
{$R *.dfm}
 
function f1(x:real):real;
begin
f1:=-2*sqr(Cos(x))+Cos(x);      //первое уравнение
end;
 
function f2(x:real;y:real):real;
begin                              //второе уравнение
f2:=Sin(2*x)-y*(Tan(x));
end;
 
procedure TForm1.Button1Click(Sender: TObject);
type mass=array [1..20] of real;
var a,h:real;
xe,ye:mass;         //объявление переменных
i:integer;
begin
 
 
 
begin
a:=strtofloat(Edit1.Text);  //ввод начального значения и шага
h:=strtofloat(Edit3.Text);
 
xe[1]:=a;
ye[1]:=f1(xe[1]);
 
for i:=1 to 11 do
begin
ye[i+1]:=ye[i]+h-2*ye[i]); //алгоритм Эйлера для уравнения 1-ого порядка
xe[i+1]:=xe[i]+h*f1;
f1book1.NumberRC[i,1]:=ye[i]; //вывод в F1book
end;
 
 end;
 
end;
 
procedure TForm1.Button2Click(Sender: TObject);
type mass=array [1..20] of real;
var a,h:real;
x:mass;
 xe,ye:mass;
   i:integer;
begin
a:=strtofloat(Edit1.Text);
h:=strtofloat(Edit3.Text);
x[1]:=a;
for i:=1 to 11 do
begin
Chart1.Series[0].AddXY(x[i],f1(x[i])) ; //рисование графика
x[i+1]:=x[i]+h;
end;
 
xe[1]:=a;
ye[1]:=f1(xe[1]);
 
for i:=1 to 11 do
begin
ye[i+1]:=ye[i]+h*f2(xe[i],ye[i]);
xe[i+1]:=xe[i]+h;
 
Chart1.Series[1].AddXY(xe[i],ye[i]) ;  //рисование графика
end;
end;
 
end.