<html><body><div style="color:#000; background-color:#fff; font-family:Courier New, courier, monaco, monospace, sans-serif;font-size:10pt"><div id="yiv1236998981"><div><div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); font-size: 10pt; font-family: "arial", "helvetica", sans-serif"><div>Tiago,<br> </div><div><br></div><div>No seu exemplo, o volume seria:</div><div><br></div><div><font face="Courier New, courier, monaco, monospace, sans-serif">x <- y <- z <- 0:100<br>r <- 30<br>xyz <- expand.grid(x,y,z)<br>dados <- cbind(xyz,as.numeric(sqrt(rowSums((xyz-50.5)^2))<=r))<br>colnames(dados) <- c("x","y","z","var")<br><br>mean(dados$var)*100^3 # Aproximação<br>4*pi*30^3/3 # Exato</font></div><div><br></div><div><br></div><div>A idéia aqui é que a proporção de pontos que caem dentro da esfera é aproximadamente a razão entre o volume da esfera e
do cubo. Essa aproximação pode ser melhorada aumentando-se o número de pontos da malha. </div><div><br></div><div>Não tinha visto na sua postagem original, mas se vc já tem
uma variável indicadora com essa informação, o mesmo princípio se aplica.</div><div><br></div><div><br></div><div>Att.,</div><div>Rubem</div><div style="font-size:10pt;"> </div> </div></div></div></div></body></html>