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

2744557306

1. a=0;! I7 A* i4 a8 _# K, L2 M
   
2. b=1;: ~) K% f; Q2 w1 c4 d# b
   
3. N=40;
   1 d2 }, U: H6 [: F, i! U
4. af=1;
   8 ~. W$ m' H1 h9 [) `6 D2 ?
5. Euler(a,b,N,af);& Q# Q\" W7 I  e3 r
   
6. %Taylor(a,b,N,af);4 f/ m* y' x) U1 }/ ~& \9 l0 |& o4 e
   
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
7 Z' v2 S8 ]2 g5 \2 k1 K: Q  }2.N：时间步数。
7 u% B6 ]; A: P& A: A+ B3.af：步长调整系数。
% `& J2 s2 O9 S8 e# w1 N4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。