SIR Model with Demography

f[x_,y_]:=ρ*(1-x)-A*x*y;
g[x_,y_]:=y*(A*x-1-ρ);

A=4;
ρ=1;
plot1=StreamPlot[{f[x,y],g[x,y]},{x,-.1,1},{y,-.1,1},PlotRange->{{-.01,1},{-.01,1}},FrameLabel->{"x","y"},LabelStyle->Directive[Black, Large]];
plot2=ContourPlot[f[x,y]==0,{x,0,1},{y,0,1},ContourStyle->{Thick,Blue}];
plot3=ContourPlot[g[x,y]==0,{x,-.1,1},{y,-.1,1},ContourStyle->{Thick,Red}];
Show[plot1,plot2,plot3]