无穷级数求和循环150次出现NaN

求级数和∑((2n-1)!!/(2n)!!)^2时，n为1到正无穷，循环达到150次次出现NaN是怎么回事，我看了一下级数是个递减数列

这是程序
for x=1:inf
y(x)=(dfac(2x-1)/dfac(2x))^2;
end
sum(y)

其中,双阶乘的程序如下：
function y1=dfac(n)
y1=1;
for i=n:-2:1
y1=y1*i;
end

xiayi

【以下回答由 GPT 生成】

这个问题在循环达到150次时出现NaN是因为在计算 dfac(2 * x) 时的 x 值太大，导致阶乘结果超出了浮点数范围，所以计算结果为无穷大（Infinity）。当我们用无穷大去除以无穷大时，结果将会是 NaN。

解决这个问题的一种方式是改进 dfac 函数，使用对数运算来计算阶乘。这可以避免超出浮点数范围的问题。修改后的代码如下：

import math

def dfac(n):
    if n % 2 == 0:
        start = 2
        step = 2
    else:
        start = 1
        step = 2
    result = 1
    for i in range(start, n+1, step):
        result *= i
    return result

y_values = []
for x in range(1, 151):
    y = (dfac(2 * x - 1) / dfac(2 * x)) ** 2
    y_values.append(y)

sum_of_y = sum(y_values)

这样，通过使用对数运算计算阶乘，可以避免大数运算时超出浮点数范围的问题，可以正常计算出结果。

就是到后面太大了，计算得到分子分母都是无穷大，所以求出来就是NaN，表示无法用matlab中的数学表达