I want to plot a phase space for a autonomous system equation with varying parameters using Matlab. Please help me on how to go about it. The system of equation is given below:
dS/dt =b - ds-KSI/(1 +αS + βI)+ γR
dI/dt = ds - KSI/(1 +αS + βI) - (d + μ)I
dR/dt = μI - (d+ γ)R
where S, I, R are susceptibe, infective, and recovery respectively. dS/dt, dI/dt, dR/dt are total derivatives.b is recruitment rate, infact all others are parameters.