Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2005
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2005

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

Search the Archive

Re: Solving double integral without Error function

  • To: mathgroup at smc.vnet.net
  • Subject: [mg57724] Re: Solving double integral without Error function
  • From: Paul Abbott <paul at physics.uwa.edu.au>
  • Date: Mon, 6 Jun 2005 04:21:33 -0400 (EDT)
  • Organization: The University of Western Australia
  • References: <d7rk3i$bkt$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

In article <d7rk3i$bkt$1 at smc.vnet.net>,
 "Xavier Brusset" <brusset at poms.ucl.ac.be> wrote:

> I have to solve this double integral:
> Phi[ q_ ]=Integrate[Integrate[y(x-q)f[x,y],{x,0,q}],{y,2,15}]
> where
> f[ x_ , y_ ]=1/(2 Pi Sqrt[27])Exp{-1/54 [4 x^2 - 6yx - 50x + 9y^2 -30 y +
> 325]}

There is a syntax error here in your definition of f.
 
> Mathematica 5.1 on windows fumbles a lot and finally stops giving an answer
> in q but using the Error function (I am using the Statistics Continuous
> distributions module).

That is irrelevant. Integrating exponentials of quadratic functions will 
lead to Erf. In general, there is no simple closed form for a second 
(indefinite) integral.

> The problem with the Error function is that when I combine this function Phi
> into another equation, I can't find a root in q.
> Which function of Mathematica should I use which circumvents the Error
> function in solving the double integral?
> Or, alternately, which root finding function could I use to find even a
> numeric instance of q?

There is a simple alternative way of proceeding. It is straightforward 
to compute the first and second derivatives of Phi with respect to q. 
You can use this with NDSolve to compute Phi[q]. See

  http://physics.uwa.edu.au/pub/Mathematica/MathGroup/ErfIntegral.nb

Cheers,
Paul

-- 
Paul Abbott                                      Phone: +61 8 6488 2734
School of Physics, M013                            Fax: +61 8 6488 1014
The University of Western Australia         (CRICOS Provider No 00126G)    
AUSTRALIA                               http://physics.uwa.edu.au/~paul
        http://InternationalMathematicaSymposium.org/IMS2005/


  • Prev by Date: Re: Viewing algebraic results
  • Next by Date: Re: Viewing algebraic results
  • Previous by thread: Solving double integral without Error function
  • Next by thread: Problem with Hyperlinks