[英]Boundary conditions for a differential equation using sympy
我試圖指定微分方程的邊界條件。
*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
請幫助。
邊界條件通過ics
命名參數傳遞給dsolve
作為字典。
從而:
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 ^ 30x - 2x + 1
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.