Advanced Search

Question: How to compute basic reproduction number?

Question:How to compute basic reproduction number?

Yee Voon 20 Maple
solve programming
February 06 2018
0

Does anyone know how to calculate basic reproduction number using maple coding? Or by any chance, anyone know how to solve it by hand with this complicated equations?
 

> 

restart
> 

interface(imaginaryunit = j)

I

 (1)
> 

lambda := k*tau*(C*Upsilon+I)/N

k*tau*(C*Upsilon+I)/N

 (2)
> 

eqn1 := (1-p)*Pi+phi*V+delta*R-(mu+lambda+`ϑ`)*S

(1-p)*Pi+phi*V+delta*R-(mu+k*tau*(C*Upsilon+I)/N+vartheta)*S

 (3)
> 

eqn2 := p*Pi+`ϑ`*S-(lambda*`ε`+mu+phi)*V

p*Pi+vartheta*S-(epsilon*k*tau*(C*Upsilon+I)/N+mu+phi)*V

 (4)
> 

eqn3 := rho*lambda*S+rho*`ε`*lambda*V+I*(1-q)*eta-(mu+beta+chi)*C

rho*k*tau*(C*Upsilon+I)*S/N+rho*epsilon*k*tau*(C*Upsilon+I)*V/N+(1-q)*eta*I-(mu+beta+chi)*C

 (5)
> 

eqn4 := (1-rho)*lambda*S+(1-rho)*`ε`*lambda*V+chi*C-I*(mu+alpha+eta)

(1-rho)*k*tau*(C*Upsilon+I)*S/N+(1-rho)*epsilon*k*tau*(C*Upsilon+I)*V/N+chi*C-(mu+alpha+eta)*I

 (6)
> 

eqn5 := beta*C+I*q*eta-(mu+delta)*R

beta*C+q*eta*I-(mu+delta)*R

 (7)
> 

``


 

Download Equation_for_basic_reproduction_number.mwEquation_for_basic_reproduction_number.mw
Please Wait...