MathGroup Archive 2008

[Date Index] [Thread Index] [Author Index]

Search the Archive

Definite integral problem

  • To: mathgroup at
  • Subject: [mg88518] Definite integral problem
  • From: UHAP023 at
  • Date: Wed, 7 May 2008 07:08:14 -0400 (EDT)

Dear All,
	Does anybody know any tricks to do the following definite 

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)/

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?

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

  • Prev by Date: Re: Calculations with Interpolating Functions
  • Next by Date: Re: Calculations with Interpolating Functions
  • Previous by thread: Re: Trouble computing conjugates
  • Next by thread: Error messages from FunctionInterpolation (bug?)