How to recreate graphs of a nonlinear system of ODEs with a condition containing the right hand side limit.

I need help recreating the graphs from Figure 2 and 3 in this paper, they relate to the section under "Single Chemotherapy"https://www.researchgate.net/publication/295685214_Mathematical_Modelling_and_Analysis_of_the_Tumor_Treatment_Regimens_with_Pulsed_Immunotherapy_and_Chemotherapy/fulltext/56f169b908ae1cb29a3d0e26/Mathematical-Modelling-and-Analysis-of-the-Tumor-Treatment-Regimens-with-Pulsed-Immunotherapy-and-Chemotherapy.pdf

Please help. I don't know much about MATLAB so unsure where to start

Answer 1

[Navya Seelam]

You can refer to the link below

https://www.mathworks.com/matlabcentral/answers/384298-how-to-solve-and-plot-system-of-nonlinear-differential-equations

[Rakeeza Malik]

My problem is that for one of my differential equations dz/dt, the condition I'm given is in the form of a right hand limit like so:

z(t^+) = z(t) + u