我用p表示密度，用d表示微分，
变形得dt=cp*dtv/（kv^2—psl+q）
左边=cp/[kv^2（kv^2T—psl+q）]*(kv^2dT)
=cp[kv^2（kv^2T—psl+q）]dkv^2T
=cp[kv^2（kv^2T—psl+q）]d（kv^2T—psl+q）
=cp/（kv^2)dln|kv^2T—psl+q|
两边积分，得
t=cp/（kv^2)dln|kv^2T—psl+q|