我有微分方程:dP/dt = kx(t)其中k是比例常数。我试着用解来求通解,但我不知道如何解释代码中的k。下面是我的代码,技术上可以工作,但没有考虑到k。如果你能告诉我如何编辑它,我将不胜感激。
import sympy as sp
t = sp.symbols('t')
x = sp.Function('x')
deq = sp.Eq(sp.diff(x(t),t), x(t))
xsoln = sp.dsolve(deq, x(t))
sp.pprint(xsoln)
让k
成为另一个符号。SymPy假定所有符号都是相互独立的。换句话说,相对于t
,k
被自动视为常数。
>>> t, k = symbols('t, k')
>>> x = Function('x')
>>> deq = Eq(diff(x(t), t), k*x(t))
>>> dsolve(deq, x(t))
Eq(x(t), C1*exp(k*t))