I am a novice. I want to get the plot of ODE, but I can only get the axes by running the following code. Where is the error,? Hope someone can help, thank you!
ClearAll; r1 = 0.00167; r2 = 0.0028; e1 = 0.0003; e2 = 0.0005; k1 = 50000; k2 = \ 150000; h2 = 23/200; h1 = 1/16; b1 = 0.000056; b2 = 0.00001; d = 0.000114; s = 10; \[Delta]12 \ = 0.00002; \[Delta]21 = 3.5; G1 = 8; a = 3.21699*10^-6; m = 5732.3; Clear[x, y, z]; equ = NDSolve[{x'[t] == r1 x[t] (1 - (x[t] + \[Delta]12 y[t])/k1) - e1 x[t] z[ t]/(1 + h1 e1 x[t]/(1 - a x[t]) + h2 e2 y[t] m^s/(m^s + x[t]^s)), y'[t] == r2 y[t] (1 - (y[t] + \[Delta]21 x[t])/k2) - e2 y[t] z[ t] m^s/(m^s + x[t]^s)/(1 + h1 e1 x[t]/(1 - a x[t]) + h2 e2 y[t] m^s/(m^s + x[t]^s)), z'[t] == (b1 e1 x[t] + b2 e2 y[t] m^s/(m^s + x[t]^s)/(1 + h1 e1 x[t]/(1 - a x[t]) + h2 e2 y[t] m^s/(m^s + x[t]^s)) - d) z[t], x[0] == 8162, y[0] == 117894, z[0] == 6}, {x[t], y[t], z[t]}, {t, 0, 10000}, AccuracyGoal -> 30, PrecisionGoal -> 30, WorkingPrecision -> 50, MaxSteps -> 5000] Plot[Evaluate[{z[t]} /. equ], {t, 0, 400}]