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 spusti s option --pylab 

R=0.5   # polomer gule
m=2.    # hmotnost gule

rmin=0.  # dolna hranica r
thetamin=0.  # dolna hranica theta
phimin=0.  # dolna hranica phi

rho=m/((4./3.)*pi*R**3)    # hustota gule

dr=R/100.   #krok v r
dtheta=pi/100.  #krok v theta
dphi=2.*pi/100.   # krok v psi


I=0. # vynulovanie sumatora
for i in range(0,100):    #for statement, i bude z mnoziny {0,1,2,3,...,99}
	r=rmin+i*dr  # r suradnica sferickej LEGO kocky
	for j in range(0,100): 
		theta =thetamin+j*dtheta   # theta-suradnica sferickej LEGO kocky
		for k in range(0,100):
			phi=phimin+k*dphi  # phi-suradnica sferickej LEGO kocky
			I=I+rho*(r**2)*sin(theta)*dtheta*dphi*dr*??????  # tu doplnte chybajuci kvadrat vzdualenosti od osi z
			#v predchadzajucom riadku si vsimnite ako je napisany objem infinitezimalnej sferickej LEGO kocky
print(I)
print('teoreticka hodnota je')
print(2.*m*R**2/5.)