Nesse teu CMR eu obteria o volume da esfera com isso 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 Estou certo ou cometendo algum engano? À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================