Метод прогонки

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, TeEngine, Series, ExtCtrls, TeeProcs, Chart, ComCtrls;
 
type
  TForm1 = class(TForm)
    Chart1: TChart;
    Series1: TFastLineSeries;
    Panel1: TPanel;
    TrackBar1: TTrackBar;
    Label1: TLabel;
    Label2: TLabel;
    Label3: TLabel;
    Label4: TLabel;
    procedure FormCreate(Sender: TObject);
    procedure TrackBar1Change(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
 
implementation
 
{$R *.dfm}
 
const
  N=20;
  M=1000;
  L=1.;
  TT=0.1;
  H=L/N;
  TAU=TT/M;
var TETA: array [0..N, 0..M] of Extended;
function F(TETA: Extended): Extended;
const
  K0=500000000.;
  E=126000.;
  R=8.34;
  Tz=1480.;
  T0=315.;
begin
  F:=K0*(1-TETA)*EXP(-E/(R*(T0+TETA*(Tz-T0))));
end;
 
procedure Calc;
var
  I,J:integer;
  DTETA: Extended;
begin
  for I:=1 to N do TETA[I,0]:=0;
  for J:=0 to M do TETA[0,J]:=1;
  for J:=0 to M do TETA[N,J]:=0;
 
  for J:=0 to M-1 do
  for I:=1 to N-1 do
  begin
    DTETA:=
      +(TETA[I+1,J]-2*TETA[I,J]+TETA[I-1,J])/H/H
      -(TETA[I+1,J]-TETA[I-1,J])/(2*H)
      +F(TETA[I,J]);
    TETA[I,J+1]:=TETA[I,J]+DTETA*TAU;
  end;
end;
 
procedure TForm1.TrackBar1Change(Sender: TObject);
var I, J: integer;
begin
  J:=TrackBar1.Position;
  Series1.Clear;
  for I:=0 to N do
    Series1.AddXY(I*H, TETA[I, J]);
  Label2.Caption:=IntToStr(J);
  Label4.Caption:=FloatToStr(J*TAU);
end;
procedure TForm1.FormCreate(Sender: TObject);
begin
  Calc;
  TrackBar1.Max:=M;
end;
end.