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

2744557306

1. a=0;; p; V: ]' p/ Z4 T
   
2. b=1;
   * i+ R& D/ S% h( X1 v' ^
3. N=40;  t6 z! Q9 Q1 C& j3 ~
   
4. af=1;
   9 D8 l9 o# l3 h! T8 A& Y
5. Euler(a,b,N,af);- ?+ w/ Q1 r+ k5 E% M8 W; s7 z
   
6. %Taylor(a,b,N,af);3 D% V2 `- g9 O' n
   
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
+ D! l+ D, @3 E. ^; x1 ]2 l2.N：时间步数。
9 C9 K" D9 _/ P7 m4 s2 k9 M5 @+ R3.af：步长调整系数。
1 K  l: B4 A/ e& S0 @" o9 ~4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。

- S. U# s! ~8 U* T

/ X+ j; ~! y) V" }
: A' n. ?# I; K& y/ h& y