Re: Re: Integrate vs NIntegrate
- To: mathgroup at smc.vnet.net
- Subject: [mg46901] Re: [mg46877] Re: Integrate vs NIntegrate
- From: "Reza Malek-Madani" <research at usna.edu>
- Date: Sun, 14 Mar 2004 03:24:23 -0500 (EST)
- Sender: owner-wri-mathgroup at wolfram.com
This worked for me on Mathematica 4.2:
f = D[1/(1 + (1 + (a - b)^2)), a];
h[a_] := NIntegrate[f*b*(1 - b)^2, {b, 0, 1}];
FindRoot[h[a] == 0, {a, 0, 1}]
Output:
{a -> 0.397861}
Reza
-------------------------------------------------------------------------
Reza Malek-Madani, Director of Research and Scholarship
Research Office, MS 10m, Nimitz Room 17
589 McNair Road
U.S. Naval Academy
Annapolis MD 21402-5031
Phone: 410-293-2504 (FAX -2507), DSN: 281-2504
Email: research at usna.edu
--------------------------------------------------------------------------
>>> "Curt Fischer" <crf3 at po.cwru.edu> 03/12/04 11:39 PM >>>
Mukhtar Bekkali wrote:
> I am confused why NIntegrate misbehaves on such a simple function as
> mine.
>
> Here is what I have:
>
> In:
>
> f=D[1/(1+(1+(a-b)^2)),a];
> g=Integrate[f*b*(1-b)^2,{b,0,1},Assumptions->0<a<1];
> FindRoot[g==0,{a,0,1}]
>
> Out:
>
> a->0.397207
>
> However, since Integrate takes long, I tried to use NIntegrate
> instead and this is what I get
>
> In:
>
> f=D[1/(1+(1+(a-b)^2)),a];
> g:=NIntegrate[f*b*(1-b)^2,{b,0,1}];
> FindRoot[g==0,{a,0,1}]
>
> Out:
>
> a->1
>
> or, FindRoot+NIntegrate give me the upper boundary of a. If I
> abandon the secant method and turn to Newton, i.e. use
> FindRoot[g==0,{a,0.5}] instead then I get the message that Jacobian is
> singular at a=0.5 and get no solution. Perturbing the starting value
> of a does not help.
>
> What is going on here?
My guess is that NIntegrate does not do problems with symbolic
parameters,
and yet you are trying to make it calculate the integral of an
expression
involving a, with no numerical value for a defined.
--
Curt Fischer
> PS. Is there a way to get M5 to tackle the problem where:
> (1) I define some function f[x]:=NIntegrate[g[x,y],{y,0,1}], then
> (2) take the derivative of f[x] with respect to x, say h[x]:=f'[x]
> and then (3) Use FindRoot to find x such that h[x]==0