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

2744557306

1. a=0;
   8 y- q# m( X  P3 Q( M
2. b=1;  s7 C& a' L2 u/ k% {; f
   
3. N=40;
   . T) l3 m3 U( s9 X& e- T/ |
4. af=1;- T8 G& K, t& j$ a) o% F% Z
   
5. Euler(a,b,N,af);
   ( z# s! S+ A& c5 L
6. %Taylor(a,b,N,af);
   ( H2 b3 c- O+ r( @2 C
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
( O, r: c. I& u# g6 e1 _$ ~0 P2.N：时间步数。
' u# ]3 A6 l# v) t6 G% g& _3.af：步长调整系数。
4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。


" Q/ J2 ~; G# Z; E; c" n3 g