Re: Calculating a simple integral
- To: mathgroup at smc.vnet.net
- Subject: [mg131075] Re: Calculating a simple integral
- From: "Nasser M. Abbasi" <nma at 12000.org>
- Date: Mon, 10 Jun 2013 04:09:03 -0400 (EDT)
- Delivered-to: l-mathgroup@mail-archive0.wolfram.com
- Delivered-to: l-mathgroup@wolfram.com
- Delivered-to: mathgroup-outx@smc.vnet.net
- Delivered-to: mathgroup-newsendx@smc.vnet.net
- References: <kp1e82$80h$1@smc.vnet.net>
- Reply-to: nma at 12000.org
On 6/9/2013 3:26 AM, dsmirnov90 at gmail.com wrote: > If there is a way to calculate with Mathematica the following integral: > > in = -((-1 + Cos[kz])/(kz^2 (kr^2 + kz^2)^2 (kz^2 - 4 \[Pi]^2)^2)) > Integrate[in, {kz, -Infinity, Infinity}, Assumptions -> kr > 0] > > Another system calculates the same integral instantly. :) > > Thanks for any suggestions. > You did not say which version you used. I do get result from 9.01, but it took about 2 minutes -------------------------- ClearAll["`*"]; in = -((-1 + Cos[kz])/(kz^2 (kr^2 + kz^2)^2 (kz^2 - 4 \[Pi]^2)^2)); Assuming[kr > 0, Integrate[in, {kz, -Infinity, Infinity}]] -(1/(32 kr^5 \[Pi]^3 (kr^2+4 \[Pi]^2)^3))(-3 kr^7- 28 kr^5 \[Pi]^2+112 kr^2 \[Pi]^4+192 \[Pi]^6+4 kr \[Pi]^(5/2) (5 kr^2+4 \[Pi]^2) MeijerG[{{1/2,1},{}},{{-(1/2),1/2,1},{0}}, -((I kr)/2),1/2]+4 kr \[Pi]^(5/2) (5 kr^2+4 \[Pi]^2) MeijerG[{{1/2,1},{}},{{-(1/2),1/2,1},{0}},(I kr)/2,1/2]+ 8 kr^3 \[Pi]^(5/2) MeijerG[{{1/2,1},{}},{{-(1/2),1,3/2},{0}}, -((I kr)/2),1/2]+32 kr \[Pi]^(9/2) MeijerG[{{1/2,1},{}}, {{-(1/2),1,3/2},{0}},-((I kr)/2),1/2]+8 kr^3 \[Pi]^(5/2) MeijerG[{{1/2,1},{}},{{-(1/2),1,3/2},{0}},(I kr)/2,1/2]+ 32 kr \[Pi]^(9/2) MeijerG[{{1/2,1},{}},{{-(1/2),1,3/2},{0}},(I kr)/2,1/2]) -------------------------------------- But something does not seem to be right. Can't evaluate any of the MeijerG functions in the result for any kr value. For example N[MeijerG[{{1/2, 1}, {}}, {{-(1/2), 1/2, 1}, {0}}, -((I kr)/2), 1/ 2] /. kr -> 0] gives error MeijerG::hdiv: MeijerG[{{0.5,1.},{}},{{-0.5,0.5,1.},{0.}},0.,0.5] does not exist. Arguments are not consistent. >> --Nasser