       Re: a problem about NDSolve and NIntegrate

• To: mathgroup at smc.vnet.net
• Subject: [mg80200] Re: a problem about NDSolve and NIntegrate
• From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
• Date: Tue, 14 Aug 2007 07:16:53 -0400 (EDT)
• Organization: The Open University, Milton Keynes, UK
• References: <f9emog\$lpq\$1@smc.vnet.net> <f9p6g1\$s1r\$1@smc.vnet.net>

```lynn wrote:
<snip>

> Thanks for your pointing out my faults.
>
> I miss-typed it when I tried to simplize my equation.
>
> But my original problem still exist,
>
> The following is the original equation,
>
> f[ _?NumberQ] := y[ ] /. NDSolve[{
>   2 y'[ ] y[ ] +
>       y''[ ] y[ ] -
>           (y'[ ])^2 ==   y[ ], y[0.00001] == 1, y'[0.00001] == =
> 0}, y,
> { ,
>                0.00001,  }, MaxSteps -> 1000][]
>
> then f=InterpolatingFunction[{{0.00001, 5.}}, <>][ ]
>
> and I want to integrate it y[ ] with a  ,
>
> \!\(g[w_?NumberQ,  _?NumberQ] := 2\
>       NIntegrate[\  \/y[ ]=D7 \/\@\( \^2 - w\^2\), { , w,  },
> MaxRecursion -> \
> 10]\)
>
> But I can't get correct answer when ask for such as g[1,2]...

Your post does not make any sense: do not use any special character set
when posting to Usenet. (Or use the Ingolf Dahl's "SetFaceAndFont"
palette. See http://library.wolfram.com/infocenter/MathSource/5752/)

I believe that f is defined more or less as follows

In:= f[c_?NumberQ] :=
y[x] /. NDSolve[{2 y'[x] y[x] + y''[x] y[x] - (y'[x])^2 == y[x],
y[0.00001] == 1, y'[0.00001] == 0}, y, {x , 0.00001, c},
MaxSteps -> 1000][]

So you can get an interpolating function with

In:= f

Out= InterpolatingFunction[][x]

and also plot it as in

In:= Plot[Evaluate@f, {x, 0.00001, 5}]

[... graphic deleted ...]

As far as I can see in the posted mess, g does not use f, and y[x] is
still *NOT* defined outside of f (this have already been pointed out by
Jens-Peer Kuska a couple of days ago); therefore g has no knowledge of
what could possibly be y[x].

--
Jean-Marc

```

• Prev by Date: Re: Drawing Potatoes and a Problem with Piecewise
• Next by Date: Re: Version 6 "Mathematica Book" - updated