Advanced Search

Question: Facing issue in integration evaluations

Question:Facing issue in integration evaluations

Tamour_Zubair 80 Maple 2015
integration numeric
February 15 2022
2

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
Please Wait...