21/03/03 19:34:25.71 RNCA+Q6l.net
rm(list=ls())
source('toolmini.R')
library(rgl)
# xyz空間において
# 1≦(1+x^2)(1+2y^2)(1+4z^2)≦8
# を満たす点(x,y,z)全体からなる領域の体積を求めよ。
x=y=z=seq(-3,3,length.out=200)
f <- function(x,y,z){
a=(1+x^2)*(1+2*y^2)*(1+4*z^2)
1<=a & a <= 8
}
gr=expand.grid(x,y,z)
idx=mapply(f,gr[,1],gr[,2],gr[,3])
b=gr[idx,]
plot3d(b,col=4,xlab='x',ylab='y',zlab='z')
nrow(b)/nrow(gr)*6^3
binom.test(nrow(b),nrow(gr))$conf*6^3