OFFSET
1,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..100
EXAMPLE
P_3(x) = 18 - 27*x + 10*x^2; a(3) = P_3(3+1) = 70.
MATHEMATICA
P[n_][x_] = Sum[a[i]*x^i, {i, 0, n - 1}]; ecu[n_] := Table[P[n][i] == i^i, {i, 1, n}]; PP[n_][x_] := P[n][x] /. Solve[ecu[n]][[1]]; Table[PP[i][i + 1], {i, 1, 22}]
a[n_] := InterpolatingPolynomial[Table[{i, i^i}, {i, n}], n+1]; Array[a, 20] (* Giovanni Resta, Jun 18 2013 *)
PROG
(PARI) a(n)=subst(polinterpolate(vector(n, i, i^i)), 'x, n+1) \\ Charles R Greathouse IV, Nov 19 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
José María Grau Ribas, Jun 18 2013
STATUS
approved