Здравствуйте, у меня возник вопрос в процессе решения уравнения в мапле, в уравнении имеются две переменных для которых собственно и требуется создани цикла чтобы программа автоматически решила мне уравнение:

restart;
> mu_in:=2.5; mu_term:=7.5; A0:=0.15; DA:=0.15; N_amp:=5;

mu_in := 2.5


mu_term := 7.5


A0 := 0.15


DA := 0.15


N_amp := 5

> N_mu:=4;

N_mu := 4

> D_mu:=(mu_term-mu_in)/N_mu;

D_mu := 1.250000000

> for i to N_mu do
> mu:=mu_in+D_mu*(i-1):
> eq[i]:=((mu+1)*(M*(M+mu)-A^2))/(((mu+M)^2-A^2)*sqrt(M^2-A^2))=1;
> for j to N_amp do
> A:=A0+DA*(j-1);
> M[i,j]:=fsolve(eq[i],M);
> end do;
> end do;

программа выдает лишь само уравнение с подставленными неизвестными не могу понять где ошибся