再论自然数全加和 - 欧拉伽马常数

此处具体算一下欧拉伽马常数是怎么得到的。

在,

中, 之间的1个单位,被平分成 份,如果周期是 ,但这个平分的每一份都取边界值的倒数为结果。而在,

过程中,每一个单位也被平分成 份,周期也是 ,但是这平分的每一份则取对应数值的倒数为结果。所以,对于积分形式,

来说,整个区间被分成,

份,每一份的单位是

而对于求和形式,

来说,整个区间被分成,

份,每一份的单位都是1。相比较而言,积分形式的对应的求和份数更多,细节更细,

在每个 的区间,积分形式为,

将内部精度提升到

求和形式为,

的每个部分都提升精度到

以最高精度计算面积的差值是,

每一个微小面积,

以最高精度计算的总面积差为,

这里的内层分的更细,外层没有那么细。回到微小面积,

以最高精度计算的总面积差为,

这个公式用程序验证,得到很好的结果。

cs 复制代码
    static double GetGamma(long S = 10000)
    {
        var q = S * S;
        var s = 0.0;
        for (var n = 1L; n <= S; n++)
        {
            for (var k = 1L; k <= q; k++)
            {
                s += k / (1.0 * n * n * q + 1.0 * n * k);
            }
        }
        return s / q;
}
r=0.5771647389914232
γ=0.5772156649015329
(γ-r)/γ=(0.5772156649015329-0.5771647389914232)/0.5772156649015329=8.82268×10^(-5)
 
相关推荐
_清歌14 小时前
DSpark 深度解读:DeepSeek-V4 如何用「半自回归」把推理速度提升 85%
算法
统计实现局14 小时前
SVD 的三步走:双对角化、Givens 收敛、排序
算法
躬行见万象14 小时前
《VLA 系列》UniLab 强化训练 | G1 机器人 |复现
算法
统计实现局14 小时前
对称不定分解(Bunch-Kaufman):为什么 Cholesky 不够用
算法
统计实现局14 小时前
dqrsl 拆解:拿着 QR 结果能算出哪 5 种东西
算法
统计实现局14 小时前
为什么 Cholesky 求逆比 Gauss-Jordan 快一倍——行列式溢出防护详
算法
To_OC1 天前
LC 994 腐烂的橘子:人人都说是 BFS 入门题,我却写了三遍才过
javascript·算法·leetcode
金銀銅鐵1 天前
[Python] 扩展欧几里得算法
python·数学·算法
To_OC1 天前
LC 200 岛屿数量:经典 DFS 入门题,我第一次写居然连方向都搞错了
javascript·算法·leetcode
金銀銅鐵1 天前
n^5 和 n 的个位数是否总相等?
python·数学