set.seed(123)
X2=rnorm(1000,mean=2.0,sd=0.1) # generación de números aleatorios de X
Y2=rnorm(1000,mean=3.0,sd=0.2) # generación de números aleatorios de Y
Z2=data.frame(X2,Y2) # generación de matriz de X,Y
A2=apply(Z2,1,prod) # área de la placa A=XY
mediaA=mean(A2) # media del vector de áreas
varianzaA=var(A2) # varianza del vector de áreas
B2=as.numeric(A2>5.9 & A2<6.1) # generación de variable de 0,1. Con valor de 1 cundo se cumple la condición y cero en otros casos
Pro3c=sum(B2)/1000 # calculo de la probabilidad
hist(A2, las=1, main =" ") # histograma del valor de las ár

mediaA
## [1] 6.023568
varianzaA
## [1] 0.2735994