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

2744557306

1. a=0;
     m- }: D0 n+ Q/ Y
2. b=1;
   . x+ L5 f4 r: {: @& K7 n
3. N=40;
   7 F' C# n. [2 M8 V% _2 @' r+ |
4. af=1;9 [6 q3 e; E! W+ s5 E
   
5. Euler(a,b,N,af);
   4 ]: u4 S- S& q& O3 [, e
6. %Taylor(a,b,N,af);, D8 i. O7 F6 p! V% a# h
   
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
8 Q% K) ]: k% j' |0 c4 R' N. q2.N：时间步数。
3.af：步长调整系数。
; e3 k. M; D) k4 ?  r* c; W4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。

! b' F% X: W6 B

6 T" J) A1 r" H; X7 e/ `