shaibuo

10 Reputation

One Badge

4 years, 46 days

MaplePrimes Activity


These are replies submitted by shaibuo

@acer 

restart;
with(plots);
p := 10; omega := 0.3e-2; q := 0.234e-1; phi := .2; beta := .1; alpha := .2; sigma := 0.4e-2; tau := 0.2e-1; theta := .13; Omega := 0.23e-1; mu := 0.2e-1; k := 0.18e-1; eta := 0.2e-1; delta := 0.1e-1; unprotect(gamma); gamma := 0.2e-1; rho := 0.2e-1; PI := 0.4e-1; xi := 0.57e-1; a := 0.11e-1; l := 0.73e-1; chi := 0.95e-1; nu := 0.8e-1; unprotect(Zeta); Zeta := 0.26e-1; varphi := 0.2e-1;
funcs := diff(S(t), t) = delta*V(t)+PI*C__Sy(t)+xi*C__Asy(t)+varphi*R(t)+(1-p)*omega-(mu+(q+(1-q)*lambda)*phi+chi*lambda+nu)*S(t), diff(V(t), t) = p*omega-rho*lambda*V(t)+nu*S(t)-(delta+mu+kappa*lambda*eta+(1-kappa)*lambda*eta)*V(t), diff(C__Sy(t), t) = q*phi*S(t)+Omega*C__Sy(t)+eta*kappa*lambda*V(t)+l*i(t)-(PI+theta+mu+beta+sigma)*C__Sy(t), diff(C__Asy(t), t) = (1-q)*phi*S(t)+theta*C__Sy(t)+eta*lambda*(1-kappa)*V(t)+alpha*i(t)-(xi+Omega+mu+alpha+tau)*C__Asy(t), diff(i(t), t) = rho*lambda*V(t)+alpha*C__Asy(t)+beta*C__Sy(t)+chi*lambda*S(t)-(mu+Zeta+gamma+a+l)*i(t), diff(R(t), t) = gamma*i(t)+tau*C__Asy(t)+sigma*C__Sy(t)-(mu+varphi)*R(t), S(0) = 200, V(0) = 100, C__Sy(0) = 80, C__Asy(0) = 50, i(0) = 70, R(0) = 40;
p1 := dsolve({funcs}, type = numeric, abserr = 2.4);
Warning, The use of global variables in numerical ODE problems is deprecated, and will be removed in a future release. Use the 'parameters' argument instead (see ?dsolve,numeric,parameters)
p1i := odeplot(p1, [t, S(t)], 0 .. 120, numpoints = 100, labels = ["Time(Days)", "Susceptible  population"], labeldirections = [horizontal, vertical], style = line, axes = boxed, color = red);
Warning, cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up
 

Page 1 of 1