How to | 微分代数方程式を解く方法
How to | 微分代数方程式を解く方法
Wolfram言語の微分方程式を解くための関数は,ユーザが予め処理しなくてもよい適切なアルゴリズムを自動的に選択して,多くの種類の微分代数方程式に適用できるようになっている.そのような方程式の1つに微分代数方程式がある.
代数的制約条件
を持つ
で与えられる微分代数方程式を設定する:
dae = {x'[t] - y[t] == 0, x[t] + 2y[t] == 0}DSolve を使って
について
および
の方程式を解く:
DSolve[dae, {x[t], y[t]}, t]DSolveの第1引数に含めることによって,初期条件
および
を 含めることもできる:
dae2 = {x''[t] == y[t], x[t] + y[t] == 0, x[Pi] == 1, x'[Pi] == 0}これでまた DSolveを使って解を求めることができる:
DSolve[dae2, {x[t], y[t]}, t]