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]]