Mathematica新特性Inactive, 求解复杂微分方程
Inactive阻止函数的计算, 求解微分方程有奇效
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]]