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

2744557306

1. a=0;
   \" D5 D( R9 j- J0 n$ ~( L3 E
2. b=1;/ ^0 g6 W5 [* E2 y8 a
   
3. N=40;5 u# p9 s3 a+ n$ l
   
4. af=1;
   - o6 X+ s- E/ |- o4 Q0 i
5. Euler(a,b,N,af);! o/ |: g. Y! d. ~# ^
   
6. %Taylor(a,b,N,af);
   ! H$ M: G: ]* T1 |/ ^\" }! m6 B
7. %RK(a,b,N,af);

复制代码

1.a 和 b：定义了求解的时间范围 ([a, b])。
3 a" U/ v/ k6 u* b9 }2.N：时间步数。
3.af：步长调整系数。
# A: u% l1 _& i3 r7 h. h% o! l4.Euler(a, b, N, af)：调用 Euler 方法进行求解。其他两个方法 Taylor 和 RK 被注释掉了。


2 e6 D% V0 H% h
% x9 V' l1 U' Z" [) a( ]9 k  f