clear all;
close all;

[t,y] = ode45(@twomassrhs,[0 100],[0 1 0 0]);

plot(t,y(:,1));
hold on;
plot(t,y(:,3),'r');
%axis([80 100 -150 150]);