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 

U0 = 20  # napatie baterie
R=10000  # odpor v ohmoch
C=0.0001 # kondenzator vo faradoch (100 mikrofarad)
t=empty(1001) # nenaplnene array casovych momentov indexovanych cislami od 0 po 1000
Q=empty(1001) # nenaplnene array hodnot naboja indexovanych cislami od 0 po 1000
dt=5./1000  # maly krok v case, 5 sekund rozdelene na 1000 intervalov
t[0]=0  # pociatocny cas
Q[0]=0  # pociatocny naboj

for i in range(0,1000):    #for statement, i bude z mnoziny {0,1,2,3,...,999}
	t[i+1]=t[i]+dt  #toto je posun v case o dt, novy casovy okamih ma index i+1
	Q[i+1]= Q[i]+dt*(U0/R-Q[i]/C/R)   #sem vlozte vypocet noveho naboja  po casaovom kroku dt, novy naboj bvude mat index i+1
					#novy naboj vyjadrite plolmocou stareho naboja platneho  pre cas s indexom i teda stary naboj ma index i
plot(t,Q)     # vykresli graf naboja v zavislosti na case
show()