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

2744557306

1. a=0;
   3 L8 }8 f& t$ \  x3 \
2. b=1;
   + f$ Y  S8 l$ D0 d/ @2 X4 j0 `\" h
3. N=40;! L1 r4 ?5 @/ I6 L
   
4. af=1;* v: c: }\" d/ F: `+ b& Q* ^\" V
   
5. Euler(a,b,N,af);
   7 e2 u& ?9 `2 _# U9 E: ~, E3 D* U# q
6. %Taylor(a,b,N,af);
   # A4 f) ]5 F) j( S) M& ^
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
2.N：时间步数。
3.af：步长调整系数。
4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。
5 T( f+ q' j+ I1 \  ?; |6 o) p

# B) y) X5 }% \8 `1 W( a

% u1 I8 J  w7 @% ~$ P- ~