function test_ode
T = 1; y_0 = [1,2];
[t_vec,y_vec] = ode45(@f,[0,T],y_0);
plot(t_vec,y_vec(:,1),'-r'); hold on;
plot(t_vec,y_vec(:,2),'-b'); hold off;

function dy = f(t,y)
A = [-2,0;0,-5]; dy = A*y;