library(lattice)
x <- 0:100
y <- 0:100
z <- 0:100
r <- 30
xyz <- expand.grid(x,y,z)
dados <- cbind(xyz,as.numeric(sqrt(rowSums((xyz-50.5)^2))<=r))
str(dados)
colnames(dados) <- c("x","y","z","var")
cloud(z ~ x * y, data = dados[dados$var==1,], groups=var, screen=list(x = -90, y = 70),
xlim = c(0,100), ylim = c(0,100), zlim = c(0,100))
# qual volume do cubo?
100^3
# qual volume da esfera?
table(dados$var) # esta é uma subestimativa
# quanto do volume do cubo é da esfera?
table(dados$var)/100^3