login
Numerator of the polynomial A_l(x) = sum_{d=1..l-1} x^(l-d)/d for index l=2n+1 evaluated at x=9.
1

%I #11 Jun 05 2016 23:34:03

%S 171,27819,11267049,12776837121,1034923809573,922117114354533,

%T 970989321415598469,31460054013865485891,43320494377092775505339,

%U 333351204231728907635493393,27001447542770041518585314553

%N Numerator of the polynomial A_l(x) = sum_{d=1..l-1} x^(l-d)/d for index l=2n+1 evaluated at x=9.

%C For denominators see A145626. For general properties of A_l(x) see A145609.

%t m = 9; aa = {}; Do[k = 0; Do[k = k + m^(2 r + 1 - d)/d, {d, 1, 2 r}]; AppendTo[aa, Numerator[k]], {r, 1, 25}]; aa (* _Artur Jasinski_ *)

%t a[n_,m_]:=Integrate[(m-x^n)/(m-x),{x,0,1}]+(m^n-m)Log[m/(m-1)]

%t Table[9 a[2 n, 9] // Simplify // Numerator, {n,1,25}] (* _Gerry Martens_ , Jun 04 2016 *)

%Y Cf. A145609 - A145640.

%K frac,nonn

%O 1,1

%A _Artur Jasinski_, Oct 14 2008

%E Edited by _R. J. Mathar_, Aug 21 2009