Метод Рунге-Кутта 4-го порядка.

[tsx1]

Добрый вечер. Столкнулся с проблемой решения системы уравнений.
Скачал программку, которая считает систему диф.уравнений, пытался изменить на свою, выдал ошибку.

Runge.rar

Необходимо изменить на

RP[0] := ((1.0-k1)*a*X[1]*Xq*(N-X[1]*Xq))/((b+Er+Es)-(1-k1+k2)*(Xq+Xr+Xs)-Xss-Eq-Er-Es);
  //means dx0/dt = ((1-k1)a*x1*Xq*(N-x1*Xq))/((b+ Er+Es)-(1-k1+k2)
  //*(Xq+Xr+Xs)-Xss-Eq-Er-Es)                                                               //
  RP[1] := Eq-e1*X[1];
  //means dx1/dt = Eq-e1*x1
  RP[2] := Er-e2*X[2];
  //means dx2/dt = Er-e2*x2
  RP[3] := Es-e3*X[3];
  //means dx3/dt = Es-e3*x3

Это система уравнений

  k1:=0.2;
    k2:=0.34;
    a:=3.2*(1/10000000);
    b:= 0.04;
    N:=100000;
    e1:=0.02;
    e2:=0.005;
    e3:=0.005;
    randomize ;
    Xq:=random(10000);
    Xr:=random(10000);
    Xs:=random(10000);
    Xss:=random(10000);
    Er:=random(10000);
    Eq:=random(10000);
    Es:=random(10000);

Коэффициенты.