MathGroup Archive 2007

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

Search the Archive

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


  • Prev by Date: Re: fastest way to add up a billion numbers
  • Next by Date: Re: The Mathematica equivalent of Cons (as in Lisp)?
  • Previous by thread: Plot3D
  • Next by thread: Re: Plot3D