作为更大模拟的一部分,我用 Python 编写了以下函数:
#!/usr/bin/python
counter = 1
while (counter < 10000):
oldpa = .5
t = 1
while (t < counter):
newpa = ((oldpa * t) + 1) / (t + 1)
t = t + 1
oldpa = newpa
counter = counter + 1
print str(counter) + "\t" + str(oldpa)
然后,我开始用 C 重写模拟,以便它运行得更快(同时也给自己一个花时间学习 C 的借口)。这是我的上述函数的 C 版本。
#include <stdio.h>
main()
{
int counter, t;
float oldpa, newpa;
counter = 1;
while ( counter < 10000 )
{
oldpa = .5;
t = 1;
while ( t < counter )
{
newpa = ((oldpa * t) + 1) / (t + 1);
t = t + 1;
oldpa = newpa;
}
counter = counter + 1;
printf("%d\t%f\n", counter, oldpa);
}
}
现在,有趣的事情来了。当我运行 Python 函数时,结果收敛到 0.999950,但是当我运行 C 函数时,它收敛到 0.999883。就我的模拟而言,这种差异实际上可以忽略不计,但我仍然想知道为什么我得到不同的结果