E=10000000; N=0 Q=0 for (i in 1:E) { nombreT <- sample(0:9,6,replace=TRUE) uniqueT = unique(nombreT) L=length(uniqueT) nombre=0 uni=0 for (j in 1:6) { nombre=nombre+10**(j-1)*nombreT[7-j] } for (j in 1:L) { uni=uni+10**(j-1)*uniqueT[L+1-j] } #cat('nT=',nombreT,' n=',nombre,' uT=',uniqueT,' u=',uni,'\n'); if (nombreT[6]%%5==0){N=N+1}; #cat('unique=',unique(nombreT),' nombre=',nombreT,'\n'); if (uni==nombre) { ##cat('Q=',Q,'\n'); Q=Q+1}; } cat('P(n)=',N/E,'\t P(q)=',Q/E,'\n');