Advanced Search

Question: help with numerical answer to maple expression

Question:help with numerical answer to maple expression

tintin 128 Maple
input
May 02 2007
0
Hello, I have the following maple expression (sorry, it is somewhat lengthy)

e_x0 :=

4*ep^2*sin(m*Pi/b*y)*sin(n*Pi/c*z)*J0*(x1-x0)*(-cos(n*Pi/c*z1)+cos(n*Pi/c*z0))*(-cos(m*Pi/b*y1)+cos(m*Pi/b*

y0))/((2*m^2*Pi^4*c^2*n^2*b^2+n^4*Pi^4*b^4+8*Pi^4*m^4*c^4*Q^2+16*Pi^4*m^4*c^4*Q^4+16*k^4*b^4*c^4*Q^4+8*m^2*

Pi^2*c^4*k^2*b^2*Q^2-32*Pi^2*n^2*b^4*Q^4*k^2*c^2+8*n^2*Pi^2*b^4*k^2*c^2*Q^2+m^4*Pi^4*c^4-32*Pi^2*m^2*c^4*Q^

4*k^2*b^2+16*Pi^4*m^2*c^2*Q^2*n^2*b^2+8*Pi^4*n^4*b^4*Q^2+32*Pi^4*m^2*c^2*Q^4*n^2*b^2+16*Pi^4*n^4*b^4*Q^4)/b

^4/c^4/Q^4)^(1/2)/a/m/Pi^2/n;


and the input parameters are

a := 0.01; b := 0.03; c := 0.08;

ep := 2; l := 1; m := 1; n := 1; Q := 1.0e4; k := 5.339e9*2*Pi/3e8;

x0 := 0.002; x1 := 0.004; y0 := 0.002; y1 := 0.004; z0 := 0.002; z1 := 0.004;

x := 0.004; y := 0.006; z := 0.006;

J0 := 200;



Using evalf(e_x0), Maple give the numerical answer as 0.000272511...

However, other programs (Matlab, Mathematica, etc) give answer as 0.000272946...
Does anyone know how the discrepancy come about? Thanks!

-Dan
Please Wait...