[t,y]=ode45(@aquabuoy,[0 60],[0 0 0 0]);

plot(t,y(:,1),'r*-',t,y(:,2),'b*-');
set(gca,'FontSize',20)
legend('Mass 1 motion', 'Mass 2 motion');
xlabel('time');
ylabel('vertical motion of mass center of gravity');