clienti <- c(1,2); pcl <- c(0.4,0.6)
durata <- c(50,60); pdu <- c(0.5,0.5)
nprove <- 100000
X <- sample(clienti,size=nprove,replace=TRUE,prob=pcl)
n.1 <- sum(X==1)
S.1 <- sample(durata,size=n.1,replace=TRUE,prob=pdu)
n.2 <- sum(X==2)
out <- sample(durata ,size=n.2*2,replace=TRUE,prob=pdu)
ris <- matrix(out,ncol=2)
S.2 <- apply(ris,1,sum)
S <- c(S.1,S.2)
frequenze(S)
mean(S)
var(S)