範例:求解一階 ODE 系統
使用解題指令群與
odesolve 函數,求解一階常微分方程式系統。
1. 定義解題區間的端點
2. 使用導數及一組初始條件,定義問題:
3. 在一個區間內繪製解:
使用 Rkadapt 函數
使用
Rkadapt 函數求解相同的微分方程式系統:
1. 定義決定任何解點 (t,Y) 之導數值向量的函數:
2. 定義 ODE 求解器的其他引數:
a. 自變數的初始值
b. 初始函數值的向量
c. [t0, t1] 上的解值數字
3. 使用 Rkadapt 函數求解解題矩陣。
4. 萃取自變數值
5. 萃取第一個解題函數值
6. 萃取第二個解題函數值
7. 萃取第三個解題函數值
8. 繪製這三個解題:
| odesolve 解題的繪圖與 Rkadapt 解題的繪圖幾乎相同。 |