Re: Plot3D
- To: mathgroup at smc.vnet.net
- Subject: [mg73912] Re: Plot3D
- From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
- Date: Sat, 3 Mar 2007 00:57:44 -0500 (EST)
- Organization: The Open University, Milton Keynes, UK
- References: <es910n$2po$1@smc.vnet.net>
yjkingma at telus.net wrote: > For the system below I would like to generate a 3D plot of i[t] when the p= > arameter b takes > values from say .2 to 3. Please help. > jack > a==.0002; b==1; c==1; d==.002; > > so==NDSolve[{s'[t]====-a*s[t]*i[t]-s[t]*i[t]*d*UnitStep[t-b]*(1-Exp[-.1(t-b= > )]), > > i'[t]Sa*s[t]*i[t]-c*i[t],r'[t]Sc*i[t], > > s[0]S10000,i[0]S1,r[0]====0},{s[t],i[t],r[t]}, > > {t,0,10},MaxSteps=AE1000,StartingStepSize=AE.001]; > > f==Plot[Evaluate[{s[t],i[t],(10000-.5 r[t])}/.so],{t,0,10},PlotPoints=AE500= > ,GridLines=AEAutomatic,PlotStyle=AE{{Thickness[.008],RGBColor[1,0,0]},{Thic= > kness[.008],RGBColor[0,1,0]},{Thickness[.008],RGBColor[0,0,1]}},Frame=AETru= > e] > Could you, please, post some syntacticly correct code? Regards, Jean-Marc