sympyを用いて微分方程式の初期値問題を解きたい

前提

pythonで微分方程式「x′′ + 2x′ + 3x = te^(-t)cost, x(0) = 1, x′(0) = 2」をsympyを用いて計算するコードを作成しようとしていましたが、実行してもエラーすらも出てこないのでどの部分が悪いのかわからず困っていた。

実現したいこと

pythonで微分方程式「x′′ + 2x′ + 3x = te^(-t)cost, x(0) = 1, x′(0) = 2」をsympyを用いて実装したい。

発生している問題・エラーメッセージ

エラーまたは計算結果も何も出ないため、どの部分が間違っているのかわからず困っています。おそらくネイピア数「e」について行に間違いがあるのではないかと思い、試行錯誤しましたがだめでした。

該当のソースコード

python

import sympy as sym t = sym . symbols ('t')e = sym . symbols ('e')x = sym . Function ('x')(t)eq = sym . Eq ( sym . diff (x , t , 2) + 2*sym . diff (x , t) + 3* x , t*e**(-t)*sym . cos (t ))print ( sym . dsolve (eq , ics ={ x . subs (t ,0) :0 , sym . diff (x ,t ,1) . subs (t ,0) :1}) )

補足情報(FW/ツールのバージョンなど)

spyder(python3.8)

コメントを投稿

0 コメント