Re: Performance improvement needed - Help.
- To: mathgroup at smc.vnet.net
- Subject: [mg63734] Re: Performance improvement needed - Help.
- From: "Tonybony" <aneves at gmail.com>
- Date: Thu, 12 Jan 2006 04:09:27 -0500 (EST)
- References: <dpiotv$l11$1@smc.vnet.net><dptdpr$l1$1@smc.vnet.net> <dq03go$ri5$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Dear Paul, The numeric interal I am trying to evaluate is of the type: I1[n_,m_,X_] := NIntegrate[Sqrt[Cos[a]] Exp[-Sin[a]^2 c^2] (LegendreP[n, m, Cos[a]] (m^2 BesselJ[m, X Sin[a]](1 - Cos[a])/(X Sin[a]) + m Cos[a] BesselJ[m+1, X Sin[a]]) + (m BesselJ[m, X Sin[a]]/(X Sin[a]) - BesselJ[m + 1, X Sin[a]])(n + m)(n - m + 1)Sin[a] LegendreP[n,m - 1, Cos[a]]), {a, 0, 1.1}, AccuracyGoal -> 6, Compiled -> False]; This is another problem that I would like to address, due do the oscilatory nature of Legendre for high 'n', how could I compute this integral safelly, say for n=40 m=1? Thanks, Antonio