dsolve コマンドにオプション numeric および method=mebdfi を指定すると、DAE 系の数値解を Modified Extended Backward Differentiation Equation Implicit 公式を使用して求めることができます。これはスティッフ法のため、スティッフ問題 を効率的に処理できます。通常の ODE 問題を解く場合も使用できますが、標準的な ODE ソルバを使用することをお勧めします (dsolve[numeric,IVP] を参照してください)。

mebdfi 方法で使用できるオプションには以下の種類があります。

'output' および 'known'

'output' オプションは dsolve からの出力を指定し、known オプションはユーザー定義の既知の関数を指定します。これらのオプションについては、dsolve[numeric] で説明します。

'startinit' および 'optimize'

これらのオプションは、数値計算の方法と動作を制御します。'startinit' オプションの使用法は dsolve[numeric,IVP] で説明するように IVP 問題と同様です。また、'optimize' オプションについては、dsolve[numeric] を参照してください。

'maxfun'

このオプションは求められる解を得るための手順の最大数を指定します。この直接ソルバがもっとも類似するのは dsolve[maxfun] です。デフォルトではこのオプションは無効になっています。

'abserr' および 'relerr'

これらのオプションは求解で目的とする精密さを指定します。詳細は、dsolve[Error_Control] で説明します。mebdfi のデフォルト値は、 と  です。

'minstep'、'maxstep'、および 'initstep'

これらのオプションは、計算手法で使用される手順のサイズに関してさらに詳しい制御を行います。詳細は 、dsolve[Error_Control] で説明します。デフォルトでは、 は、計算によって決定され、 は無効になっています。また、 が成り立ちます。 と指定した場合、この値が、低い次数で計算される最初の数ステップを実行できるだけの小さい適正なサイズであることを確認してください。

'maxord'

このオプションは、一連の計算で使用する公式の最大次数を 2 から 8 までの整数で指定します。デフォルト値は  8 ですが、難解な問題では、この値を  か  に指定する必要がある場合もあります。低い次数の方が安定しますが精度は低くなります。

数値計算を行うと mebdfi プロシージャの条件エラーに関するエラーメッセージが返されることがあります。

plots パッケージの odeplot 関数を使用すると結果をプロットすることができます。