Euler Taylor Runge-Kutta 方法求解常微分方程

2744557306

1. a=0;
   # _) |1 ]1 V4 U* V\" `
2. b=1;5 N1 P\" s% H* g1 T/ Z- w$ j: j
   
3. N=40;
   6 t7 _6 I! H1 z+ X& i5 h/ \
4. af=1;
   7 c! w, @9 I! {1 e\" \( `
5. Euler(a,b,N,af);
     p% ~3 y# `$ I5 b4 w4 e, d
6. %Taylor(a,b,N,af);
   8 w6 p2 `* q% g$ F
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
" P% s' K: ^% c2.N：时间步数。
; G8 R, @" F/ {, G# |3.af：步长调整系数。
1 E/ l2 H% F4 |3 G4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。
) p) Y9 V3 m- ~6 y$ \- D


8 V. Q5 O- T8 F9 Y