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

2744557306

1. a=0;
   : ]$ ?5 k  P, g# R7 X' N
2. b=1;+ ]- {\" K' u  G9 a; r4 K& ~: u) ~' c
   
3. N=40;/ J6 c& `2 e2 M# t
   
4. af=1;& T% G\" }! M6 [$ A. D
   
5. Euler(a,b,N,af);  {' d; y+ ]& `; P. j$ l- {, V
   
6. %Taylor(a,b,N,af);
   5 p$ o4 ^2 G% _5 q& s+ x
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
! E1 g1 A( t( e' ~6 s, B% R9 K: ?2.N：时间步数。
3.af：步长调整系数。
' R6 }4 S5 V! J1 X( ?9 s4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。



1 m* Y9 a/ ?8 S/ {/ @$ j8 ]% @