res=load('harm.dat');

figure('Position',[10,10,900,900])

subplot(3,1,1)
plot(res(:,1),res(:,2))
xlabel('aeg(s)');ylabel('koordinaat');
grid on

subplot(3,1,2)
plot(res(:,1),res(:,3))
xlabel('aeg(s)');ylabel('kiirus')
grid on

subplot(3,1,3)
plot(res(:,1),res(:,6))
hold on
plot(res(:,1),res(:,7))
xlabel('aeg(s)');ylabel('energies');
legend ("etot",'etotcorr')
grid on