Deferensial Heun

O... ya, ne da lagi yang penyelesaian PDB nya pake metode Heun, tapi yang ne menurut versi saya.....


var Form1: TForm1;

implementation

{$R *.dfm}


procedure TForm1.Button1Click(Sender: TObject);
var X,Y,H,AY,BY,CY,DY:Real;
I:Integer;

begin X:=StrToFloat(Edit1.Text);
Y:=StrToFloat(Edit2.Text);
H:=StrToFloat(Edit3.Text);
I:=1;
DY:=Y;
repeat
Y:=DY;
//----------------slope titik awal
AY:=4*exp(0.8*X)-0.5*Y;
//----------------pers prediktor
BY:=Y+AY*H;
//----------------slope titik akhir
CY:=4*Exp(0.8*(X+H))-0.5*BY;
//----------------pers korektor
DY:=Y+(AY+CY)*H/2;
//----------------tampilkan
StringGrid1.Cells[0,I]:=FloatToStr(X);
StringGrid1.Cells[1,I]:=FloatToStr(Y);
StringGrid1.Cells[2,I]:=FloatToStr(AY);
StringGrid1.Cells[3,I]:=FloatToStr(BY);
StringGrid1.Cells[4,I]:=FloatToStr(CY);
StringGrid1.Cells[5,I]:=FloatToStr(DY);
X:=X+H;
Inc(I);
until abs((DY-Y)/DY)<=0.01
end;

end.

enjoy the code..........