Question: Matlab to Maple

November 03 2012 AliKhan 10
Maple 15

1

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)

Please Wait...