MathGroup Archive 2004

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

Search the Archive

Re: Another Integrate error


This works correctly in version 5.1 (MS Windows).

Steve Luttrell

"Ray Koopman" <koopman at sfu.ca> wrote in message 
news:comgts$9gn$1 at smc.vnet.net...
> This is a bivariate normal probability density:
>
> In[1]:= f[x_,y_,r_] = Exp[(x^2 + y^2 - 2*r*x*y)/(r^2 - 1)] /
>                      (Pi Sqrt[1 - r^2]);
>
> It should integrate to 1:
>
> In[2]:= Integrate[f[x,y,r],{x,-Infinity,Infinity},{y,-Infinity,Infinity}]
> Out[2]= 0
>
> In[3]:= Integrate[f[x,y,r],{x,-Infinity,Infinity},{y,-Infinity,Infinity},
>                  Assumptions -> {Element[r,Reals], -1 < r < 1}]
> Out[3]= 1
>
> In this case I happened to know that the integal should be 1,
> but what about more obscure cases?
> When can I believe a result from Integrate?
>
> In[4]:= $Version
> Out[4]= 5.0 for Mac OS X (November 19, 2003)
> 



  • Prev by Date: Re: .WAV file headers
  • Next by Date: Re: Another Integrate error
  • Previous by thread: Re: Another Integrate error
  • Next by thread: Re: Another Integrate error