dS[c] := Gamma*(P*In[a]+R[a]+S[a])-rho*S[c]-S[c]*beta[h]*In[v]/(S[c]+In[c]+R[c]+S[a]+In[a]+R[a]+In[m])-mu[h]*S[c];
dIn[c] := S[c]*beta[h]*In[v]/(S[c]+In[c]+R[c]+S[a]+In[a]+R[a]+In[m])-rho*In[c]-gamma*In[c]-mu[h]*In[c];
dR[c] := gamma*In[c]-rho*R[c]-R[c]*mu[h];
dS[a] := rho*S[c]-S[a]*beta[h]*In[v]/(S[c]+In[c]+R[c]+S[a]+In[a]+R[a]+In[m])-mu[h]*S[a];
dIn[a] := rho*In[c]+S[a]*beta[h]*In[v]/(S[c]+In[c]+R[c]+S[a]+In[a]+R[a]+In[m])-gamma*In[a]-mu[h]*In[a];
dR[a] := gamma*In[a]+rho*R[c]-R[a]*mu[h];
dIn[m] := Gamma*In[a]*(1-P)-mu[m]*In[m];
dS[v] := Gamma[v]-S[v]*beta[v]*(In[c]+In[a]+In[m])/(S[c]+In[c]+R[c]+S[a]+In[a]+R[a]+In[m])-mu[v]*S[v];
dIn[v] := S[v]*beta[v]*(In[c]+In[a]+In[m])/(S[c]+In[c]+R[c]+S[a]+In[a]+R[a]+In[m])-mu[v]*In[v];
solve({dIn[a]=0, dIn[c]=0, dIn[m]=0, dIn[v]=0, dR[a]=0, dR[c]=0, dS[a]=0, dS[c]=0, dS[v]=0}, {In[a], In[c], In[m], In[v], R[a], R[c], S[a], S[c], S[v]});
Warning, solutions may have been lost
↧