Question: How to solve error occur in DEtools?

Attractor.mw
Hello, anybody could help with this error?

restart

interface(imaginaryunit = j);

I

(1)

with(DEtools):

k := .5; 1; tau := .95; 1; mu := 0.1e-1; 1; pi := 116.1; 1; vartheta := 0.8e-2; 1; phi := 0.25e-2; 1; epsilon := 0.2e-2; 1; rho := 0.5e-1; 1; beta := 0.115e-1; 1; chi := 0.598e-2; 1; q := .5; 1; eta := .2; 1; delta := .1; 1; alpha := 0.57e-1; 1; p := .2; 1; Upsilon := 1.2

.5

 

.95

 

0.1e-1

 

116.1

 

0.8e-2

 

0.25e-2

 

0.2e-2

 

0.5e-1

 

0.115e-1

 

0.598e-2

 

.5

 

.2

 

.1

 

0.57e-1

 

.2

 

1.2

(2)

lambda := k*tau*(C(t)*Upsilon+I(t))/(S(t)+V(t)+C(t)+I(t)+R(t)):

DEplot3d({d*C(t)/dt = rho*lambda*S(t)+rho*`ε`*lambda*V(t)+(1-q)*eta*I(t)-(mu+beta+chi)*C(t), d*I(t)/dt = (1-rho)*lambda*S(t)+(1-rho)*`ε`*lambda*V(t)+chi*C(t)-(mu+alpha+eta)*I(t), d*R(t)/dt = beta*C(t)+q*eta*I(t)-(mu+delta)*R(t), d*S(t)/dt = (1-p)*pi+phi*V(t)+delta*R(t)-(mu+lambda+`ϑ`)*S(t), d*V(t)/dt = p*pi+`ϑ`*S(t)-(lambda*`ε`+mu+phi)*V(t)}, {C(t), I(t), R(t), S(t), V(t)}, t = 0 .. 300, stepsize = .1, [[S(0) = 8200, V(0) = 2800, C(0) = 200, I(0) = 210, R(0) = 200]], linecolour = t, axes = BOXED, scene = [S, V, C, I, R])

Error, (in DEtools/DEplot/CheckInitial) the 'number' option must be specified before initial conditions

 

``


 

Download Attractor.mw

 

Please Wait...