% AXIFY : script to label the axes of the harmonic oscillator plot

axis([0 10*pi -.5 2]); 
legend('position x', 'velocity xdot'); 
xlabel('time');
title('Harmonic Oscillator'); 
shg
