Problem with combination NDSolve NIntegrate
- To: mathgroup at smc.vnet.net
- Subject: [mg72739] Problem with combination NDSolve NIntegrate
- From: Paul Schneider <paul.schneider at wu-wien.ac.at>
- Date: Wed, 17 Jan 2007 06:44:44 -0500 (EST)
I am using version 5.2
NIntegrate tries to evaluate the integrand symbolically. I try to
circumvent this by using ?NumberQ but without success. Please try the
following toy example:
trial[a_?NumberQ] :=
NDSolve[{ti'[t] == ti[t] , ti[0] == a} , ti, {t, 0, 10}];
NIntegrate[ (ti[5] /. trial[a]), {a, 0, 10}]
Plot[ (ti[5] /. trial[a]), {a, 0, 10}]
NIntegrate does not evaluate to a number and throws a warning, while I
get a perfect plot from Plot.
Does anybody know what to do?
Cheers,
Paul