#!/home/misha/venv/bin/python
import math as  m
import numpy as np
import cmath as cm
import matplotlib.pyplot as plt

a=np.loadtxt("dim1_rkf45.dat")
print(np.max(a[:,1]))
plt.subplot(2,2,1)
plt.plot(a[:,0],a[:,1])
plt.legend(["H"])
plt.xlabel("t,sek")
plt.title("H(t),m")
plt.grid()


plt.subplot(2,2,2)
plt.plot(a[:,0],a[:,2])
plt.legend(["V"])
plt.xlabel("t,sek")
plt.title("V(t),m/s")
plt.grid()


plt.subplot(2,2,3)
plt.plot(a[:,0],a[:,3],a[:,0],a[:,4],a[:,0],a[:,5])
plt.legend(["at","aa","atot"])
plt.xlabel("t,sek")
plt.grid()

plt.subplot(2,2,4)
plt.plot(a[:,0],a[:,5])
plt.legend(["Etot"])
plt.xlabel("t,sek")
plt.grid()

plt.show()