n <- 1e6; dadi <- trunc(runif(n, min=1, max=7))+trunc(runif(n, min=1, max=7))
n1=0; for(i in 1:n) if(dadi[i]>7) n1=n1+1; print(n1/n)
#[1] 0.416937
# Per avere l'andamento delle uscite (10, 100, ... 10^6 uscite):
for(k in 1:6) {tot=10^k; n1=0; for(i in 1:tot) if(dadi[i]>7) n1=n1+1; print(n1/tot)}
#[1] 0.4
#[1] 0.45
#[1] 0.399
#[1] 0.4103
#[1] 0.4149
#[1] 0.416937
#
# Per avere una (facile) rappresentazione grafica della distribuzione delle
# uscite carico:
source("http://macosa.dima.unige.it/r.R")
# Quindi:
noClassi=1; Istogramma(dadi, 1.5,12.5, 1)
          
#[i lati della griglia distano il   1 %] 
#  Frequenze e frequenze percentuali: 
#27801, 55540, 82972, 111333, 138900, 166620, 139240, 111049, 83185, 55674, 27686
#2.78,5.55,8.3,11.13,13.89,16.66,13.92,11.1,8.32,5.57,2.77
#  Per altre statistiche usa il comando   altrestat() 
altrestat()
          
#   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
#      2       5       7       7       9      12 
#    I pallini marroni sono 5° e 95° percentile 
#           Il pallino rosso è la media