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

2744557306

1. a=0;
   ( I- c! Y3 c3 ~: D, ^
2. b=1;/ `0 ~/ F8 p3 R: I
   
3. N=40;; Y( e* p) _1 {$ _* ^2 [
   
4. af=1;
   7 ^8 R% l: }- c, U/ F# K. V# S
5. Euler(a,b,N,af);
   3 r, {5 K: [* P/ M\" {, e4 v* b
6. %Taylor(a,b,N,af);* a; c2 I1 m6 }, W' @& I\" N
   
7. %RK(a,b,N,af);

复制代码

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

* K0 e$ r, ?4 G% S" g7 Q! Z
% ]1 q$ g8 v' Z
3 b3 Y9 n' j8 k3 Q/ |$ V- ^. j