function(n.otazok=50,n.ziakov=40000) { # n.otazok ... pocet otazok # n.ziakov ... pocet ziakov vysledky<-rep(0,n.ziakov) for(i in 1:n.ziakov){ ziak<-runif(1) # Aky dobry je ziak: 0...retard, 1...genius body<-0 # Kolko ziak dostane bodov for(j in 1:n.otazok){ obtiaznost<-j/n.otazok #Aka tazka je uloha: 0...1+1=?, 1...dokazte NP!=P if(runif(1)<0.6+0.6*(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 }