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

  

 

posted on 2017-02-21 13:58  大宝pku  阅读(723)  评论(0编辑  收藏  举报

导航