c/c++语言开发共享如何在C中实现连续分数的自然对数？

中央循环是混乱的。 返工。 也不需要递归。 只需先计算最深的术语，然后逐步解决问题。

 double contFragLog(double z, int n) { double zz = z*z; double cf = 1.0; // Important this is not 0 for (int i = n; i >= 1; i--) { cf = (2*i -1) - i*i*zz/cf; } return 2*z/cf; } void testln(double z) { double y = log((1+z)/(1-z)); double y2 = contFragLog(z, 8); printf("%e %e %en", z, y, y2); } int main() { testln(0.2); testln(0.5); testln(0.8); return 0; } 

产量

 2.000000e-01 4.054651e-01 4.054651e-01 5.000000e-01 1.098612e+00 1.098612e+00 8.000000e-01 2.197225e+00 2.196987e+00 

[编辑]

根据@MicroVirus的提示，我发现double cf = 1.88*n - 0.95; 比double cf = 1.0;更好地工作double cf = 1.0; 。 使用更多的术语时，使用的值会产生较小的差异，但是良好的初始cf需要较少的术语才能获得良好的答案，尤其是对于|z| 接近0.5。 当我研究0 < z <= 0.5可以在这里做更多的工作。 @MicroVirus对2*n+1建议可能接近我的建议，因为n是一个n 。

这是基于反向计算并注意到CF[n]的值随着n增加而增加。 我很惊讶“种子”值似乎不是一个很好的整数方程式。

这是使用递归的问题的解决方案（如果有人感兴趣）：

 #include  #include  /* `i` is the iteration of the recursion and `n` is just for testing when we should end. 'zz' is z^2 */ double recursion (double zz, int i, int n) { if (!n) return 1; return 2 * i - 1 - i * i * zz / recursion (zz, i + 1, --n); } double contFragLog (double z, int n) { return 2 * z / recursion (z * z, 1, n); } void testln(double z) { double y = log((1+z)/(1-z)); double y2 = contFragLog(z, 8); printf("%e %e %en", z, y, y2); } int main() { testln(0.2); testln(0.5); testln(0.8); return 0; } 

输出与上面的解决方案相同：

 2.000000e-01 4.054651e-01 4.054651e-01 5.000000e-01 1.098612e+00 1.098612e+00 8.000000e-01 2.197225e+00 2.196987e+00