MathGroup Archive 1992

[Date Index] [Thread Index] [Author Index]

Search the Archive

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





  • Prev by Date: Re: Academic pricing for Mma
  • Next by Date: window size
  • Previous by thread: 3D graphics
  • Next by thread: GraphicsArray