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

2744557306

1. a=0;) k\" q/ B3 W. q2 k0 W3 u3 O
   
2. b=1;; v\" D9 f, q0 {- a! j& E8 N7 z
   
3. N=40;' _: j1 C6 B+ I+ l) ^
   
4. af=1;\" l  g; z, w/ _. ]! p
   
5. Euler(a,b,N,af);
   : F0 J/ C& W# v% T% b  N
6. %Taylor(a,b,N,af);# p& b: V1 T( u  l1 Q- q' Z' X
   
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
2.N：时间步数。
3.af：步长调整系数。
4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。
$ j8 z- h6 |$ P& s# j


5 G7 J4 |  f0 e: W8 c, Y
- J* c- f4 F5 i