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

时间:2023-03-09 03:24:23
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]]

  

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