Can you do the graph in MatLab and share the code

Consider the following model for the propagation of a disease through a population as =-qs +p(N S-1), at = qsi rl. Here p gives the rate at which recovered people become susceptible again. Assume r/q<.>