Value a(6) in doubt for a = A101877

[Maximilian Hasler]

your maple calculation should give a fraction and not a decimal.
(do the sum using fractions [takes 4 ms in PARI] and not decimals ;
or, if using decimals, you should do the summation in an intelligent
manner (at least start with the smallest ones, or better, reduce it
pairwise or so))
I cannot imagine the fraction one gets is not correct, and you can
check that the conversion to decmal I gave, is correct:

(09:28) gp > 10^frac(log(denominator(%69))/log(10))
%74 = 2.7500177025316476300406310946451800750
(09:28) gp > 10^frac(log(numerator(%69))/log(10))
%75 = 1.6560046469767282799471827882184497196
(09:28) gp > 10*%/%74
%76 = 6.0217963159008818116821418831900358175


On Jan 15, 2008 9:02 AM, Rainer Rosenthal <r.rosenthal at web.de> wrote:
> Maximilian Hasler wrote:
>
> [ a certain PARI/GP calculation ]
>
> I cannot confirm your PAR/GP result. I copy-pasted your list into Maple
> and got
>
> ReciprocalSum := S -> add(1/S[i],i=1..nops(S));
> ReciprocalSum(
>            [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18,
>             19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,
>             36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52,
>             53, 54, 55, 56, 57, 58, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70,
>             72, 74, 75, 76, 77, 78, 81, 82, 84, 85, 86, 87, 88, 90, 91, 92, 93,
>             94, 95, 96, 98, 99, 100, 102, 104, 105, 106, 108, 110, 111, 114, 115,
>             116, 117, 119, 120, 121, 122, 123, 124, 126, 128, 129, 130, 132, 133,
>             134, 136, 138, 140, 141, 143, 145, 147, 148, 150, 152, 153, 155, 156,
>             159, 161, 162, 164, 168, 170, 171, 174, 177, 180, 182, 183, 184, 185,
>             186, 187, 188, 190, 192, 196, 201, 203, 204, 205, 207, 209, 210, 212,
>             215, 217, 221, 222, 228, 230, 232, 238, 242, 245, 246, 247, 248, 250,
>             252, 253, 255, 258, 259, 261, 266, 268, 272, 276, 280, 282, 285, 287,
>             290, 294, 295, 299, 301, 304, 305, 306, 310, 315, 319, 322, 323, 324,
>             328, 329, 341, 344, 345, 348, 354, 360, 363, 368, 370, 372, 375, 376,
>             377, 380, 384, 387, 402, 406, 407, 413, 414, 424, 430, 434, 435, 444,
>             451, 460, 465, 469]
> );
>
> with result 5.768... which is less than 6.
>
> (Thanks for taking back the "achieve a(n)" misinterpretation. So I was able to
> cancel my related response. But please use the a = A101877 in the subject line,
> which I forgot in my first posting.)
>
> Cheers,
>
> Rainer Rosenthal
> r.rosenthal at web.de
>
>
>
>
>



-- 
Maximilian F. Hasler (Maximilian.Hasler(AT)gmail.com)