B_:={y, y2, S, I, W, R}$
FOR EACH EL_ IN B_ DO DEPEND EL_,T$

B1_:={bi, gam, mu, bw, k, xi, a}$
 %NUMBER OF STATES
NX_:=4$
 %NUMBER OF INPUTS
NU_:=0$
 %NUMBER OF OUTPUTS
NY_:=2$

C_:={df(S, t) = -bi*S*I - S*mu - S*bw*W + R*a + mu,
df(I, t) = bi*S*I - gam*I + S*bw*W - mu*I,
df(W, t) = xi*I - xi*W,
df(R, t) = gam*I - R*mu - R*a,
y = k*I,
y2 = S + R + I}$
FLAG_:=1$
SHOWTIME$
DAISY()$
SHOWTIME$
END$
