Re: Integrate vs NIntegrate
- To: mathgroup at smc.vnet.net
- Subject: [mg89144] Re: [mg89129] Integrate vs NIntegrate
- From: DrMajorBob <drmajorbob at att.net>
- Date: Tue, 27 May 2008 07:13:31 -0400 (EDT)
- References: <9489155.1211802216995.JavaMail.root@m08>
- Reply-to: drmajorbob at longhorns.com
The fundamental theorem of calculus frequently does not apply when, as in
your example, Integrate is free to return a Complex-valued function.
(A better question is why you put in all those extraneous stars to inhibit
copy/pasting the code.)
Despite that, here's a slightly different approach where the FTC does
apply:
q = -1/2; a = 0; b = 3;
h[x_] = (1 + x^3)^q;
f[x_] = Integrate[h[t], {t, 0, x}]
If[x >= -1, (x^3)^(1/3) Hypergeometric2F1[1/3, 1/2, 4/3, -x^3],
Integrate[1/Sqrt[1 + t^3], {t, 0, x}, Assumptions -> x < -1]]
Res1 = N[f[b] - f[a]]
Res2 = NIntegrate[h[x], {x, a, b}]
1.65267
1.65267
Contrast that definition of f with yours:
Integrate[h[x], x]
(1/(3^(1/4) Sqrt[1 + x^3]))2 (-1)^(
1/6) Sqrt[-(-1)^(1/6) ((-1)^(2/3) + x)] Sqrt[
1 + (-1)^(1/3) x + (-1)^(2/3) x^2]
EllipticF[ArcSin[Sqrt[-(-1)^(5/6) (1 + x)]/3^(1/4)], (-1)^(1/3)]
If you plot both versions from a to b (use PlotRange->All to be sure),
you'll find that your version is complex on part of the range (between 2
and 3), but mine is real.
Bobby
On Mon, 26 May 2008 05:23:35 -0500, Armen Kocharyan
<armen.kocharyan at gmail.com> wrote:
> *q=-1/2; a=0; b=3;*
> *h[x_]=(1+x^3)^q;*
> *f[x_]=Integrate[h[x],x];*
> *Res1=N[f[b]-f[a]]*
> *Res2=NIntegrate[h[x],{x,a,b}]*
--
DrMajorBob at longhorns.com