Error in implemented function

Question

Nik joung 2017-11-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/366977-error-in-implemented-function

评论： Nik joung 2017-11-15

采纳的回答： Michal

Hello,

I have the function written above. It was calculated by WolframAlpha and the result is 0.

I've tried to implement the formula in Matlab with the code below:

 tmp=0;
for s=0:24
    if (48-s >= 24)
        bin1 = nchoosek(48-s,24);
    else
        bin1 = 0;
    end
    if (47-s >= 24)
        bin2 = nchoosek(47-s,24);
    else
        bin2 = 0;
    end
    tmp = tmp + (-1)^s * nchoosek(24,s) * (bin1-bin2);  
end

But my code outputs (variable tmp) the result -2 instead of 0. I am not seeing, where the mistake is.

I would appreciate, if someone could help me.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Michal 2017-11-14

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/366977-error-in-implemented-function#answer_291073

编辑：Michal 2017-11-14

The double precision (default numerical values class in Matlab) is not enough for this kind of computation!!! You need to use symbolic math toolbox, to get multi-precision accuracy. Wolfram alpha (with Mathematica engine) using this approach by default.

Just a simple code using matlab symbolic toolbox: tmp = 0 !!!

tmp=sym(0);
syms s
for s=0:24
    bin1 = nchoosek(sym(48)-s,sym(24));
    bin2 = nchoosek(sym(47)-s,sym(24));
    aux = (-1)^s * nchoosek(sym(24),s) * (bin1-bin2);
    tmp = tmp + aux;
end
tmp