Re: exercise in double number arithmetic

The last version I posted gives the same results as WolframAlpha to
C(n,k) for n up to 130.
130 65 coef.calc               95067625827960698145584333020095113100 ok
WolframAlpha gives   C(130,65)=95067625827960698145584333020095113100
For n>130, my results are wrong.
Notice that I used M=100 the number of columns of coeffs array, so one
have to take M=130 as max value for M.
And we have
130 90 coef.calc             5334728207404302438396453128740400  ok
WolframAlpha gives C(130,90)=5334728207404302438396453128740400
When we take M=130,
100 100 coef.calc 1  ok   (true)
130 128 coef.calc 8385  ok  (true)
130 129 coef.calc 130  ok (true)
130 130 coef.calc 2  ok (false) this is the extreme case, must be 1.
And after this case begin the wrong results.
Ahmed