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

2744557306

1. a=0;
   ' n: ]3 S3 r0 u' Y9 @4 M0 e
2. b=1;/ k# |% v# _# M
   
3. N=40;
   ' g) Z, I! c  n, [0 }) X: l, ]- M
4. af=1;- i' e8 K! f/ b& }1 j\" ]+ ^1 v4 c2 [
   
5. Euler(a,b,N,af);
   $ T3 e! @' B! m7 h8 i( |7 {3 I
6. %Taylor(a,b,N,af);
   - G0 Q7 R) K, f3 b/ ?: h\" i
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
2.N：时间步数。
  M) U2 u0 }9 B# L: G3.af：步长调整系数。
4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。
$ {2 a3 Y, C, X  I: [9 C

7 a6 `! b5 m' h9 x) z
8 H! H  N* w+ q0 P! [, d5 X6 b: C