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

2744557306

1. a=0;: @9 g/ p/ Y) `6 Z0 N# K
   
2. b=1;; }' b/ G) V8 s( Q; s+ K7 F+ a8 L
   
3. N=40;
   & D6 a\" A$ m- h/ Q' N& ]
4. af=1;
   . |' }1 x- N& Z5 }' C( V
5. Euler(a,b,N,af);
   5 @# g! |) Q) p# }1 t$ F
6. %Taylor(a,b,N,af);\" [) k! N% R$ {9 Q, b8 ]% i
   
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
  t7 Z% I. G4 z: [- P4 r2.N：时间步数。
3.af：步长调整系数。
" g! B! u: ?+ b7 B* s8 j# o$ F, s$ e4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。
& `/ k( V& r; f9 Z/ {8 h2 d

% K" e1 }8 l$ e1 T2 C& ?% ^9 ~

7 I, e5 W" {* Q( ]% E% g0 R. x