lin2009 发表于 2011-1-17 16:56:17

结果应该为0的求极限运算的问题

本帖最后由 lin2009 于 2011-1-17 16:58 编辑

limit(Q(n), n = infinity) = 0,但只是第一次运算能正确显示为0,后面运行的结果都不等于0。大家看看这是这么回事?
d := proc (n) options operator, arrow; (5/3)*(5*sqrt(6)*(50-20*sqrt(6))^n-12*(50-20*sqrt(6))^n-5*sqrt(6)*(50+20*sqrt(6))^n-12*(50+20*sqrt(6))^n)*sqrt(6)/((50-20*sqrt(6))^n-(50+20*sqrt(6))^n) end proc;

Q := proc (n) options operator, arrow; 10^n*(product(1/d(k), k = 1 .. n)) end proc;

evalf(Q(1120));

limit(Q(n), n = infinity);
应为0,但显示如下结果:(第一次运行时等于0)
signum((1/10)^n)*signum(6^((1/2)*n))*signum(product(((50-20*sqrt(6))^k-(50+20*sqrt(6))^k)/(5*sqrt(6)*(50-20*sqrt(6))^k-12*(50-20*sqrt(6))^k-5*sqrt(6)*(50+20*sqrt(6))^k-12*(50+20*sqrt(6))^k), k = 1 .. n))*infinity
页: [1]
查看完整版本: 结果应该为0的求极限运算的问题