Plot3D
- To: mathgroup at smc.vnet.net
- Subject: [mg73858] Plot3D
- From: yjkingma at telus.net
- Date: Fri, 2 Mar 2007 06:29:55 -0500 (EST)
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]