how to plot output of NDSolve
- To: mathgroup at yoda.physics.unc.edu
- Subject: how to plot output of NDSolve
- From: Roque Donizete de Oliveira <oliveria at engin.umich.edu>
- Date: Tue, 15 Dec 92 18:10:54 -0500
Can anyone tell me how to get the output of NDSolve in a form
suitable for plotting ? For example, I want to plot X[t]
from t=d to t=2*d, where X[t] is the solution of something like:
In[71]:= d = 5
In[72]:= NDSolve[{ode1,ode2,ode3,ode4,
X[d] == x[d], Y[d] == y[d], Lambda[d] == lambda[d], Z[d] == lambda'[d]},
{X[t],Y[t],Lambda[t],Z[t]},{t,d,2 d}]
Out[72]= {{X[t] -> InterpolatingFunction[{5., 10.}, <>][t],
> Y[t] -> InterpolatingFunction[{5., 10.}, <>][t],
> Lambda[t] -> InterpolatingFunction[{5., 10.}, <>][t],
> Z[t] -> InterpolatingFunction[{5., 10.}, <>][t]}}
In[73]:= Out[%72] /. t->5
Out::intm: Machine-size integer expected at position 1 in
Out[X[t] -> InterpolatingFunction[{5., 10.}, <>][t]].
Out::intm: Machine-size integer expected at position 1 in
Out[Y[t] -> InterpolatingFunction[{5., 10.}, <>][t]].
Out::intm: Machine-size integer expected at position 1 in
Out[Lambda[t] -> InterpolatingFunction[{5., 10.}, <>][t]].
General::stop: Further output of Out::intm
will be suppressed during this calculation.
Out[73]= {{Out[X[5] -> 69.7962], Out[Y[5] -> 30.2038],
> Out[Lambda[5] -> 0.755095], Out[Z[5] -> 0.713493]}}
The final numerical result is fine. But all these messages ...
Thanks.
Roque