Wolfram Cloud Document

ClearAll["Global`*"];

f1[S_,e_,I1_,I2_,R_]:=b*(S+R)-d*S-g*S*

S+e+I1+I2+R

K

-bcd*S*(I1+I2);

f2[S_,e_,I1_,I2_,R_]:=-d*e-g*e*

S+e+I1+I2+R

K

+bcd*S*(I1+I2)-scd1*e;

f3[S_,e_,I1_,I2_,R_]:=-d*I1-g*I1*

S+e+I1+I2+R

K

+(1-f)*scd1*e-scd2*I1;

f4[S_,e_,I1_,I2_,R_]:=-d*I2-g*I2*

S+e+I1+I2+R

K

+f*scd1*e-acd*I2;

f5[S_,e_,I1_,I2_,R_]:=-d*R-g*R*

S+e+I1+I2+R

K

+scd2*I1;

Expand[{f1[S,e,I1,I2,R]0,f2[S,e,I1,I2,R]0,f3[S,e,I1,I2,R]0,f4[S,e,I1,I2,R]0,f5[S,e,I1,I2,R]0}]

bR+bS-dS-bcdI1S-bcdI2S-

egS

K

-

gI1S

K

-

gI2S

K

-

gRS

K

-

g

2

S

K

0,-de-

2

e

g

K

-

egI1

K

-

egI2

K

-

egR

K

+bcdI1S+bcdI2S-

egS

K

-escd10,-dI1-

egI1

K

-

g

2

I1

K

-

gI1I2

K

-

gI1R

K

-

gI1S

K

+escd1-efscd1-I1scd20,-acdI2-dI2-

egI2

K

-

gI1I2

K

-

g

2

I2

K

-

gI2R

K

-

gI2S

K

+efscd10,-dR-

egR

K

-

gI1R

K

-

gI2R

K

-

g

2

R

K

-

gRS

K

+I1scd20

NSolve[%,{S,e,I1,I2,R},Reals]