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)