mathematica一直运行没错误，大家帮忙看一下

上官

Clear[Am, As, Aa, \[Alpha], \[Rho], \[Theta]m, \[Theta]s, \
9 ^" b* A9 k/ k3 @\[CurlyPhi]m, \[CurlyPhi]s, \[Epsilon]]
\[Gamma]a = 0.1; \[Gamma]m = 0.15; \[Gamma]s =
1 - \[Gamma]a - \[Gamma]m;
4 k3 m  c: x- S3 Y; _4 l\[Epsilon] = 0.04; \[Alpha] = 0.3; \[Rho] = 0.04;
\[Theta]m = 0.75; \[Theta]s = 0.9;
1 @+ V& ?. Z8 D1 z; MgRate = 0.02;
2 |4 Z7 ]* s: }, VAm = (gRate + \[Rho])/\[Alpha]; Ba = 4; Bm = 1; Bs = 2.5;
ps = Bm/Bs; pa = Bm/Ba;
\[Delta] = 0.03;
B = \!\(TraditionalForm\`\*
FractionBox[
RowBox[{
5 H+ O% n! x* X0 zRowBox[{
RowBox[{
StyleBox["(",
SpanMinSize->1.,
& e- U5 b# ?( J1 f# T4 I/ cSpanMaxSize->1.],
RowBox[{"1", "\[Minus]", "\[Alpha]"}],
StyleBox[")",
SpanMinSize->1.,
SpanMaxSize->1.]}], "gRate"}], "+", "\[Rho]"}],
2 |$ q2 Z7 M4 x! e, ]9 Q      "\[Alpha]"] \[Minus] \[Delta]\);
cap = 10;
* P$ f) ^9 f) h4 j5 a7 s; Q/ ycsp = (pa*cap)/ps;
D = ((1 \[Minus] \[Alpha])*
+ e" |6 A, j: Z! c& S- o* N    gRate + \[Rho] - \[Alpha]*\[Delta])/(\[Rho] + gRate);
\[CurlyPhi]m = 0.1; \[CurlyPhi]s = 0.1;
Print["*** Initial Values ***"]
E0 = 1.5;
; C4 \- ?$ E- nK0 = E0/B;
hm0 = 0.25; hs0 = 0.25;(* initial values *)
: t! l2 ^  o5 c) |. ]0 ~\[Eta]m0 = hm0/K0; \[Eta]s0 = hs0/K0;
xm0 = (B*\[Gamma]m^\[Epsilon]*
# K2 S. Y0 j# e3 O. i   hm0^(\[Theta]m*(1 - \[Epsilon])))/(\[Gamma]a^\[Epsilon]*pa^(
    1 - \[Epsilon]) + \[Gamma]m^\[Epsilon]*
1 U+ y! Q! o* o- v    hm0^(\[Theta]m*(1 - \[Epsilon])) + \[Gamma]s^\[Epsilon]*(ps*
      hs0^\[Theta]s)^(1 - \[Epsilon]));
xs0 = (B*\[Gamma]s^\[Epsilon]*(ps*
3 S: [" P: S6 ]1 f- \+ ^     hs0^\[Theta]s)^(1 - \[Epsilon]))/(\[Gamma]a^\[Epsilon]*pa^(
    1 - \[Epsilon]) + \[Gamma]m^\[Epsilon]*
    hm0^(\[Theta]m*(1 - \[Epsilon])) + \[Gamma]s^\[Epsilon]*(ps*
      hs0^\[Theta]s)^(1 - \[Epsilon]));
Print["\[Eta]_{m,0}=" <> ToString[\[Eta]m0],
", \[Eta]_{s,0}=" <> ToString[\[Eta]s0],
", x_{m,0}=" <> ToString[xm0], ", x_{s,0}=" <> ToString[xs0]]
TT = 100;(* end time *)
  T- p" g# ]* @9 Z: e+ U(* Solve differential equations *)
Sol = NDSolve[{xs'[t] = (1 - \[Epsilon])*
! u2 L8 A) t. @2 E9 H4 {+ l     xs[t]*(   (1 - xs[t]/
         B)*(\[Theta]s*\[CurlyPhi]s*(xs[t]/(ps*\[Eta]s[t]) - 1) -
         xm[t]/B \[Theta]m*\[CurlyPhi]m*(xm[t]/\[Eta]m[t] - 1))),
   xm'[t] == (1 - \[Epsilon])*
5 P" A# }$ m/ Q     xm[t]*(   (1 - xm[t]/
: l# x2 s& F- d4 ]" U1 b( y          B)*\[Theta]m*\[CurlyPhi]m*(xm[t]/\[Eta]m[t] - 1) -
       xs[t]/B*\[Theta]s*\[CurlyPhi]s*(xs[t]/(ps*\[Eta]s[t]) -
          1) ), \[Eta]m'[
     t] == \[CurlyPhi]m*
) ~. N" u! U+ n& P6 x3 i' b: O      xm[t] - (\[CurlyPhi]m + gRate)*\[Eta]m[t], \[Eta]s'[
     t] == \[CurlyPhi]s*xs[t]/ps - (\[CurlyPhi]s + gRate)*\[Eta]s[t],
) w) S4 p6 t( K" A   K'[t] == gRate*K[t], hm[t] == \[Eta]m[t]*K[t],
   hs[t] == \[Eta]s[t]*K[t],
) h( z4 o2 p: y. U- o- O& u. S   Sa[t] == (\[Gamma]a^\[Epsilon]*(pa)^(1 - \[Epsilon]))/(\[Gamma]a^\
\[Epsilon]*pa^(1 - \[Epsilon]) + \[Gamma]m^\[Epsilon]*
4 K; n5 d3 z2 t3 K3 X! T" U( N       hm[t]^(\[Theta]m*(1 - \[Epsilon])) + \[Gamma]s^\[Epsilon]*(ps*
; R: ?6 q7 D1 f0 K& q         hs[t]^\[Theta]s)^(1 - \[Epsilon])) + (\[Gamma]m^\[Epsilon]*
/ ]7 n1 \/ @9 Z) S7 [      hm[t]^(\[Theta]m*(1 - \[Epsilon]))*pa*
      cap)/((\[Gamma]a^\[Epsilon]*pa^(
         1 - \[Epsilon]) + \[Gamma]m^\[Epsilon]*
! ^3 P" E) i; V: V' r- B         hm[t]^(\[Theta]m*(1 - \[Epsilon])) + \
\[Gamma]s^\[Epsilon]*(ps*hs[t]^\[Theta]s)^(1 - \[Epsilon]))*k (t)*
      xm (t)),
   Sm[t] == (\[Gamma]m^\[Epsilon]*
     hm[t]^(\[Theta]m*(1 - \[Epsilon])))/(\[Gamma]a^\[Epsilon]*pa^(
      1 - \[Epsilon]) + \[Gamma]m^\[Epsilon]*
' G- H$ G9 ?. r, R( J/ q- A( o      hm[t]^(\[Theta]m*(1 - \[Epsilon])) + \[Gamma]s^\[Epsilon]*(ps*
1 w9 ^; K- ?1 `! S4 ?        hs[t]^\[Theta]s)^(1 - \[Epsilon])),
   Ss[t] == (\[Gamma]s^\[Epsilon]*(ps*
        hs[t]^\[Theta]s)^(1 - \[Epsilon]))/(\[Gamma]a^\[Epsilon]*pa^(
       1 - \[Epsilon]) + \[Gamma]m^\[Epsilon]*
+ S& L7 e! E9 q9 z9 @2 m       hm[t]^(\[Theta]m*(1 - \[Epsilon])) + \[Gamma]s^\[Epsilon]*(ps*
8 x" I( o, w8 W% w/ N         hs[t]^\[Theta]s)^(1 - \[Epsilon])) - (\[Gamma]m^\[Epsilon]*
      hm[t]^(\[Theta]m*(1 - \[Epsilon]))*ps*
      csp)/((\[Gamma]a^\[Epsilon]*pa^(
         1 - \[Epsilon]) + \[Gamma]m^\[Epsilon]*
         hm[t]^(\[Theta]m*(1 - \[Epsilon])) + \
\[Gamma]s^\[Epsilon]*(ps*hs[t]^\[Theta]s)^(1 - \[Epsilon]))*k (t)*
      xm (t)), xm[0] == xm0,
   xs[0] == xs0, \[Eta]m[0] == \[Eta]m0, \[Eta]s[0] == \[Eta]s0,
9 @: I1 D' l- Y   K[0] == K0}, {xm, xs, \[Eta]m, \[Eta]s, K, hm, hs, Sa, Sm, Ss}, {t,
! r2 p4 o: D: l' j    0, TT}]
Plot[{Evaluate[Sa[t] /. Sol], Evaluate[Sm[t] /. Sol],
- l" B7 e: c  f8 \& z  Evaluate[Ss[t] /. Sol]}, {t, 0, TT}, AxesOrigin -> {0, 0},
* v, z0 Y8 p2 E/ i. P2 s7 A. i' R, h PlotRange -> {0., 0.8}, PlotStyle -> {Blue, Dashed, Dashing[{0.05}]}]
' [* `' g, Z/ q8 ^) e2 zPlot[{Evaluate[D*Sa[t] /. Sol],
+ k8 ]$ C5 B- {6 Y% q! z/ D8 ^  Evaluate[(D*Sm[t] + (\[Alpha]*(gRate + \[Delta]))/(\[Rho] +
       gRate)) /. Sol], Evaluate[D*Ss[t] /. Sol]}, {t, 0, TT},
3 j$ J3 U% e) j; }: q AxesOrigin -> {0, 0}, PlotRange -> {0., 0.8},
& N# `0 y( Z4 R. y3 y" F PlotStyle -> {Blue, Dashed, Dashing[{0.05}]}]




Set::wrsym: Symbol D is Protected.
9 q6 g3 w1 V+ c- R! y
+ B& L/ p8 [( w3 V0 F) WNDSolve::deqn: Equation or list of equations expected instead of 0.96 (1-6.66667 xs[t]) xs[t] (-0.5 xm[t] (-1+xm[t]/\[Eta]m[t])+0.09 (-1+(2.5 xs[t])/\[Eta]s[t])) in the first argument {0.96 (1-6.66667 xs[t]) xs[t] (-0.5 xm[t] (-1+xm[t]/\[Eta]m[<<1>>])+0.09 (-1+(2.5 xs[t])/\[Eta]s[<<1>>])),<<13>>,K[0]==10.}.
4 w  O) s6 A0 b' ?6 q5 N
1 F. O- t9 e  ~0 ~3 R2 P! m* @
( f2 H7 L) G1 Z: B/ X" C" d7 o* A
, Y1 w- I  c6 z+ Q5 ~/ \- \
" S0 A, Y& _# \/ q1 f3 E! r; U