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

2744557306

1. a=0;
   ! e2 X7 j) V! s
2. b=1;
   2 f) ~  p* v! d7 L& d2 y1 s' ^
3. N=40;$ ]2 P8 \9 D2 B0 U& \! w+ z6 t
   
4. af=1;- J8 T0 {. E1 I( O( J$ V  k8 X
   
5. Euler(a,b,N,af);9 v4 k( F' ~' ]1 t: l4 A
   
6. %Taylor(a,b,N,af);
   ) l2 r, ?. }: a; q( c. P3 r% `
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
2.N：时间步数。
$ |) D  o  l2 D4 a5 H3.af：步长调整系数。
4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。