MATLAB爱好者论坛-LabFans.com - 查看单个帖子

wxuesen · 2012-11-13, 00:13

function c=com(k)
global k
syms x;
c=limit('x+6000!/k!/(6000-k+1)!',x,0);

global k
syms y;
s=0;
for k=941:1059
s=s+limit('y+com(k).*(1/6)^k.*(5/6)^(6000-k)',y,0);
end
s
这是我自己编的程序，一运行就死机了，请问是哪里出问题了？
供大家参考，帮忙改进下这个程序。
谢谢了。

2012-11-13, 00:13	#2
wxuesen 初级会员注册日期: 2012-04-15 帖子: 2 声望力: 0	回复: 关于求组合数及阶乘中大数的算法。 function c=com(k) global k syms x; c=limit('x+6000!/k!/(6000-k+1)!',x,0); global k syms y; s=0; for k=941:1059 s=s+limit('y+com(k).(1/6)^k.(5/6)^(6000-k)',y,0); end s 这是我自己编的程序，一运行就死机了，请问是哪里出问题了？供大家参考，帮忙改进下这个程序。谢谢了。