Fibonacci even numbers. Code doesn't work for large numbers

Question

Maria K 2020-7-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/563372-fibonacci-even-numbers-code-doesn-t-work-for-large-numbers

评论： Maria K 2020-7-17

I have writen a code to calculate for every d numbers all the fibonacci numbers that are even, but it doesn't seem to work for large d such as 100 or 150. Can you tell me what i did wrong? My code is this

function y = evenFibo(d)

%y = x;

fibonacci(1) = 0; fibonacci(2) = 1;

for i = 3 : d +1

fibonacci(i) = fibonacci(i-1) + fibonacci(i-2) ;

end

fibonacci = nonzeros(fibonacci);

even = [];

for i = 1 : d

if rem (fibonacci(i),2) == 0 == 1 %einai artios

even = [even, fibonacci(i)];

fibonacci(i)

end

y= numel(even);

end

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

请先登录，再进行评论。

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

Answer 1

John D'Errico 2020-7-11

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/563372-fibonacci-even-numbers-code-doesn-t-work-for-large-numbers#answer_464417

编辑：John D'Errico 2020-7-11

Why are you surprised?

log2(fibonacci(100))
ans =
       68.263

You cannot represent an integer exactly in double precision if it is greater than 2^53-1. That appears to happen around n == 79.

log2(fibonacci(79))
ans =
       53.684

So while there is nothing overtly wrong with your code, what is wrong is your assumption that you can use double preciion arithmetic to compute the result for large index.

If you seriously need to compute that result for larger indices, then you need to use a tool like the symbolic toolbox, or perhaps my own VPI toolbox, which can handle numbers with more than roughly 16 decimal digits. Or, you could write your own code, which might be less difficult than you may think. (Of course, I don't know how hard you think that is. This is indeed a conundrum.)

Finally, many such questions really don't ask to compute an actual sum, but they often ask for something like a remainder with respect to some given integer modulus. In fact, that is very doable in double precision arithmetic. In the case of your code, it seems you are actually computing and saving the numbers themselves, so it would not help here.