Abdi Sukamto Thesis

In[2]:=

Manipulate@ sol = Quiet @ NDSolve@8 s @tD b − β Hs@tD ∗ i@tDL − b s@tD, i @tD α r@tD − Hγ + bL i@tD, r @tD γ i@tD − Hγ + bL r@tD, s@0D a2, r@0D a3, i@0D a1<, 8s, i, r<, 8t, 0, Tfin<, MaxSteps → ∞D; Quiet @ TableForm@88TableForm@8 8ParametricPlot3D@Evaluate@8s@tD, i@tD, r@tD< ê. solD, 8t, 0, Tfin<, PlotPoints → 500, AxesLabel → 8s, i, r<, ImageSize → 200, FaceGrids → All, AxesEdge → 880, 0<, 81, 0<, 80, 0<<, ImagePadding → 25, BoxRatios → 81, 1, 1<, Boxed → True, BoxStyle → Directive@[email protected], FaceGridsStyle → Directive@[email protected], Darker@Gray, 0.1DD, PlotLabel → "Hi,s,rL phase plot"D<
2

abdiok.nb

length of simulation HtimeL

Hi,f L phase plot

i time plot 0.10

parameters b a

s 1.00 0.98 0.96 0.94 0.92 0.90 i 0.00 0.02 0.04 0.06 0.08 0.10

0.08 0.06 0.04 0.02

b

5

10

15

20

s time plot

g

Hi,cL phase plot

1.00 Out[2]=

initial conditions iH0L sH0L

r

0.98 0.015 0.010 0.005 0.000

0.96 0.94 0.92

rH0L

i 0.02 0.04 0.06 0.08 0.10

5

10

15

20

r time plot

Hs,rL phase plot r

0.015 0.015 0.010 0.005

0.010 0.005

0.90 0.92 0.94 0.96 0.98 1.00 5

10

15 20

s