MathGroup Archive 1997

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

Search the Archive

Re: Integration bugs in Mathematica 3.0

  • To: mathgroup at smc.vnet.net
  • Subject: [mg9276] Re: [mg9213] Integration bugs in Mathematica 3.0
  • From: Elias Saab <mathes at math.missouri.edu>
  • Date: Mon, 27 Oct 1997 02:46:47 -0500
  • Sender: owner-wri-mathgroup at wolfram.com

Did you try this one?
In[5]:= 2 Integrate[Sqrt[R^2-x^2-y^2],{x,-R,R},
        {y,-Sqrt[R^2-x^2],Sqrt[R^2-x^2]}]        //Timing

                            3
                      4 Pi R
Out[5]= {1.98 Second, -------}
                         3
The above example was executed in 2.2 on an SGI 133 MHZ. I tried it on a
Pentium II 300MHZ with Mathematica 3.01 and nothing happen in about 1
minute. I told Wolfram about this problem even before the release of
3.0.  It did not seem that they did anything about it. Elias Saab



On Fri, 24 Oct 1997, Carlos A. Felippa wrote:

> In a recent FEM project assigned as HW, about one third of the Integrate
> results (216 out of roughly 600)  were wrong.   Extracted sample:
> 
> fu=(2+Cos[phiQ])/Sqrt[6-2*Cos[phiQ]*(-2+Cos[theta])-4*Cos[theta]];
> BQ=Integrate[fu,{theta,0,Pi}]; 
> Print["BQ=",BQ//InputForm];
> fu=fu/.{phiQ->0};
> BQ=Integrate[fu,{theta,0,Pi}]; 
> Print["BQ=",BQ//InputForm];
> 
> Results:
> 
> BQ=0   (* wrong *)
> BQ=3*EllipticK[-3]  (* demonstrates bug *)
> 
> The students were not amused.  2.2 does it right, but most of them did
> not have access to it on campus machines.
> 
> Further, Integrate takes 2-5 times longer than in 2.2.
> 



  • Prev by Date: Re: Newbie question: big matrix calculations
  • Next by Date: Re: Re: Horse Race Puzzle (fwd)
  • Previous by thread: Re: Integration bugs in Mathematica 3.0
  • Next by thread: Statistics`ContinousDistributions` (integrating over a UniformDistribution)