範例:Stiffb 與 Stiffr 的增廣雅可比
使用
Jacob 函數計算常微分方程式 (ODE) 的增廣雅可比矩陣,然後將其提供作為求解器
Stiffb 與
Stiffr 的輸入。
1. 定義四個未知數系統:
2. 定義系統中的常數。
3. 定義初始值。
4. 定義對應至系統右側的向量函數 D(x,y)。
6. 呼叫 Stiffb 與 Stiffr 函數:
| 傳回的矩陣包含 5 欄,對應至數據點及四個未知數的解。 |
7. 從傳回的 Stiffb 與 Stiffr 矩陣萃取這四個未知數的解:
8. 從下列兩個函數繪製並比較每個未知數傳回的解:
此繪圖顯示兩個 ODE 求解器傳回相同的解。