分からない問題はここに書いてね462

>>258
モンテカルロ法でやってみたら１になったけどあってる？
> # Monte Carlo
> N=1e7
> x=runif(N,-1,0)
> y=runif(N,-4,4)
> r=sqrt(x^2+y^2)
> theta=atan(y/x)
> s=diff(range(x))*diff(range(y))
> s*mean(r<4*(1+cos(theta)))/diff(range(y))
[1] 1