我无法找到这个微分方程的特定解决方案。
from sympy import *
m = float(raw_input('Mass:\n> '))
g = 9.8
k = float(raw_input('Drag Coefficient:\n> '))
v = Function('v')
f1 = g * m
t = Symbol('t')
v = Function('v')
equation = dsolve(f1 - k * v(t) - m * Derivative(v(t)), 0)
print equation
对于m = 1000而k = .2则返回
Eq(f(t), C1*exp(-0.0002*t) + 49000.0)
这是正确的，但我希望在v（0）= 0时解决的等式应该返回
Eq(f(t), 49000*(1-exp(-0.0002*t))