графика в Делфи.

[pionerka]

Здравствуйте,ребята. Я зашла в тупик, 2 месяца уже решаю эту задачу.Силы на исходе.Пожалуйста, помогите как сможите. Иначе остаюсь на осень.
Надо решить на делфи задачу по химии. Все формулы есть уже в коде, я не могу понять только одно.Почему при данных значениях(см. рис) задача не работает(температуру нужно вводить обязательно в Кельвинах).
Сама задача.
Составить модель реакции разложения первого порядка типа А->В+С, протекающей в абиатическом режиме с тепловым эффектом, если известна начальная концентрация и температура.Константа скорости реакции от температуры приведены.
t 50 55 60 65 70 72.5 78 84
к,ч 1,4 1,6 2,2 3,0 4,0 4,7 7,2 9,4
Главное что бы графики строился, когда 1.Cb и Cc равны 0 и 2.когда не равны 0.
Безымянный.JPG

procedure TForm1.Button1Click(Sender: TObject);
//const k:array [1..8] of real=(1.4,1.6,2.2,3.0,4.0,4.7,7.2,9.4);
Const
  dt=0.3;
  l=0.3;
  S=45;
  V=8;
  H=0.4;
Var
  k:real;
  Ca:real;
  Cb:real;
  Cc:real;
  W:real;
  t,Temp:real;
  Tx:real;
begin
  Series1.Clear;
  Series2.Clear;
  Series3.Clear;
  Series4.Clear;
  k:=StrToFLoat(Edit1.Text);
  Ca:=StrToFLoat(Edit2.Text);
  Cb:=StrToFLoat(Edit3.Text);
  Temp:=StrToFLoat(Edit4.Text);
  Cc:=StrToFLoat(Edit5.Text);

  t:=0;
  While Ca>0 do
  begin
      Tx:=Temp+k*Temp*V*H*Ca/(l*S);
      W:=V*H*(Ca/dt)+l*S*(Tx-Temp);
      Ca:=Ca-k*Ca*Cb*t;
      Cb:=Cb+k*Ca*Cb*t;
      Cc:=Cc+k*Ca*Cb*t;
      t:=t+dt;
      series1.AddXY(t,Ca,'',clred);
      series2.AddXY(t,Cb,'',clgreen);
      series3.AddXY(t,Cc,'',clblack);
      series4.AddXY(t,W,'',clblue);
  end;
end;

end.

сама задача в арх
sos.rar

[VVM]

У вас при значениях как на картинке получается бесконечный цикл. Эта формула (при этих значениях): Ca:=Ca-k*Ca*Cb*t; не изменяет значение переменной Са => бесконечный цикл.

[pionerka]

значит надо сделать цикл не бесконечным?

[RUSt88]

нужно пересмотреть формулы

[was3110]

математическая модель процесса с ошибкой. Проверьте формулы. Или попробуйте посчитать вручную. Cb=0 выступает как сомножитель второго слагаемого, и поэтому не будут меняться ни Са, ни Сс, ни Tx, ни W.
W=1983 на протяжении всего цикла и поэтому в график не попадает (>900). Tx=479 на протяжении всего цикла (графика нет). Са..Сс остаются равны начальным значениям. Чтож вы хотите от компьютера? Я не химик, а математик. Если объясните модель (все переменные что значат) разберусь и поправлю. Сайт указан в профиле.

[pionerka]

хорошо.попробую разобраться