from pylab import *
from numpy import *  #je dobre zvyknutis si aj na tento import
#uvedene dva importy v neinteraktivne spustanych .py suboroch robia nieco podobne
# ako ked sa interaktivny ipython spustkik s option --pylab 

R=0.5   # polomer gule
m=2.    # hmotnost gule
xmin=-0.5  # dolna hranica x-suradnice
ymin=-0.5  # dolna hranica y-suradnice
zmin=-0.5  # dolna hranica z-suradnice

rho=m/????    # opravte tak, aby rho bola hustota gule teda hmotnost deleno obnjem gule

d=0.01   # hrana malej LEGO kocky
#velka kocka opisana guli ma hranu 1, takze dokopy je 100*100&*100 LEGO kociek

I=0. # vynulovanie sumatora pre vypocet momentu zotrvacnosti

for i in range(0,100):    #for statement, i bude z mnoziny {0,1,2,3,...,99}
	x=xmin+i*d  # x-suradnica LEGO kocky
	for j in range(0,100): 
		y=????   # doplnte:   y-suradnica LEGO kocky
		for k in range(0,100):
			z=????  # doplnte:   z-suradnica LEGO kocky
			#teraz budeme testovat, ci lego kocka (x,y,z) lezi vnutri gule
			if (x**2+????)>R**2:    #doplnte a opravte aj !!!! CHYBU !!!!!
				# ked sme tu, znamena to, ze LEGO kocka (x,y,z)je vnutri gule
				I=I+rho*???? #doplnte, aby tu bol prispevok jednej LEGO kocky k momentu zotrvacnosti
print(I)
print('teoreticka hodnota je')
print(2.*m*R**2/5.)				
