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

2744557306

1. a=0;
   \" B0 a3 {: X' X$ _9 p, d+ U& e
2. b=1;, R3 X4 I: v! L\" s: e
   
3. N=40;+ b4 T. B7 g! ~6 \0 @* Z
   
4. af=1;& S) I* ~( J9 k3 a# ?5 R& h
   
5. Euler(a,b,N,af);
   ) h- Q# B7 V* Y\" L, ~. @, q& \( p
6. %Taylor(a,b,N,af);
   & l% E- K, _% ?; {+ x8 }6 m& b. j
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
2.N：时间步数。
3.af：步长调整系数。
0 L( n4 q* N/ n4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。
, Y1 b% W: F. K0 R

3 R& L. H  n. D

, N8 |+ h7 W  C1 L2 G4 j