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

2744557306

1. a=0;
   & J& _7 t( x, E\" ]8 g\" f
2. b=1;
   ! U+ o; }1 _/ e0 [
3. N=40;; [: R\" ?; ?: c
   
4. af=1;6 h: N' r* y/ P$ F# _
   
5. Euler(a,b,N,af);
   4 O( p' e' Z6 j/ v/ m# l) g, m
6. %Taylor(a,b,N,af);
   9 U. m) H\" n- q' n2 |  L
7. %RK(a,b,N,af);

复制代码

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


9 E7 u2 p7 U2 G0 O& w