Definite integral problem
- To: mathgroup at smc.vnet.net
- Subject: [mg88518] Definite integral problem
- From: UHAP023 at alpha1.rhbnc.ac.uk
- Date: Wed, 7 May 2008 07:08:14 -0400 (EDT)
Dear All,
Does anybody know any tricks to do the following definite
integral?
J0part=(1/Sqrt[const/(Rvt^2 + Rx^2*Sin[theta]^2)^3] +
(3*Dc*Sec[theta]^2*(Rvt^2 + Rx^2*Sin[theta]^2)^3)/
(const*Rx^2))^(-1)
IntJ0term =
Integrate[J0part, {theta, 0, Pi/2},
Assumptions -> {Rx > 0, Rvt > 0, const > 0, Dc > 0}]
The integrand seems to be well behaved -- temporarily making some simple
numerical substitutions eg. Rx = 1; Rvt = 2; const = 3; Dc = 4; enables it
to plotted over the range of the integration limits and also integrated
with NIntegrate.
Plot[J0part, {theta, 0, Pi/2}] reveals a simple smooth curve.
J0part can be split into two partial fractions, neither of which
will symbolically integrate.
Curiously, different versions of Mathematica and trivially different
versions of the integrand get different responses from Integrate[] -- some
give up almost immediately, other rack up hours of CPU time without
producing a result.
Any ideas?
Thanks
Tom Crane.
Ps. I am running Mathematica 4.0.
Pps. The email address in the message header is a spam trap only.
--
Tom Crane, Dept. Physics, Royal Holloway, University of London, Egham Hill,
Egham, Surrey, TW20 0EX, England.
Email: T.Crane at rhul dot ac dot uk
Fax: +44 (0) 1784 472794