MathGroup Archive 2004

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

Search the Archive

Re: Animating Solutions of NDSolve with respect to Initial Conditions

  • To: mathgroup at smc.vnet.net
  • Subject: [mg52543] Re: Animating Solutions of NDSolve with respect to Initial Conditions
  • From: Peter Pein <petsie at arcor.de>
  • Date: Wed, 1 Dec 2004 05:58:20 -0500 (EST)
  • References: <cohjd6$1of$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

Narasimham G.L. wrote:

> To see all circularly bent shapes of a strip fixed at one end it is
> animated as follows:
> 
> << Graphics`Animation`
> MovieParametricPlot[{Sin[t*s]/t , (1 - Cos[t*s])/t}, {s, 0, 2 Pi},
> {t, -1,1}, Frames -> 20, Axes -> False, AspectRatio -> Automatic, 
> PlotRange -> {{-3, 7}, {-3, 7}}];
> 
> Like the above and unlike the example below where solutions are
> superimposed in Show mode in a single frame, I like to animate an
> NDSolve output to varoious Boundary Conditions to see their effect
> dynamically in separate frames. How can this be done?
> 
> yvar = y /. 
> First /@ (NDSolve[{y'''[t] + y[t] == 0, y[0] == #1, y'[0] == #2, 
> y''[0] == #3}, y, {t, 0, 2 Pi}] & @@@ {{5, 1, -2}, {5, 
> 1, -1}, {5, 1, 0}, {5, 1, 1}, {5, 1, 2}});
> Plot[Evaluate[#[t] & /@ yvar], {t, 0, 5}];
> 
You nearly got it:

yvar = y /.
       First /@ (NDSolve[{y'''[t] + y[t] == 0, y[0] == #1, y'[0] == #2,
                   y''[0] == #3}, y, {t, 0, 2 Pi}] & @@@
             Table[{5, 1, p}, {p, -2, 2, .1}]);
Plot[#, {t, 0, 5}, PlotRange -> {-40, 10}] & /@
    Evaluate[#[t] & /@ yvar];

-- 
Peter Pein
10245 Berlin


  • Prev by Date: Re: FourierTransform
  • Next by Date: Re: Animating Solutions of NDSolve with respect to Initial Conditions
  • Previous by thread: Re: FourierTransform
  • Next by thread: Re: Animating Solutions of NDSolve with respect to Initial Conditions