function(n.otazok=50,n.ziakov=40000,obtiaznost=-0.5,sigma=1.2) 
{
# n.otazok ... pocet otazok
# n.ziakov ... pocet ziakov

vysledky<-rep(0,n.ziakov)
for(i in 1:n.ziakov){
  print(i)
  ziak<-rnorm(1,sd=sigma)  # Aky dobry je ziak: 0...priemerny
  body<-0 # Kolko ziak dostane bodov
  for(j in 1:n.otazok){
     if(runif(1)<exp(ziak-obtiaznost)/(1+exp(ziak-obtiaznost))) body<-body+1
     }
   vysledky[i]<-body
   }

vysledky<-100*vysledky/n.otazok
hist(vysledky,main=mean(vysledky),breaks=2*(0:51)-1,prob=TRUE,col="gray",ylim=c(0,0.025))
curve(dnorm(x,mean=mean(vysledky),sd=sd(vysledky)),from=0,to=100,add=TRUE,col="blue",lwd=2)
vysledky
}