Re: Numerical integration
- To: mathgroup at smc.vnet.net
- Subject: [mg79675] Re: [mg79621] Numerical integration
- From: Bob Hanlon <hanlonr at cox.net>
- Date: Wed, 1 Aug 2007 05:12:22 -0400 (EDT)
- Reply-to: hanlonr at cox.net
Version 6 appears to have a problem with this integral. $Version 6.0 for Mac OS X x86 (32-bit) (June 19, 2007) expr1 = Integrate[Log[1 + z^2]*(BesselJ[1, z]^2/z), {z, 0, Infinity}] (1/2)*HypergeometricPFQ[{1/2}, {1, 2}, 1]*(1 - 2*EulerGamma + Log[4]) expr2 = expr1 // FunctionExpand // Simplify (-(1/2))*(BesselI[0, 1]^2 - BesselI[1, 1]^2)* (-1 + 2*EulerGamma - Log[4]) expr1 // N 0.790559 expr2 // N 0.790559 NIntegrate[Log[1 + z^2]*(BesselJ[1, z]^2/z), {z, 0, Infinity}] 0.869058 Bob Hanlon ---- dimitris <dimmechan at yahoo.com> wrote: > In[3]:= > Integrate[Log[1 + z^2]*(BesselJ[1, z]^2/z), {z, 0, Infinity}] > N[%, 10] > > Out[3]= > MeijerG[{{1/2}, {}}, {{0, 0}, {-1}}, 1]/(2*Sqrt[Pi]) > > Out[4]= > 0.8732180258611361020606751916`10. > > On another CAS I took, > > convert("Integrate[Log[1 + z^2]*(BesselJ[1, z]^2/z), {z, 0, > Infinity}]",FromMma,evaluate); > evalf(%,20); > > 1/2 > 1/2 (2 Pi BesselI(0, 1) BesselK(0, 1) > > 1/2 / 1/2 > + 2 Pi BesselK(1, 1) BesselI(1, 1)) / Pi > / > > 0.87321802586113613925 > > Both CAS return the same symbolic result. > [An interesting challenge is to simplify the MeijerG > output of Mathematica to that of the other CAS] > > I want to check this symbolic result with NIntegrate. > I have "played around" with the options but I could get > "more closely" than > > In[16]:= > NIntegrate[Log[1+z^2]*(BesselJ[1,z]^2/z),{z, > 0, },MaxRecursion\[Rule]18]//InputForm > > >From In[16]:= > \!\(\* > RowBox[{\(NIntegrate::"slwcon > "\), \(\(:\)\(\ \)\), "\<\"Numerical integration > converging too slowly; suspect > one of the following: singularity, value of the integration being > 0, \ > oscillatory integrand, or insufficient WorkingPrecision. If your > integrand is \ > oscillatory try using the option Method->Oscillatory in NIntegrate. \ > \\!\\(\\*ButtonBox[\\\"More...\\\", ButtonStyle->\\\"RefGuideLinkText\\ > \", \ > ButtonFrame->None, ButtonData:>\\\"NIntegrate::slwcon\\\"]\\)\"\>"}]\) > > >From In[16]:= > \!\(\* > RowBox[{\(NIntegrate::"ncvb"\), \(\(:\)\(\ \)\), "\<\"NIntegrate > failed > to converge to prescribed accuracy > after \\!\\(19\\) recursive bisections in \\!\\(z\\) near \\!\\(z\ > \) = \ > \\!\\(35857.55603944016`\\). \\!\\(\\*ButtonBox[\\\"More...\\\", \ > ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ > ButtonData:>\\\"NIntegrate::ncvb\\\"]\\)\"\>"}]\) > > Out[16]//InputForm= > 0.8732193803103058 > > I would be very happy if someone pointed a tactic for perfroming > satisfactory numerical integration with mathematica for this integral. > I use mathematica 5.2 but you can use Mathematica 6 as well! > > The integral arose in another forum. There it was pointed out that > the performance of Mathematica 6 is buggy as regards numerical > integration. > > I look forward to seeing any replies. > > Greetings from burning Greece! > > Dimitris > > PS > > See here > > http://groups.google.gr/group/sci.math.symbolic/browse_thread/thread/57af36ff6f540a0d/a076ffbc412f974a?hl=el#a076ffbc412f974a > > for above mentioned thread. > >