2,2

restart; nmax:=17: jn:=nmax: im:=nmax: Omega(0):=1: for n from 1 to nmax do for j from 1 to jn do cfn1(1, j):=1 end do: for i from 2 to im do cfn1(i, 1):=0 end do: for j from 2 to jn do for i from 2 to im do cfn1(i, j):=cfn1(i-1, j-1)*(j-1)^2+cfn1(i, j-1) end do end do: Omega(n):= (sum((-1)^(k+n+1)*(bernoulli(2*k)/(2*k))*cfn1(n-k+1, n), k=1..n))/(2*n-1)! end do: for n from 1 to nmax do d(n):=(Omega(n)*2^(2*n-1)) end do: for n from 2 to nmax do Zc(n-1):= d(n-1)*2/((2*n-1)*(n-1)) end do: c(1):=denom(Zc(1)): for n from 1 to nmax-1 do c(n+1):= lcm(c(n)*(n+1)*(2*n+3)/2, denom(Zc(n+1))); p(n+1):=c(n) end do: y(1):=Zc(1): for n from 1 to nmax-2 do y(n+1):= Zc(n+1)-((2*n+2)/(2*n+3))*y(n) end do: for n from 2 to nmax do ZETA(n, 1):= p(n)*y(n-1) end do: a:=n-> ZETA(n, 1): seq(a(n), n=2..nmax);

A160474 is the Zeta triangle.

Sequence in context: A015271 A099397 A093251 this_sequence A135531 A038777 A022079

Adjacent sequences: A160472 A160473 A160474 this_sequence A160476 A160477 A160478

easy,sign

Johannes W. Meijer (meijgia(AT)hotmail.com), May 24 2009