import numpy

data = numpy.loadtxt('sun_bz')
# sloupce:
# 1 - vln. delka lam [nm]
# 2 - spektrum slunce f_lam [erg/s/cm2]
# 3 - propustnost XYZ
# 4 - propustnost Johnson

s = 0
for l in range(data[:,0].size):
    f = data[l,1] * data[l,2]
    s = s + f

print(s)

# obrazeni dat
#
# import matplotlib.pyplot as plt
# plt.plot(data[:,0],data[:,1])