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

2744557306

1. a=0;% L; b- s3 e8 ?. v/ k& m9 u
   
2. b=1;
   : q2 O+ L1 ]3 p8 F7 g5 Z
3. N=40;
   # s1 z- B  v; @6 @/ Q2 O# ~) X8 R; `
4. af=1;
   ; w\" i& s  T+ _: B4 d/ Y\" d
5. Euler(a,b,N,af);
   & P. i6 n% e5 [. `
6. %Taylor(a,b,N,af);& T  A7 J4 r1 p, r4 H9 H
   
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
2.N：时间步数。
; P6 \7 r2 w- h; M( c/ b3.af：步长调整系数。
4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。
- k% B3 Z( u- |$ l% j


- D% E) q4 k# P
1 [. \, t1 d9 q( k