コマンド差分シンボルD2は、二階微分を表すD、で表すことができ、D3は第3次微分を表し、。
厳密解を探します
1.微分方程式
R = dsolve（ 'eqn1'、 'eqn2'、...、 'COND1'、 'COND2'、...、 'VAR'）。
次のように彼は説明した：eqniはI番目の微分方程式を表し、I THコンディ初期条件を表し、varはtと、デフォルトでは、微分方程式の独立変数を表します。
>> dsolve（ 'のDy = 3 * X ^ 2'、 'Y（0）= 2'、 'X'） ANS =
X ^ 3 + 2
2.微分方程式
>> [X、Y] = dsolve（ 'Dxの= Y'、 'D2Y-のDy = 0'、 'X（0）= 2'、 'Y（0）= 1'、 'のDy（0）= 1' ） X = EXP（T）+ 1 、Y = EXP（T）
3.微分方程式を解く
>> [X、Y] = dsolve（ 'Dxを+ 5 * X + Y =のEXP（T）'、 'DY-X-3 *のY = 0'、 'X（0）= 1'、 'Y（0 ）= 0' 、 'T'） のx =の EXP（T×（15 ^（1/2） - 1））* ^（1/2（15） - 4）*（（13 * 15 ^（1/2 ））/ 330 - EXP（2 * T - 15 ^（1/2）* T）*（15 ^（1/2）/ 165 + 1/22）+ 1/22） - EXP（-t×（15 ^（1/2）+ 1））*（EXP（2 * T + 15 ^（1/2）* T）*（15 ^（1/2）/ 165 - 1/22）+（15 ^（1 / 2）*（15 ^（1/2） - 13））/ 330）*（15 ^（1/2）+ 4） Y =の EXP（-t×（15 ^（1/2）+ 1）） *（EXP（2 * T + 15 ^（1/2）* T）*（15 ^（1/2）/ 165 - 1/22）+（15 ^（1/2）*（15 ^（1 / 2） - 13））/ 330）+ EXP（T×（15 ^（1/2） - 1））*（（13 * 15 ^（1/2））/ 330 - EXP（2 * T - 15 ^ （1/2）* T）*（15 ^（1/2）/ 165 + 1/22）+ 1/22） >> ezplot（x、y）は
その差をプロットezplot
横軸としてXとプロット（x、y）は、yがプロットされ、縦軸であります
Y1、Y2 ...要素にX横座標値にプロット（X、Y1、X、Y2、...）は、プロットされた曲線の複数の座標
xのプロット、式中Yが（X）式の知られている、またはそのようなY = Fとして形成されています
そしてezplot陰関数は、グラフィックスを描画することで、Yは、F（x）はこの関数の、データ準備なしexplotは、関数グラフを直接描か=ようなパターンを書き込むことができないフォームF（X、Y）= 0であります
近似解
ODEソルバ
| ソルバー | 質問タイプ | 精度 | ときに使用します |
|---|---|---|---|
ode45 |
非剛体 | で | ほとんどの場合、あなたは最初のソルバー試してみてください |
ode23 |
低いです | 緩い許容範囲の問題についてや適度な剛性ケースに |
|
ode113 |
低〜高 | 計算オーバヘッドが大量に必要とする厳しい制限を伴うまたはODE関数の場合の許容誤差の問題のために、 |
|
ode15s |
剛性 | 中低〜 | 場合は |
ode23s |
低いです | 問題のエラーのより寛大なマージンのために、
品質行列がある場合、それは定数行列でなければなりません。 |
|
ode23t |
低いです | ちょうど適度な剛性の問題について、あなたは数値減衰が使用、ではない解決する必要があります
|
|
ode23tb |
低いです | そして |
|
ode15i |
完全隐式 | 低いです | 完全暗黙発行F（T、Y、Y 'のための ）= 0 と差分インデックス微分代数方程式（DAE）、使用 |
微分方程式1. 
インライン（）は簡単な言葉で機能を定義するために使用され、関数がインラインを使用して定義され
、bについて、XがYを代入することによって得ることができます
>> F =インライン（ '* X + B'、 'A'、 'B'、 'X'）。
>> F（1,2,3）
ANS =
5
研究ルンゲクッタ原則。。以上であるために