import numpy as np


alap = np.double(1)    # alapharmonikus

sum_sq_tot = alap**2   # total harm. négyzetösszeg (alap+harm)
sum_sq_harm  = 0       # harmonikus négyzetösszeg (csak harm)

for i in range(3,10000,2):
    sum_sq_tot += (alap*np.double(1)/np.double(i))**2
    sum_sq_harm += (alap*np.double(1)/np.double(i))**2

print ('{:3d} {:8.2f} {:8.2f}'.format(i, 100*np.sqrt(sum_sq_harm)/np.sqrt(alap**2), 100*np.sqrt(sum_sq_harm)/np.sqrt(sum_sq_tot)))
