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

2744557306

1. a=0;+ h9 p5 f0 W. H& e$ H1 a
   
2. b=1;
   7 w# S0 n\" r: w
3. N=40;
   # i% u& P9 ?6 y& P0 S
4. af=1;5 ]3 N4 ]7 J/ e/ L3 S) T/ C( V
   
5. Euler(a,b,N,af);& l\" w- V% ~. l4 T, a  k
   
6. %Taylor(a,b,N,af);
   ; f3 a2 q4 O8 v, v
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
( Y  ?: E/ _3 H6 }7 F8 S' p  \) C, q2.N：时间步数。
3.af：步长调整系数。
4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。
% y' S7 t$ T8 U6 f* @7 m( c
0 z8 {. k) {9 k, J) ]
/ l; s7 w& [3 b  D4 j& r0 }5 C# @
4 a/ S% f$ j) r0 m8 y