#!/usr/bin/python
import numpy as np
import matplotlib.pyplot as plt

res=np.loadtxt("harm.dat")
#******** visual ************************************
fig=plt.figure(figsize=(10,10))
p1=plt.subplot(121)
p1.set_xlabel("t")
p1.set_ylabel("v")
p1.plot(res[:,0],res[:,2],color="red",label="x")
p1.grid(color="black")
p1.legend(loc='upper right')

p2=plt.subplot(122)
p2.set_xlabel("x")
p2.set_ylabel("v")
p2.plot(res[:,1],res[:,2],color="black",label="kiirus")
p2.grid(color="black")
p2.legend(loc='upper right')

plt.show()