Моделирование в Делфи

[Бриит]

unit Unit1;

interface

uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, TeEngine, Series, TeeProcs, Chart, StdCtrls, ExtCtrls;

type
  TForm1 = class(TForm)
    Panel1: TPanel;
    Button1: TButton;
    Chart1: TChart;
    Label1: TLabel;
    Label2: TLabel;
    Edit1: TEdit;
    Edit2: TEdit;
    Button2: TButton;
    Label3: TLabel;
    Label4: TLabel;
    Edit3: TEdit;
    Edit4: TEdit;
    Label5: TLabel;
    Label6: TLabel;
    Edit5: TEdit;
    Edit6: TEdit;
    Label7: TLabel;
    Label8: TLabel;
    Edit7: TEdit;
    Edit8: TEdit;
    Label10: TLabel;
    Edit9: TEdit;
    Series2: TPointSeries;
    Series1: TPointSeries;
    Timer1: TTimer;
    Edit10: TEdit;
    procedure Button1Click(Sender: TObject);
    procedure Button2Click(Sender: TObject);
    procedure Timer1Timer(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  Form1: TForm1;
  N10,N20:real;
  N1,N2:array [0..1000] of real;
  t:real;
  r1,r2,k1,k2,a12,a21,dt:real;
  i:integer;
implementation

{$R *.dfm}

procedure TForm1.Button1Click(Sender: TObject);
{const
  r1=2; r2=2;
  k1=200; k2=200;
  a12=0.5; a21=0.5;
  dt=0.1; }

begin
  N10:= Random(100)+5;
  N20:= Random(150)+10;
   r1:=strtofloat(edit3.text);
   r2:=strtofloat(edit4.text);
   k1:=strtofloat(edit5.text);
   k2:=strtofloat(edit6.text);
   a12:=strtofloat(edit7.text);
   a21:=strtofloat(edit8.text);
   dt:=strtofloat(edit9.text);

  N1[0]:=N10; N2[0]:=N20;
  t:=0;   i:=0;
  Timer1.Enabled:=true;
  //hile (abs(n1[i]-n1[i-1])>2) or (abs(n2[i]-n2[i-1])>2) do begin
 {   repeat
    series1.AddXY(t,N1[i],'',clred);
    series2.AddXY(t,N2[i],'',clgreen);
    inc(i);
    t:=t+dt;
    N1[i]:=r1*N1[i-1]*(k1-N1[i-1]-a12*N2[i-1])/k1;
    N2[i]:=r2*N2[i-1]*(k2-N2[i-1]-a21*N1[i-1])/k2;
    Edit10.Text:= FloatToStr(t);
  until {(abs(n1[i]-n1[i-1])<0.5) or (abs(n2[i]-n2[i-1])<0.5)  t>10  }
end;

procedure TForm1.Button2Click(Sender: TObject);
begin
series1.Clear;
series2.Clear;
end;

procedure TForm1.Timer1Timer(Sender: TObject);
begin
    series1.AddXY(t,N1[i],'',clred);
    series2.AddXY(t,N2[i],'',clgreen);
    inc(i);
    t:=t+dt;
    N1[i]:=r1*N1[i-1]*(k1-N1[i-1]-a12*N2[i-1])/k1;
    N2[i]:=r2*N2[i-1]*(k2-N2[i-1]-a21*N1[i-1])/k2;
    Edit10.Text:= FloatToStr(t);
    Timer1.Interval:=1000;
    if t>10 then  Timer1.Enabled:=False;

end;

end.

вот и все так кажись

[Бриит]

это прикольно бегающие точки