Make a single plotting function that plots all $4$ cell populations on a logarithmic scale, and caps all solutions to ${10}^{-6}$ or larger. It should also have a reasonable plot window size and automatically save to a filename of your choice. Plot $P$ in magenta, $L$ in blue, $C$ in red, and $M$ in black. Use grid lines. Make sure to create a legend. Rescale time to years $($approximating as $365$ days$).$

Write and use your code In python.

Part $2$

Describe what you see in dynamics of this new model. How does it differ from the model from Lecture $4($with "default" parameter values$)?$ What do you think drives the differences? In particular, point out differences in $M.$

For simplicity, call our original model the $"$PCM model," and call your new one the $"$PLCM model."