fortran66 https://blog.hatena.ne.jp/fortran66/ fortran66のブログ https://fortran66.hatenablog.com/ 至点 www.youtube.com Euler の γ 収束が遅い（～O(1/n) ） Euler の γ ですが、引くlog(n) を定義からずらして log(n+1/2) にすると 定義より早く O(1/24n2) で収束するのが草。 計算による比較 小さい項からの和 program euler_gamma implicit none integer, parameter :: kd = kind(0.0d0) real(kd) :: s integer :: i, n = 10**7 s = 0.0_kd do i = n, 1, -1 s = s + 1.0_kd / i end d… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Ffortran66.hatenablog.com%2Fentry%2F2022%2F06%2F22%2F023902" title="【寝言】夏至 - fortran66のブログ" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://m.media-amazon.com/images/I/41mUrsUo3TL._SL500_.jpg Hatena Blog https://hatena.blog 2022-06-22 02:39:02 rich https://fortran66.hatenablog.com/entry/2022/06/22/023902 1.0 100%