Re: Calculating a simple integral

• To: mathgroup at smc.vnet.net
• Subject: [mg131141] Re: Calculating a simple integral
• From: "Brambilla Roberto Luigi (RSE)" <Roberto.Brambilla at rse-web.it>
• Date: Fri, 14 Jun 2013 04:59:35 -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: <7586978.94875.1370767085206.JavaMail.root@m06> <20130613063855.8617D6A70@smc.vnet.net>

```Using Mathematica 5.1 or 8  I obtained:

in = -((-1 + Cos[z])/(z^2 (r^2 + z^2)^2 (z^2 - 4 =A6=D0^2)^2));

*numerically*
r = 0.1;
NIntegrate[in, {z, -Infinity, Infinity}]
0.50377

*analytically*
f1[x_] := MeijerG[{{1, 1,3/2},{}},{{1, 1, 3/2}, {0, 1/2}}, x, 1/2]
f2[x_] := MeijerG[{{1, 1,3/2},{}},{{1, 3/2, 2}, {0, 1/2}}, x, 1/2]

A = (3r^7+28r^5Pi^2) + 16Pi^(5/2) (5r^2+4Pi^2) f1[ I r/2] + 32Pi^(5/2) (r^2+4Pi^2) f2[ I r/2];

B = 32r^5Pi^3 (r^2+4Pi^2);

Im[A/B]

0.50377

Why to hate Meijer functions with parameter 1/2?

Roberto Brambilla

-----Messaggio originale-----
Da: Dr. Wolfgang Hintze [mailto:weh at snafu.de]
Inviato: gioved=A8=AC 13 giugno 2013 8.39
A: mathgroup at smc.vnet.net
Oggetto: Re: Calculating a simple integral

On 11 Jun., 08:23, Andrzej Kozlowski <akozlow... at gmail.com> wrote:
> No, it's similar to:
>
> Integrate[(1 -
>     Cos[x])/(x^2*(x^2 - 4*Pi^2)^2), {x, -Infinity, Infinity}]
>
> 3/(32*Pi^3)
>
> On 10 Jun 2013, at 10:11, djmpark <djmp... at comcast.net> wrote:
>
>
>
> > Doesn't this have a singularity at 2 Pi that produces
> > non-convergence? It's similar to:
>
> > Integrate[1/x^2, {x, \[Epsilon], \[Infinity]}, Assumptions ->
> > \[Epsilon] > 0]
>
> > 1/\[Epsilon]
>
> > That diverges as epsilon -> 0.
>
> > Are you sure you copied the integral correctly?
>
> > David Park
> > djmp... at comcast.net
> >http://home.comcast.net/~djmpark/index.html
>
> > From: dsmirno... at gmail.com [mailto:dsmirno... at gmail.com]
>
> > 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.

Sorry, but I made indeed a calculation error!
Correcting it the partial fraction decomposition leads to Dmitry's result.
Furthermore, calculating first the indefinite integral and then taking limits leads to a false result.
Direct calculation of the integral leads to MeierG functions which are useless because we cannot enter any numerical value.
So, rather than provding the correct result Mathematica comes up with different false result depending on the method used, and we cannot tel which one is correct without "research" work.
Summarizing, I need to restate my criticism of Mathematica with respect to integration (I'm using version 8).

Regards,
Wolfgang

RSE SpA ha adottato il Modello Organizzativo ai sensi del D.Lgs.231/2001, in=
forza del quale l'assunzione di obbligazioni da parte della Societ=A8=A4 av=
viene con firma di un procuratore, munito di idonei poteri.
RSE adopts a Compliance Programme under the Italian Law (D.Lgs.231/2001). Ac=
cording to this RSE Compliance Programme, any commitment of RSE is taken by=
the signature of one Representative granted by a proper Power of Attorney.=
Le informazioni contenute in questo messaggio di posta elettronica sono ris=
ervate e confidenziali e ne e' vietata la diffusione in qualsiasi modo o for=
ma. Qualora Lei non fosse la persona destinataria del presente messaggio, La=
invitiamo a non diffonderlo e ad eliminarlo, dandone gentilmente comunicazi=
one al mittente. The information included in this e-mail and any attachments=
are confidential and may also be privileged. If you are not the correct rec=
ipient, you are kindly requested to notify the sender immediately, to cancel=
it and not to disclose the contents to any other person.

```

• Prev by Date: Re: Calculating a simple integral
• Next by Date: Re: Removing list elements in x by list y
• Previous by thread: Re: Calculating a simple integral
• Next by thread: Re: Calculating a simple integral