restart; _local(gamma); _local(GAMMA); _local(Pi); _local(pi); _local(I); _local(D);
I
Warning, The imaginary unit, I, has been renamed _I
b := .45; mu[t] := 1.112; delta := .181; rho := 0.2e-1; beta[k] := .123; sigma := 0.9e-1; alpha := 0.312e-1; gamma := 0.14e-2; phi := .24; xi := .134; A[h] := .576; k[1] := 0.1e-2; beta[c] := 0.1e-1; mu[c] := 0.19e-2; eta := .557; z := .636; phi[c] := 0.57e-1;
0.45
1.112
0.181
0.02
0.123
0.09
0.0312
0.0014
0.24
0.134
0.576
0.001
0.01
0.0019
0.557
0.636
0.057
ODE1 := diff(B(T), T) = rho*b-mu[t]*B(T)-delta*B(T);
d
--- B(T) = 0.0090 - 1.293 B(T)
dT
ODE2 := diff(C(T), T) = delta*B(T)-mu[t]*C(T)+1-rho*b-beta[k]*C(T)*H(T)+sigma*C(T);
d
--- C(T) = 0.181 B(T) - 1.022 C(T) + 0.9910 - 0.123 C(T) H(T)
dT
ODE3 := diff(D(T), T) = beta[k]*C(T)*H(T)-(mu[t]+alpha+gamma)*C(T)-phi*C(T);
d
--- D(T) = 0.123 C(T) H(T) - 1.1446 C(T) - phi C(T)
dT
ODE4 := diff(E(T), T) = alpha*D(T)-xi*E(T)-mu[t]*E(T);
d
--- E(T) = 0.0312 D(T) - 1.246 E(T)
dT
ODE5 := diff(F(T), T) = phi*D(T)+xi*E(T)-mu[t]*F(T)-sigma*C(T);
d
--- F(T) = phi D(T) + 0.134 E(T) - 1.112 F(T) - 0.09 C(T)
dT
ODE6 := diff(G(T), T) = (1-k[1])*A[h]-beta[c]*G(T)*H(T)-mu[c]*H(T)-eta*J(T);
d
--- G(T) = 0.575424 - 0.01 G(T) H(T) - 0.0019 H(T) - 0.557 J(T)
dT
ODE7 := diff(H(T), T) = k[1]*A[h]+beta[c]*G(T)*H(T)-mu[c]*H(T)-z*H(T);
d
--- H(T) = 0.000576 + 0.01 G(T) H(T) - 0.6379 H(T)
dT
ODE8 := diff(J(T), T) = z*H(T)-(phi[c]+eta)*J(T);
d
--- J(T) = 0.636 H(T) - 0.614 J(T)
dT
ans := dsolve({ODE1, ODE2, ODE3, ODE4, ODE5, ODE6, ODE7, ODE8, B(0) = B0, C(0) = C0, D(0) = D0, E(0) = E0, F(0) = F0, G(0) = G0, H(0) = H0, J(0) = J0}, numeric, output = listprocedure);
Warning, The use of global variables in numerical ODE problems is deprecated, and will be removed in a future release. Use the 'parameters' argument instead (see ?dsolve,numeric,parameters)
[T = proc(T) ... end;, B(T) = proc(T) ... end;,
C(T) = proc(T) ... end;, D(T) = proc(T) ... end;,
E(T) = proc(T) ... end;, F(T) = proc(T) ... end;,
G(T) = proc(T) ... end;, H(T) = proc(T) ... end;,
J(T) = proc(T) ... end;]
B0 := 100; C0 := 90; D0 := 45; E0 := 38; F0 := 100; G0 := 45; H0 := 50; J0 := 20;
100
90
45
38
100
45
50
20
ans := dsolve([ODE1, ODE2, ODE3, ODE4, ODE5, ODE6, ODE7, ODE8, B(0) = B0, C(0) = C0, D(0) = D0, E(0) = E0, F(0) = F0, G(0) = G0, H(0) = H0, J(0) = J0], numeric, output = listprocedure);
Error, (in f) unable to store 'HFloat(450.486)-HFloat(90.0)*phi' when datatype=float[8]