seed<-20000; n<-10; draw<-function(n,generator) { result<-c(); for(i in 1:seed) { result<-c(result,sum(generator(n))) } mean<-mean(result); sd<-sd(result); result<-(result); x=seq(from=min(result),to=max(result),length=100); print(min(result)) y<-dnorm(x,sd=sd(result),mean=mean(result)) hist(col="blue",breaks =20,result,ylim=c(0,max(y)),xlab="Wartość",ylab="Częstość",freq=FALSE,main=paste("Liczebność próby: ",n,sep=" ")); lines(x,y,lwd=2,col="red"); } #Rozkład jednostajny na [0,1] par(mfrow=c(2,3)) draw(1,runif); draw(2,runif); draw(200,runif); #Rozkład wykładniczy l=1; draw(1,rexp); draw(2,rexp); draw(200,rexp);