[go: up one dir, main page]

login
Column 1 of A102054, the matrix inverse of A060083 (Euler polynomials).
5

%I #9 Mar 30 2012 18:51:51

%S 1,2,1,4,-13,142,-1931,36296,-893273,27927346,-1081725559,50861556172,

%T -2854289486309,188475382997654,-14467150771771043,

%U 1277417937676246672,-128570745743431055281,14632875988040732946106,-1869882665740777942166543,266593648798424693540514836

%N Column 1 of A102054, the matrix inverse of A060083 (Euler polynomials).

%C 1-a(n+1) equals the n-th partial sum of the Genocchi numbers (A001469).

%F a(n) = 1 - Sum_{k=1, n} A001469(k) for n>0, with a(0)=1.

%F This sequence's twin numbers are given in A133135. - _Paul Curtz_, Aug 07 2008

%o (PARI) {a(n)=local(M=matrix(n+2,n+2));M[1,1]=1;if(n>0,M[2,1]=1;M[2,2]=1); for(r=3,n+2, for(c=1,r,M[r,c]=if(c==1,M[r-1,1], if(c==r,1,M[r,c]=M[r-1,c]-((matrix(r-1,r-1,i,j,M[i,j]))^-1)[r-1,c-1])))); return(if(n==0,1,M[n+2,2]))}

%Y Cf. A001469, A060083, A102054, A102056.

%K sign

%O 0,2

%A _Paul D. Hanna_, Dec 28 2004