Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
1997
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 1997

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

Search the Archive

Re: Iterating NDSolve

  • To: mathgroup at smc.vnet.net
  • Subject: [mg7145] Re: [mg7117] Iterating NDSolve
  • From: Andrei Constantinescu <constant at athena.polytechnique.fr>
  • Date: Tue, 13 May 1997 01:57:45 -0400 (EDT)
  • Sender: owner-wri-mathgroup at wolfram.com

 Hi Darmiiss,

 I did something similar some time ago using Map. Take a look 
 a the code below and may be it helps:

 tmax = 20;

 sol = Map[
     NDSolve[  { l'[t] == B[r[t], l[t]] , r'[t] ==  0.35 (# -  r[t] - l[t]) , 
        l[0] == 0, r[0] == 0 }, 
        {r[t],l[t]}, {t, 0, tmax} ] &, 
       Union[ Table[ i, {i, 2.95, 3.0, 0.002} ], Table[ i, {i, 3.40, 3.45, 0.002} ] ]
     ]
 traj = Map[ 
  ParametricPlot[ {#[[1,1,2]], #[[1,2,2]]}  , {t, 0, tmax}, 
               DisplayFunction -> Identity] &, sol ]
 b0 = ContourPlot[ B[r,l], {r,0.01,5}, {l,0,5}, 
       Contours -> {0} , ContourShading -> False,
       ContourStyle -> {Hue[0]},
       PlotPoints -> 25,
       DisplayFunction -> Identity]
 w0 = ContourPlot[ W[r,l], {r,0.01,4}, {l,0,4}, 
       Contours -> Range[10] , ContourShading -> False,
       ContourStyle -> {Hue[0.5]},
       PlotPoints -> 25,
       DisplayFunction -> Identity]
 Show[ b0,w0,Sequence[traj],  DisplayFunction -> $DisplayFunction ,
      FrameLabel -> {"r - vitesse d'avancement", "l - lit des debris" , "" , "" } ,
      DefaultFont -> {"Helvetica-Bold", 12} ]

 Cheers,

   Andrei

______________________________________________________________________
  Andrei Constantinescu               constant at athena.polytechnique.fr   
                                                                       
  LMS Ecole Polytechnique                  tel:   (33)-1-69.33.33.30
  91128 PALAISEAU cedex - FRANCE           fax:   (33)-1-69.33.30.26    



  • Prev by Date: Re: Iterating NDSolve
  • Next by Date: Re: f + g
  • Previous by thread: Re: Iterating NDSolve
  • Next by thread: Re: Graphics`Graphics3D`Shadow