TA的每日心情 | 奋斗 9 小时前 |
---|
签到天数: 2384 天 [LV.Master]伴坛终老
|
clc;
& P+ z% o( U7 p# Qclear;& m0 z6 t. K* o* s
y(1)=1;
* S' [3 ]- g5 ]8 S! G0 f! @1 V9 Kh=0.1;5 y, ^0 |5 s( m0 r' x* X
ta=0;
$ @7 g" T% B3 g; g* Xtb=1;& o* k! G6 a. O' _1 X r1 D' C
n=round((tb-ta)/h);# S5 B; n$ d, Y2 P# A; R
t=ta;
* D0 y2 ?& j" D5 q& ~! @for k=1:n
6 R+ a. \3 K- A9 t; I/ f y(k+1)=y(k)+h*(-y(k));
0 f/ j% a5 \, O! u8 N/ ~& f/ l t=[t,t(k)+h];
4 c/ J+ {6 Z r& k: wend7 g# @* t1 Q* n& |/ N. Y
y
2 U" }$ L3 d1 V0 eplot(t,y); |
|