Tamour_Zubair

80 Reputation

One Badge

2 years, 61 days

MaplePrimes Activity


These are questions asked by Tamour_Zubair

I want to apply rk-2 for the following system of odes. I want to know what is process of rk-2 method maple used? Is it possible that we can see the complete process maple used?

Download Question.mw

 

Hello Sir
Hope you are fine. I am facing issues in integration evaluations. I have attched code and pic as well. Kindly guide me. I am waiting for your kind response.
Thanks

question.mw

restart; printlevel := 4; kernelopts(version)

{--> enter Terminate, args =

<-- exit Terminate (now at top level) = }
{--> enter ModuleUnload, args =
<-- exit ModuleUnload (now at top level) = }
{--> enter OnUnload, args =
<-- exit OnUnload (now at top level) = }

 

`Maple 2015.0, X86 64 WINDOWS, Feb 17 2015, Build ID 1022128`

(1)

with(LinearAlgebra):

L[time] := 1:

ff11 := cos(1.0*k*(p+q+x+t)):

NULL

JJx11 := QQ*(int(int(p*ff11/`&gamma;&gamma;Exact`, p), q))/m;

175631174533.479692645444566411*(int(int(0.273299999999999999999999999999e-21*p*cos(.50*p+.50*q+.50*x+.50*t)/(p^2+q^2+0.746928900000000000000000000000e-43)^(1/2), p), q))

(2)

evalf(IntegrationTools:-Expand(%));

Warning,  computation interrupted

 

``

``

int(int(p*ff11/`&gamma;&gamma;Exact`, p = 0 .. 1), q = 0 .. 1); evalf(IntegrationTools:-Expand(%))

int(int(0.273299999999999999999999999999e-21*p*cos(.50*p+.50*q+.50*x+.50*t)/(p^2+q^2+0.746928900000000000000000000000e-43)^(1/2), p = 0 .. 1), q = 0 .. 1)

 

0.273299999999999999999999999999e-21*(Int(Int(p*cos(.50*p)*cos(.50*q)*cos(.50*x)*cos(.50*t)/(p^2+q^2+0.746928900000000000000000000000e-43)^(1/2)-1.*p*cos(.50*p)*cos(.50*q)*sin(.50*x)*sin(.50*t)/(p^2+q^2+0.746928900000000000000000000000e-43)^(1/2)-1.*p*cos(.50*p)*sin(.50*q)*sin(.50*x)*cos(.50*t)/(p^2+q^2+0.746928900000000000000000000000e-43)^(1/2)-1.*p*cos(.50*p)*sin(.50*q)*cos(.50*x)*sin(.50*t)/(p^2+q^2+0.746928900000000000000000000000e-43)^(1/2)-1.*p*sin(.50*p)*sin(.50*q)*cos(.50*x)*cos(.50*t)/(p^2+q^2+0.746928900000000000000000000000e-43)^(1/2)+p*sin(.50*p)*sin(.50*q)*sin(.50*x)*sin(.50*t)/(p^2+q^2+0.746928900000000000000000000000e-43)^(1/2)-1.*p*sin(.50*p)*cos(.50*q)*sin(.50*x)*cos(.50*t)/(p^2+q^2+0.746928900000000000000000000000e-43)^(1/2)-1.*p*sin(.50*p)*cos(.50*q)*cos(.50*x)*sin(.50*t)/(p^2+q^2+0.746928900000000000000000000000e-43)^(1/2), p = 0. .. 1.), q = 0. .. 1.))

(3)

``

Download question.mw

1 2 3 4 Page 4 of 4