aa79d79f-3b2f-4ca8-b20f-137c0e8ca024 42 Reg Curry 4 5 false false false 0 0 decimal

&[DATE] &[TIME] - &[FILENAME]

COMPARISON OF SIR FINITE DIFFERENCE AND ODE RUNGE-KUTTA SOLVER

γ 114 / :

Nus 330106^* : 3.30 108^*

R.o 7.0 :

β R.o γ * :

lastj 1000 :

δt 110 / :

COUPLED ODEs BEING SOLVED

FINITE DIFFERENCE

Ans I 1 el 1 : R 1 el 0 : S 1 el Nus : 31 line j 2 lastj range S j el S j 1 - el β I j 1 - el * S j 1 - el * Nus / δt * - : R j el R j 1 - el γ I j 1 - el * δt * + : I j el I j 1 - el β I j 1 - el * S j 1 - el * Nus / γ I j 1 - el * - (δt * + : 31 line for S R I 31 mat 21 line 21 line :

Ans1 Ans eval :

S Ans1 1 el :

I Ans1 3 el :

RUNGE KUTTA ODE SOLVER

R Ans1 2 el :

S rows 1000

RK t S t I t R 3 1 sys 100 100 Rkadapt :

t S' β t I * t S * Nus / - ≡

0 S Nus ≡

Days 1 lastj range :

Data Days 10 / I augment :

Data2 Days 10 / R augment :

t R' γ t I * ≡

0 R 0 ≡

Data3 Days 10 / S augment :

t I' β t I * t S * Nus / γ t I * - ≡

0 I 1 ≡

I max 1.9209 108^*

S.RK RK 2 col :

I.RK RK 3 col :

R.RK RK 4 col :

SRK RK 2 col :

IRK RK 3 col :

RK rows 101

RRK RK 4 col :

S.RK RK 1 col SRK augment :

I.RK RK 1 col IRK augment :

R.RK RK 1 col RRK augment :

Plot Data Data2 Data3 I.RK R.RK S.RK 61 sys :

Plot

R.RK 702 el 2.6203 108^*

Data I max 2 findrows 51.4 1.9209 108^* 12 mat

Data2 702 el 2.8488

I.RK I.RK max 2 findrows 50 1.908 108^* 12 mat

CASES.FD Days 10 / I R + augment :

CASES.RK RK 1 col IRK RRK + augment :

Plot2 CASES.FD CASES.RK 21 sys :

Plot2