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

2744557306

1. a=0;% T% \8 A# G) h\" l) G; i
   
2. b=1;
   ; V% E$ P% c$ y. w2 Q$ F7 z6 d* `
3. N=40;
   & H% k0 y# B0 t/ Q1 f$ g/ m8 K4 p  A
4. af=1;) W% g% X/ l% W9 X
   
5. Euler(a,b,N,af);
   / C6 ~; Z' R$ a+ V( v; D
6. %Taylor(a,b,N,af);
   + b3 b0 V- \5 i1 e
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
2.N：时间步数。
5 ?; y. q* b6 F% K5 U6 j9 g9 e3.af：步长调整系数。
( n5 Y: b& y! u7 n: Q# }1 G* ]4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。



: D& W# p9 }6 t$ ~# s5 b* p% B
. _: `3 f7 D$ v+ q( k4 ?% d* \- H