Maple. Система диф. уравнений. Вложенные for-циклы.

[Mokusko]

Здравствуйте.
Нужно решить систему диф. уравнений и построить график зависимостей Ir & Is от \Theta & z, но вот никак пока не получается. Решить пытаюсь в Maple.
Система такая:
$\frac{dI_R(z)}{dz}$=-\sum_{i=1}^{M}(\Gamma(0,\Theta_i)\cdot\frac{I_R(z)\cdot(I_S)_i(z)}{I})-\alpha\cdot I_R(z)-\sum_{i=1}^{M}f(\Theta_i)\cdot I_R(z)

$\frac{d(I_S)_i(z)}{dz}$=\sum_{j=0,i \ne j}^{M}(\Gamma(\Theta_i,\Theta_j)\cdot \frac{(I_S(z))_j\cdot (I_S(z))_i}{I})-\alpha\cdot (I_S)_i+f(\Theta_i)\cdot I_R(z)

I=I_R+\sum_{i=1}^{M}(I_S)_i

I_R(z=0)=1, I_S(z=0)=0

Зависимость Г от углов известна. М=50 или любое другое число. i,j – целые. f(theta)=const для начала.
Я пока решала упрощенный вариант, где и Г константа. Пробовала делать через sum, но ничего не получалось. Перешла на двойной цикл for по і и j. Но в результате оно все равно делает что-то не то и, соответственно, ничего не считает (численным методом).
Это имеющийся кусок кода. Здесь I_R(z), I_S[i](z) - x(z), y[i](z)

\alpha := .3; G := 3.5; L := 6; f := 1.1;
for i from 1 by 1 to 50 do
I_0 := x(z)+y[i](z); ICon := x(0) = 1, y[i](0) = 0;
for j from 1 by 1 to 50 do
if i \ne j then
d1 := diff(x(z), z) = -G\cdot x(z)\cdot y[i](z)/I_0-\alpha\cdot x(z);
d2 := diff(y[i](z), z) = G\cdot y[i](z)\cdot y[j](z)/I_0-\alpha\cdot y[i](z);
dsys := {d1, d2};
$F := dsolve({ICon, op(dsys)}, [x(z), y[i](z)], numeric);$
end if;
end do;
end do;

Maple выдает ошибку Error, (in dsolve/numeric/process_input) unknown y[2] present in ODE system is not a specified dependent variable or evaluatable procedure, но я не понимаю, почему.

Можете что-нибудь посоветовать?

[Deggial]

Тема перемещена из форума «Помогите решить / разобраться (М)» в форум «Околонаучный софт»
Причина переноса: в соответствующий раздел.

Mokusko, программный код оформляется тегом code

[one man]

y[j]

А в каком месте Вы задаёте эти значения?