<div><font face="trebuchet ms,sans-serif">Nesse teu CMR eu obteria o volume da esfera com isso</font></div><div><font face="trebuchet ms,sans-serif"><br></font></div><div><div><font face="'courier new', monospace">library(lattice)</font></div>
<div><font face="'courier new', monospace">x <- 0:100</font></div><div><font face="'courier new', monospace">y <- 0:100</font></div><div><font face="'courier new', monospace">z <- 0:100</font></div>
<div><font face="'courier new', monospace">r <- 30</font></div><div><font face="'courier new', monospace">xyz <- expand.grid(x,y,z)</font></div><div><font face="'courier new', monospace">dados <- cbind(xyz,as.numeric(sqrt(rowSums((xyz-50.5)^2))<=r))</font></div>
<div><font face="'courier new', monospace">str(dados)</font></div><div><font face="'courier new', monospace">colnames(dados) <- c("x","y","z","var")</font></div><div>
<font face="'courier new', monospace">cloud(z ~ x * y, data = dados[dados$var==1,], groups=var, screen=list(x = -90, y = 70),</font></div><div><font face="'courier new', monospace">      xlim = c(0,100), ylim = c(0,100), zlim = c(0,100))</font></div>
<div><font face="'courier new', monospace"><br></font></div><div><font face="'courier new', monospace"># qual volume do cubo?</font></div><div><font face="'courier new', monospace">100^3</font></div>
<div><font face="'courier new', monospace"># qual volume da esfera?</font></div><div><font face="'courier new', monospace">table(dados$var) # esta é uma subestimativa</font></div><div><font face="'courier new', monospace"># quanto do volume do cubo é da esfera?</font></div>
<div><font face="'courier new', monospace">table(dados$var)/100^3</font></div><div style="font-family:'trebuchet ms',sans-serif"><br></div></div><font face="trebuchet ms,sans-serif"><div>Estou certo ou cometendo algum engano?</div>
<div><br></div><div>À disposição.</div><div>Walmes.</div><br clear="all"></font><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">fone: (+55) 41 3361 3573</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>