ПОЖАЛУЙСТА, Помогите разобраться в готовой программе и найти ошибку

unit TeploObm;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, Spin, ExtCtrls;
 
type
  Res=array [1..501] of real;
 
type
  TForm1 = class(TForm)
    Timer1: TTimer;
    PaintBox1: TPaintBox;
    SpinEdit1: TSpinEdit;
    SpinEdit2: TSpinEdit;
    Label1: TLabel;
    Label2: TLabel;
    Button1: TButton;
    Button2: TButton;
    Label3: TLabel;
    Label4: TLabel;
    Label5: TLabel;
    Button3: TButton;
    Label6: TLabel;
    procedure PaintBox1Paint(Sender: TObject);
    procedure Button1Click(Sender: TObject);
    procedure Button2Click(Sender: TObject);
    //
    function Fu(omega,ti:real):real;
    function Progonka(a,b,c,f: Res):Res;
    procedure MasNul(var a,b,c,f:Res);
    procedure FonPolya;
    procedure Timer1Timer(Sender: TObject);
    procedure Button3Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
const
  dx=0.01; l=5;
  dt=0.1;
  a2=10; // k=, c=, ro=
 
var
  Form1: TForm1;
  i,tx: Integer;
  str: Boolean;
  omega,n: Integer;
  alfa,t,ti,ui:real;
  u0,u,a,b,c,f: Res;
 
implementation
 
uses Unit2;
 
{$R *.dfm}
 
function TForm1.Fu(omega, ti:real):real;
begin
   Fu:=20+10*sin(omega*ti);
end;
 
function TForm1.Progonka(a,b,c,f:Res):Res;
var gam,bet,pr: Res;
    fi: real;
begin
   bet[1]:=0; gam[1]:=0;
   for i:=1 to n do
   begin
      fi:=b[i];
      if i>1 then fi:=fi-a[i]*bet[i];
      if fi=0 then break;
      fi:=1/fi;
      if i<n then bet[i+1]:=c[i]*fi;
      gam[i+1]:=(f[i]+a[i]*gam[i])*fi;
   end;
   pr[n]:=gam[n+1];
   for i:=n downto 2 do
     pr[i-1]:=bet[i]*pr[i]+gam[i];
   Progonka:=pr;
end;
 
procedure TForm1.FonPolya;
begin
   PaintBox1.Canvas.Rectangle(0,0,650,310);
   PaintBox1.Canvas.MoveTo(30,270);
   PaintBox1.Canvas.LineTo(630,270);
   PaintBox1.Canvas.MoveTo(50,290);
   PaintBox1.Canvas.LineTo(50,20);
   PaintBox1.Canvas.MoveTo(47,25);
   PaintBox1.Canvas.LineTo(50,20);
   PaintBox1.Canvas.LineTo(53,26);
   PaintBox1.Canvas.MoveTo(625,267);
   PaintBox1.Canvas.LineTo(630,270);
   PaintBox1.Canvas.LineTo(624,273);
   PaintBox1.Canvas.TextOut(30,15,'U(t)');
   PaintBox1.Canvas.TextOut(630,270,'t');
end;
 
procedure TForm1.MasNul(var a,b,c,f: Res);
begin
   for i:=1 to 501 do
   begin
      a[i]:=0;
      b[i]:=0;
      c[i]:=0;
      f[i]:=0;
   end;
end;
 
procedure TForm1.PaintBox1Paint(Sender: TObject);
begin
   str:=False;
   FonPolya;
   n:=round(l/dx);
end;
 
procedure TForm1.Button1Click(Sender: TObject);
begin
   str:=not str;
   if str then begin
     Button1.Caption:='Стоп';
     alfa:=StrToInt(SpinEdit1.Text)*0.000001;
     omega:=StrToInt(SpinEdit2.Text);
     t:=5*sqr(alfa/(l*omega)); ti:=0;
     FonPolya;
     for i:=1 to 30 do
     begin
        PaintBox1.Canvas.MoveTo(48,270-i*8);
        PaintBox1.Canvas.LineTo(53,270-i*8);
     end;
     tx:=round(t);
     if frac(t)<0.5 then inc(tx);
     for i:=1 to tx do
     begin
        PaintBox1.Canvas.MoveTo(50+i*round(570/tx),268);
        PaintBox1.Canvas.LineTo(50+i*round(570/tx),273);
     end;
     for i:=1 to n do
       u0[i]:=20;
     PaintBox1.Canvas.MoveTo(50,110);
     PaintBox1.Canvas.Pen.Color:=clRed;
     Timer1.Enabled:=True;
   end
   else begin
     Timer1.Enabled:=False;
     Button1.Caption:='Запуск';
   end;
end;
 
procedure TForm1.Button2Click(Sender: TObject);
begin
   Form1.Close;
end;
 
procedure TForm1.Timer1Timer(Sender: TObject);
begin
   if ti>t then begin
     Timer1.Enabled:=False;
     Button1.Caption:='Запуск';
     str:=not str;
     PaintBox1.Canvas.Pen.Color:=clBlack;
   end else begin
     MasNul(a,b,c,f);
     for i:=1 to n+1 do
     begin
        if i>1 then a[i]:=-(a2*dt/(2*sqr(dx)));
        b[i]:=-(1+a2*dt/sqr(dx));
        c[i]:=-(a2*dt/(2*sqr(dx)));
        if i=1 then
          f[i]:=-(u0[i]+(a2*dt/(2*sqr(dx)))*(u0[i+1]-2*u0[i]+Fu(omega,ti)+Fu(omega,(ti+dt))))
        else f[i]:=-(u0[i]+(a2*dt/(2*sqr(dx)))*(u0[i-1]-2*u0[i]+u0[i+1]));
     end;
     u:=Progonka(a,b,c,f);
     PaintBox1.Canvas.LineTo(50+round(ti*570/tx),270-round(u[250]*240/20));
     u0:=u;
     ti:=ti+dt;
   end;
end;
 
procedure TForm1.Button3Click(Sender: TObject);
begin
Form2.Show;
form2.Image1.Picture.LoadFromFile('задание.bmp');
end;
 
end.