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

2744557306

1. a=0;
   8 v2 R$ T, D- s+ z5 S* f
2. b=1;
   9 \  G) g  T. B; J8 M' h0 o
3. N=40;
   2 s: M1 A/ ]\" `1 T
4. af=1;
   . A0 P8 _/ Y! R\" F, I
5. Euler(a,b,N,af);
   # b4 r1 T0 [# `. s
6. %Taylor(a,b,N,af);
   - w\" v% B7 M% f$ {\" B9 s
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
2.N：时间步数。
3.af：步长调整系数。
4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。
- O% U$ g1 v/ v/ J% ]4 }

) k, s9 e+ T1 i/ I
4 B! i: P8 X' f2 x% Q' h& d7 o
6 l" }3 c; N& U" @2 G% ?4 h