となり、円周率(の逆数)を得ることができます。 ※a=l と記述されるのが一般ですが、見難いため a にしています。 今回は、この問題をRで実装し、シミュレートして確認します。 実装 Rコード buffon.needle = function(n = 100, sec = 0.1, l = 1.0, d = 1.0) { x <- y <- xd <- yd <- phi <- ctr <- numeric(0) sum <- 0; PI <- rep(NA, n) layout(matrix(c(1, 2), 2)) for (j in 1:n) { # Draw parallel lines xlim <- c(-0.5 * l, 1.5 * l); ylim <- c(0, 2 * d) plot(1, xlim = xlim, ylim = ylim, type = 'n', xl
![ビュフォンの針をRのanimationで可視化モンテカルロ法 - Qiita](https://cdn-ak-scissors.b.st-hatena.com/image/square/33a5285878deff178886232fae16b535426f69d4/height=288;version=1;width=512/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn.qiita.com%252Fassets%252Fpublic%252Fadvent-calendar-ogp-background-7940cd1c8db80a7ec40711d90f43539e.jpg%3Fixlib%3Drb-4.0.0%26w%3D1200%26mark64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTk3MiZoPTM3OCZ0eHQ9JUUzJTgzJTkzJUUzJTgzJUE1JUUzJTgzJTk1JUUzJTgyJUE5JUUzJTgzJUIzJUUzJTgxJUFFJUU5JTg3JTlEJUUzJTgyJTkyUiVFMyU4MSVBRWFuaW1hdGlvbiVFMyU4MSVBNyVFNSU4RiVBRiVFOCVBNiU5NiVFNSU4QyU5NiVFMyU4MyVBMiVFMyU4MyVCMyVFMyU4MyU4NiVFMyU4MiVBQiVFMyU4MyVBQiVFMyU4MyVBRCVFNiVCMyU5NSZ0eHQtYWxpZ249bGVmdCUyQ3RvcCZ0eHQtY29sb3I9JTIzM0EzQzNDJnR4dC1mb250PUhpcmFnaW5vJTIwU2FucyUyMFc2JnR4dC1zaXplPTU2JnM9NTAwMGJiMDY2NjI4NjU2NzIyZGU5MDBiMGNlYjRkM2Q%26mark-x%3D120%26mark-y%3D96%26blend64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZoPTc2Jnc9OTcyJnR4dD0lNDBrYW5lc2hpbiZ0eHQtY29sb3I9JTIzM0EzQzNDJnR4dC1mb250PUhpcmFnaW5vJTIwU2FucyUyMFc2JnR4dC1zaXplPTM2JnR4dC1hbGlnbj1sZWZ0JTJDdG9wJnM9NjQyMjc4ZDBmOTRkYTEyYzhhZDg0YTQwOTdjMWVlYzU%26blend-x%3D120%26blend-y%3D445%26blend-mode%3Dnormal%26txt64%3DaW4g5qCq5byP5Lya56S-44Ko44Km44Os44Kr%26txt-width%3D972%26txt-clip%3Dend%252Cellipsis%26txt-color%3D%25233A3C3C%26txt-font%3DHiragino%2520Sans%2520W6%26txt-size%3D36%26txt-x%3D134%26txt-y%3D546%26s%3D53ebd9482b9f03cea238fa4962194ae5)