% myode45_1.m
% Using ode45 to solve
% dy/dt = -y - 5*exp(-t)*sin(5*t), y(0) = 1
myf = @(t,y) -y-5*exp(-t)*sin(5*t);
y_exa = @(t) exp(-t).*cos(5*t);
tspan1 = [0,3]; y0 = 1;
[t1,y1] = ode45(myf, tspan1, y0);
ye = y_exa(t1);
tspan2 = linspace(0,3,100);
[t2,y2] = ode45(myf, tspan2, y0);
plot(t1,ye, t1, y1, '*--', t2, y2,'o--');
xlabel t, ylabel y(t)
err_max1 = max( abs(y1-ye ) )
err_max2 = max( abs(y2-y_exa(t2) ) )
err_max1 =
2.8991e-004
err_max2 =
2.8741e-004
|
| LAST UPDATE: 2009.11.19 - 15:01 |
|
메일보내기
이름검색
Fox and Rabbit
ode_2nd_2
