我有以下微分方程:
我试图找出y对x的二阶导数,预期的结果如下:
我试着用同情来表达上面的意思。这就是我尝试过的:
from sympy import *
x = symbols('x')
y = Function('y')
#initial value condition: y(1) = 0.5
# first derivative represented as d1
d1 = sin(x+y(x))/x
# second derivative represented as d2
d2 = diff(d1,x)
d2的结果是:
这就是我陷入困境的地方;如何用表达式d1替换(d/dx)(y(x))术语?
谢谢你的帮助。
发布于 2022-05-01 03:45:31
您可以使用subs替换导数:
In [13]: d1.diff(x).subs(y(x).diff(x), d1)
Out[13]:
⎛ sin(x + y(x))⎞
⎜1 + ─────────────⎟⋅cos(x + y(x))
⎝ x ⎠ sin(x + y(x))
───────────────────────────────── - ─────────────
x 2
x
https://stackoverflow.com/questions/72076247
复制