忘记贴代码了
![[泪]](/static/emoticons/u6cea.png)
c1 = 0.3962;
c2 = 0.9680;
c3 = 0.4102;
c4 = 0.0220;
c5 = 0.0369;
c6 = 1.3380;
c7 = 0.2498;
c8 = 50.4122; c9 = 0.5316;
a = 2.6839*1.6839;
dh = D[h[t], t]
Kh[h_] = c8*(1 + (-c5*h)^c6)^(-c7*
c9)*(1 - (-c5*h)^(c6*c7))/(((1 + (-c5*h)^c6)^c7)^2);
dkh = D[Kh[h], h]
Solve[(c3 - c4)*((1/(1 + (-c5*h2)^c6))^c7) + c4 - c1 == 0, h2]
f[h_] = (c3 - c4)*(1/(1 + (-c5*h)^c6))^c7 + c4
df = D[f[h[t]], h[t], t]
DSolve[df == a*dkh, h[t], t]