[Thiago]

Собственно построить биффуркационную диаграмму(т.е. зависимость параметра l(так обозначил лямбду) от х). Сам дельфи осваивать начал недавно, программа получилась кривая, так что не судите строго. Решил ур. ВдП методом рунге кутта, дальше насколько понял - надо сделать 3 цикла - по l,в нем по времени, в этом еще по иксу. В программе еще 2 графика - фазовый портрет и временная реализация, но главное построить этот. Помогите, если кто может

Код:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136

unit Unit1;    interface    uses   Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,   Dialogs, TeEngine, Series, ExtCtrls, TeeProcs, Chart, StdCtrls, Buttons;    type   TForm1 = class(TForm)     LabeledEdit1: TLabeledEdit;     LabeledEdit2: TLabeledEdit;     LabeledEdit3: TLabeledEdit;     Chart1: TChart;     Chart2: TChart;     Chart3: TChart;     Series1: TPointSeries;     Series2: TPointSeries;     Series3: TPointSeries;     BitBtn1: TBitBtn;     BitBtn2: TBitBtn;     procedure BitBtn1Click(Sender: TObject);     procedure BitBtn2Click(Sender: TObject);     procedure Chart1DblClick(Sender: TObject);   private     { Private declarations }   public     { Public declarations }   end;    var   Form1: TForm1;   l,w:double;    implementation    {$R *.dfm}       function dx(x,y:double):double; begin result:=y; end;    function dy(x,y:double):double; begin result:=(l-x*x)*y-w*w*x; end;    procedure runge(var x,y:double); const dh=0.05; var k1,k2,k3,k4, l1,l2,l3,l4, q1,q2,q3,q4:double; begin k1:=dx(x,y); l1:=dy(x,y);       k2:=dx(x+dh*k1/2,y+dh*q1/2); l2:=dy(x+dh*k1/2,y+dh*q1/2);       k3:=dx(x+dh*k2/2,y+dh*q2/2); l3:=dy(x+dh*k2/2,y+dh*q2/2);       k4:=dx(x+dh*k3,y+dh*q3); l4:=dy(x+dh*k3,y+dh*q3);       x:=x+dh*(k1+2*k2+2*k3+k4)/6; y:=y+dh*(l1+2*l2+2*l3+l4)/6;    end;    procedure TForm1.BitBtn1Click(Sender: TObject); var t_stop,time,x,y:double; begin Chart1.Series[0].Clear; l:=StrToFloat(LabeledEdit1.Text); w:=StrToFloat(LabeledEdit2.Text);    t_stop:=30000;     While l<5 do   begin   x:=0.1;   y:=0.1;      time:=0;     While time<t_stop do       begin       runge(x,y);       if time>25000 then       Chart1.Series[0].AddXY(l,x);       time:=time+1;       //Application.ProcessMessages;       end;       //Chart1.Series[0].AddXY(l,x);   l:=l+0.1;  end;    end;       procedure TForm1.BitBtn2Click(Sender: TObject); begin halt; end;    procedure TForm1.Chart1DblClick(Sender: TObject); var x,z,y,t_stop1,time1: double; begin Chart2.Series[0].Clear; Chart3.Series[0].Clear; Chart1.Series[0].GetCursorValues(l,x); z:=0.1; y:=0.1; t_stop1:=2000; time1:=0.1;  while time1<t_stop1 do   begin   runge(x,y);   Chart2.Series[0].AddXY(y,z);   Chart3.Series[0].AddXY(time1,x);   Application.ProcessMessages;   time1:=time1+0.1; end; end;       end.