使用sympy微分方程的边界条件

Question

我想指定微分方程的边界条件。

*y"= 900(y - 1 + 2x) ; y(0)=5, y(2)=10* 

from sympy import * 
x=symbols('x') 
y, g = symbols('y g', cls=Function) 
diffeq = (Eq(y(x).diff(x, x) - 900*y(x) + 900, 1800*x),y(0):5,y(2)=10) 
A=dsolve(diffeq, y(x)) 
print A 

但它显示错误

diffeq = (Eq(y(x).diff(x, x) - 900*y(x) + 900, 1800*x),y(0):5,y(2)=10) 
                 ^
SyntaxError: invalid syntax 

请帮助。

Answer 1

边界条件作为字典传递给dsolve，通过ics命名参数。

这样：

from sympy import * 
x=symbols('x') 
f=symbols('f', cls=Function) 
dsolve(Eq(f(x).diff(x,x), 900*(f(x)-1+2*x)), f(x), ics={f(0):5, f(2):10}) 

您可以将最后一行粘贴到sympy live，以验证它的工作原理。答案是：

F（X）= C1 * E^-30x + C2 * E^30倍 - 2 + 1