I have a system of differential equations. I want to vary all the parameters, but the code shown below doesn't work for me. Could someone help me? I'm barely learning.
This is my code:
DynamicModule[{eqns, init, sol, H, L, F, Z, S, P, t},
Manipulate[eqns = {H'[t] =gamma1 (k1 F[t] + q1 K2 P[t])/(F[t] +
lambda ) P[t])-(1 - q2) gamma1 P[t](delta1+eta1) H[t],
L'(t)=eta1 H[t] (1- (L[t]+S[t])/(B))-L[t](tau1+epsilon1)L[t],
F'[t]=p1 tau1 L[t]-mu1 F[t],
Z'(t)=gamma2 q_2 P[t]-(delta2+eta2) Z[t],
S'[t]=eta2 Z[t] (1-(L[t]+S[t])/B))-S[t] (tau2+epsilon2),
P'[t]=p2 tau2 [t]-mu2 [t]};
init = {H[0] == H0, L[0] == L0, F[0] == F0, Z[0] == Z0, S[0] == S0,
P[0] == P0};
sol = NDSolve[{eqns, init}, {H, L, F, Z, S, P}, {t, 0, 300}];
Plot[Evaluate[{H[t], L[t], F[t], Z[t], S[t], P[t]} /. sol], {t, 0, 300},
PlotLegends -> Placed[{"H", "L", "F", "Z", "S", "P"}, Scaled[{0.9, 0.9}]]],
{{B, 1}, 0, 1200}, {{gamma1, 1}, 0, 2},{{gamma2,1},0,2},{{eta
1,1},0,2},{{eta2,1},0,2},{{delta1 ,1},0,2}, {{delta2,1},0,2},{{tau
1,1},0,2},{{tau2,1},0,2},{{epsilon1,1},0,2},{{epsilon2,1},0,2},
{{p1,1},0,2},{{p2,1},0,2},{{mu1,1},0,2},{{mu2,1},0,2},{{k1,1},0,2},
{{k2,1},0,2},{{q1,1},0,2},{{q2,1},0,2},
Delimiter, "initial conditions", {{H0, 0.5}, 0, 1}, {{L0, 0}, 0, 1}, {{F0,
1}, 0, 1}, {{Z0, 1}, 0, 1}, {{S0, 1}, 0, 1}, {{P0, 1}, 0, 1}]]
This is what appears to me:
Ctrl +A
andCtrl+C
very easy to me ? $\endgroup$