[R-br] Cálculo de Volume de Sólido Irregular

[Walmes Zeviani]

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 em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120327/ab0050e5/attachment.html>