for n in range(1,600):
Gy[n] = (Gy[n]+Gy[n-1])/2
Gy[n] = Gy[n]*0.50 + Gy[n-1]*0.50
Ax[n] = Ax[n]*0.10 + Ax[n-1]*0.90
Oxz = arcsin(Ax/9.81)
Oxz = Oxz - np.average(Oxz)
Oxz = Oxz*180/pi
Gy = Gy - np.average(Gy)
PostProcessing.identifySpikes(Oxz, t, Os, ts)
ts = ts - ts[0]
PostProcessing.periodDampingInertia(ts, m, g, d, Os, T, En, I)
Ta = np.average(T)
Tmax = np.max(T)
Ena = np.average(En)
Ia = np.average(I)
# Ia = (m*g*d*(Ta**2)*(1-Ena**2))/(2*pi)**2
Imax = np.max(I)
print(Ta, Tmax, Ena, Ia, Imax)
