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