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

2744557306

1. a=0;. c1 G7 q5 g1 O\" j4 h5 j$ S0 }) h* e
   
2. b=1;4 i# u7 T* y0 j' h$ N. l0 \  ^
   
3. N=40;! C* f' H8 d' H1 A) r, C$ ^4 V
   
4. af=1;* D; |2 c- \  B5 A: V# J) y# D( ]
   
5. Euler(a,b,N,af);% t0 D% o5 Y3 k/ ]/ L  ]6 s
   
6. %Taylor(a,b,N,af);
   - I/ `3 S1 g  E1 X/ C/ m: i
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
, f/ K8 D4 H% P, V2 u: Q2.N：时间步数。
3 Y  P* X  R* I7 k. x' J7 T- e3.af：步长调整系数。
4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。
4 c& y, t& x& d" [. v/ d( c# M


9 g9 R8 p4 d/ v% v( i: i