Mathematica新特性Inactive, 求解复杂微分方程

Inactive阻止函数的计算, 求解微分方程有奇效

?

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

Block[{Integrate = Inactive[Integrate]},  DSolve[((H - h0)^(7/5) P0 (T - c t0))/(T0 (H - h0 - m*h[t])^(7/5)) +     1/2 \[Rho]*(h'[         t] - (h[t] m (m \[Pi] r^2 - S2))/(2 (H - h0 +             h[t] m) \[Pi] r^2) + (h[            t] m^2 (2 H \[Pi] r^2 - 2 h0 \[Pi] r^2 + h[t] m*\[Pi]*r^2 -              h[t]*S2 -             2 H*S2))/(2 (H - h0 +              h[t] m)^2 \[Pi] r^2) - (m (2 H \[Pi] r^2 -              2 h0 \[Pi] r^2 + h[t] m \[Pi] r^2 - h[t] S2 -              2 H S2))/(2 (H - h0 + h[t] m) \[Pi] r^2)*h'[t]) ==    P0 + g \[Rho] (h[t] +        h0 - (h[t] m (2 H \[Pi] r^2 - 2 h0 \[Pi] r^2 +             h[t] m \[Pi] r^2 - h[t] S2 - 2 H S2))/(2 (H - h0 +             h[t] m) \[Pi] r^2) - m*h[t]) +     1/2 \[Rho] Derivative[1][h][t]^2, h[t], t]]