Now, write a function PLCM$\_$model$(Y,t,$ args$)$ that implements your updated model. It should return a vector of solutions in the order $(P,L,C,M).$ Don't forget to define your parameters and pass them as arguments to this new model function.

Part $2$

Now, solve and plot your model with the

default parameter values. Don't forget to use good labels so we can understand your plots. Plot time in units of years. Use logarithmic yscale as appropriate. Use tricks from lecture to resize your plate as necessary.