Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- ##################################################
- # Typ kodowania UTF8
- # Realizator: Kretek Kacper
- # Tytu: Symulacyjne badanie wielokrotnych sum wielko±ci losowych
- ##################################################
- N=10000 # Ilosc rzutow
- N2=100000 # Ilsc rzutow
- K=1:6 # Definiujemy zwykla szescio scienna kostke
- # Ad I.1
- rzut1=sample(K,N,replace=T)
- hist(rzut1,main = "rzut jedna kostka N", xlab="os X", ylab="os Y" )
- # Ad I.2
- s2k=function(){sample(K,1)+sample(K,1)}
- rzut2=sapply(1:N, function(j){s2k()})
- hist(rzut2,main = "rzut dwoma kostkami N", xlab="os X", ylab="os Y")
- # Ad I.3
- s3k=function(){sample(K,1)+sample(K,1)+sample(K,1)}
- rzut3=sapply(1:N, function(j){s3k()})
- hist(rzut3,main = "rzut trzema kostkami N", xlab="os X", ylab="os Y")
- # Ad I.4
- s4k=function(){sample(K,1)+sample(K,1)+sample(K,1)+sample(K,1)}
- rzut4=sapply(1:N, function(j){s4k()})
- hist(rzut4,main = "rzut czterema kostkami N", xlab="os X", ylab="os Y")
- # Ad I.5
- s5k=function(){sample(K,1)+sample(K,1)+sample(K,1)+sample(K,1)+sample(K,1)}
- rzut5=sapply(1:N, function(j){s5k()})
- hist(rzut5,main = "rzut piecioma kostkami N", xlab="os X", ylab="os Y")
- # Ad II
- # 1Kostka
- rzut1=sample(K,N2,replace=T)
- hist(rzut1,main = "rzut jedna kostka N2", xlab="os X", ylab="os Y" )
- # 2Kostka
- s2k=function(){sample(K,1)+sample(K,1)}
- rzut2=sapply(1:N2, function(j){s2k()})
- hist(rzut2,main = "rzut dwoma kostkami N2", xlab="os X", ylab="os Y")
- # 3Kostka
- s3k=function(){sample(K,1)+sample(K,1)+sample(K,1)}
- rzut3=sapply(1:N2, function(j){s3k()})
- hist(rzut3,main = "rzut trzema kostkami N2", xlab="os X", ylab="os Y")
- # 4Kostka
- s4k=function(){sample(K,1)+sample(K,1)+sample(K,1)+sample(K,1)}
- rzut4=sapply(1:N2, function(j){s4k()})
- hist(rzut4,main = "rzut czterema kostkami N2", xlab="os X", ylab="os Y")
- # 5Kostka
- s5k=function(){sample(K,1)+sample(K,1)+sample(K,1)+sample(K,1)+sample(K,1)}
- rzut5=sapply(1:N2, function(j){s5k()})
- hist(rzut5,main = "rzut piecioma kostkami N2", xlab="os X", ylab="os Y")
- # Ad III
- ZK=c(0,2,2,4,4,6) # Definiujemy zepsuta kostke
- # 1Kostka
- rzut1=sample(ZK,N,replace=T)
- hist(rzut1,main = "rzut jedna kostka ZK", xlab="os X", ylab="os Y" )
- # 2Kostka
- s2k=function(){sample(ZK,1)+sample(ZK,1)}
- rzut2=sapply(1:N, function(j){s2k()})
- hist(rzut2,main = "rzut dwoma kostkami ZK", xlab="os X", ylab="os Y")
- # 3Kostka
- s3k=function(){sample(ZK,1)+sample(ZK,1)+sample(ZK,1)}
- rzut3=sapply(1:N, function(j){s3k()})
- hist(rzut3,main = "rzut trzema kostkami ZK", xlab="os X", ylab="os Y")
- # 4Kostka
- s4k=function(){sample(ZK,1)+sample(ZK,1)+sample(ZK,1)+sample(ZK,1)}
- rzut4=sapply(1:N, function(j){s4k()})
- hist(rzut4,main = "rzut czterema kostkami ZK", xlab="os X", ylab="os Y")
- # 5Kostka
- s5k=function(){sample(ZK,1)+sample(ZK,1)+sample(ZK,1)+sample(ZK,1)+sample(ZK,1)}
- rzut5=sapply(1:N, function(j){s5k()})
- hist(rzut5,main = "rzut piecioma kostkami ZK", xlab="os X", ylab="os Y")
- # 1Kostka
- rzut1=sample(ZK,N2,replace=T)
- hist(rzut1,main = "rzut jedna kostka ZK N2", xlab="os X", ylab="os Y" )
- # 2Kostka
- s2k=function(){sample(ZK,1)+sample(ZK,1)}
- rzut2=sapply(1:N2, function(j){s2k()})
- hist(rzut2,main = "rzut dwoma kostkami ZK N2", xlab="os X", ylab="os Y")
- # 3Kostka
- s3k=function(){sample(ZK,1)+sample(ZK,1)+sample(ZK,1)}
- rzut3=sapply(1:N2, function(j){s3k()})
- hist(rzut3,main = "rzut trzema kostkami ZK N2", xlab="os X", ylab="os Y")
- # 4Kostka
- s4k=function(){sample(ZK,1)+sample(ZK,1)+sample(ZK,1)+sample(ZK,1)}
- rzut4=sapply(1:N2, function(j){s4k()})
- hist(rzut4,main = "rzut czterema kostkami ZK N2", xlab="os X", ylab="os Y")
- # 5Kostka
- s5k=function(){sample(ZK,1)+sample(ZK,1)+sample(ZK,1)+sample(ZK,1)+sample(ZK,1)}
- rzut5=sapply(1:N2, function(j){s5k()})
- hist(rzut5,main = "rzut piecioma kostkami ZK N2", xlab="os X", ylab="os Y")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement