Dear All

I trying to solve the following in MAPLE, the code works fine in MATLAB but I am not sure why it doesn't give me values in MAPLE

n:=3: alpha_p:=1: p:=2: mu0:=4*Pi*1e-7: Br:=1.12:

A1:=unapply(sin((n*p+1)*alpha_p*Pi/(2*p))/((n*p+1)*alpha_p*Pi/(2*p)),n):

A2:=unapply(sin((n*p-1)*alpha_p*Pi/(2*p))/((n*p-1)*alpha_p*Pi/(2*p)),n):

M1:=unapply((Br/mu0)*alpha_p*(A1(n)+A2(n)),n):

M2:=unapply((Br/mu0)*alpha_p*(A1(n)-A2(n)),n):

M3:=unapply(M1(n)+n*p*M2(n),n):

A3:=unapply(((n*p-1/(n*p))*M1(n)/M3(n)+1/(n*p)),n):

When I solve this in Matlab for n=1,2,3 all the values (A1(n), A2(n), M1(n), M2(n), M3(n)) agrees with MAPLE except A3(n) value.

A3(n) value for n=1,2,3 is as follows Pls help me in finding error

Matlab Maple

n=1 6.934462022512822e+15 1.12e+10

n=2 -5.551115123125783e-17 0

n=3 9.245949363350426e+15 -Float(infinity)

Looking forward to your kind replies.

Best Regards

A.Q. (Soton)