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

2744557306

1. a=0;
   9 L( _8 P0 W\" `- j' V, }
2. b=1;
   - y$ ^, O+ s8 K6 z& I: a) d0 I
3. N=40;  {5 K8 W( @6 k8 w# w5 h
   
4. af=1;
     s- i8 M: Z* j
5. Euler(a,b,N,af);
   2 j) ?0 m# J* V
6. %Taylor(a,b,N,af);
   : U! `' ^4 M* _5 ]8 X
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
2.N：时间步数。
3.af：步长调整系数。
" t# b) _; r0 a: S% ^8 g4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。
" z" Y4 k# H5 x) l5 T


- L  v% y" Q" L6 P/ G: }
! R! S" U. g$ |6 e- j. g3 Y