res=load('harm_math.dat');

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

subplot(2,1,1)
plot(res(:,1),res(:,2)/(2*pi))
xlabel('aeg(s)');ylabel('nurk');
grid on

subplot(2,1,2)
plot(res(:,2)/(2*pi),res(:,3))
xlabel('nurk(s)');ylabel('kiirus')
grid on

