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