Re: Mathematica gives bad integral ??
- To: mathgroup at smc.vnet.net
- Subject: [mg24330] Re: Mathematica gives bad integral ??
- From: zhl67 at my-deja.com
- Date: Sun, 9 Jul 2000 04:52:39 -0400 (EDT)
- References: <8k3n38$3pk@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
John,
I am confused with your post also: what you seem to have wished to do
is to perform a symbolic indefinite integration with Mathematica and
compare the result with existing text. As you have understood, there
are many equivalent forms for the same algebraic expression. Then you
seem to wish to confirm yourself with numerical tests. At this point,
however, there is some ambiguities in your description: at which
interval for x did you evaluate the DEFINITE integration? it is clear
that evaluating the numeric value of an indefinite integral is
meaningless: there is always an arbitrary integration constant in it
(which most symbolic algebra systems omit in their output).
I tested the same integration and get the following result:
In[1]:=Integrate[1/Sqrt[1-Sin[2x]],{x,0,Pi/8}]//N
Out[1]:=0.518675+0. I
In[2]:=NIntegrate[1/Sqrt[1-Sin[2x]],{x,0,Pi/8}]
Out[2]:=0.518675
In[3]:=Integrate[1/Sqrt[1-Sin[2x]],{x,0,Pi/16}]//N
Out[1]:=0.220268+0. I
In[2]:=NIntegrate[1/Sqrt[1-Sin[2x]],{x,0,Pi/16}]
Out[2]:=0.220268
Mathematica seems to give consistent result!
On the other hand, evaluating the difference of Schaum's result gives
(f[x_] is defined to be the Schaum's integral)
In[6]:=f[Pi/8]-f[0]
Out[6]:=0.518675
In[7]:=f[Pi/16]-f[0]
Out[7]:=0.220268
This also give the same result! (All these are with Ver 4.0)
L Zhao
Sent via Deja.com http://www.deja.com/
Before you buy.