Question: How to Plot a figure from Differential Equation (Implicit functions)

I have Differential equation sytems as follow:

r := 10^(-2);k := 10^12;c1 := 5*10^(-11);c2 := 10^(-13);mu := 8;lambda := 4.16;a := 2*10^3;q := 300;s0 := 3*10^5;d := 10^(-3);rho := 10^(-12);gama := 10^2;delta := 10^4;b := 5*10^6;

eq1 := {N(0) = 2*10^10, P(0) = 5*10^7, Q(0) = 0, diff(N(t), t) = r*N(t)*(1 - N(t)/k) - c1*N(t)*P(t) - mu*N(t)*Q(t)/(a + Q(t)), diff(P(t), t) = s(t) + s0 - d*P(t) + rho*N(t)*P(t)/(gama + N(t)) - c2*N(t)*P(t) - delta*P(t)*Q(t)/(b + Q(t)), diff(Q(t), t) = q(t) - lambda*Q(t)}. I replace I with P.

N(t):The number of cancer cells; P(t): The number of immune cells (Healthy T-Lymphocuyes); Q(t): The amount of chemotherapeutic agent in bloodstream; q(t): Source terms standing for chemotherapy.

My question is how to get the first two figures below based on these information: q_infinity means constant chemotherapy infusion, while s_infinity means constant immunotherpy.

I got the plot for Q(t) as follow:

Please anyone has any suggestions to help. Thank you very much.