|
|

楼主 |
发表于 2018-4-8 23:06
|
显示全部楼层
本帖最后由 elim 于 2018-4-8 20:59 编辑
回到主题. 下面给出一个无可推诿的估算:
把上述 b(n), (以及 nb(n) = n×b(n)) 的表达式输入 pari/gp, 取 a(1) = ln(1+1/2), 我们有
(1) b(n) 对 a(n) 的拟合表现是惊人的: a(n)=b(n)+o(1/n^5).
(2) 当 n > 10^10,a(n) 的计算对大多数计算机已经属老牛破车,但 b(n) 瞬间即得,且与a(n)的理论值的差别已不足为道。
(3) 计算证实 na(677760) < 2 < na(677761), na(n) 对充分大的 n, na(n) > 2。
(4) 计算及分析证实 na(1842345) 是 na(n) 的最大值。
(5) a(n) 的渐近展开式证实 lim A(n) = 2/3,
(6) a(n) 的渐近展开式证实对a(1)=ln(1+1/2), n < 10^388 时 |A(n) - 2/3| > 0.01, n > 10^389 时 |A(n) -2/3| < 0.01
至此我们用多种方式严格证明了 A(n) = n(na(n)-2)/log(n) 的原始迭代数值计算无法任意逼近其极限.
但通过数学分析得到的 a(n) 乃至 A(n) 的渐近展开可以任意逼近相应的极限,也能任意逼近序列的项. 使得原始数值计算无法实现的操作的理论结果可以被任意逼近。这是数学分析非常了不起的成就。
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有帐号?注册
x
|