1,2

A(n,k)= (3*n - 3*k + 1)A(n - 1, k - 1) + (3*k - 2)A(n - 1, k); a(n)=diagonal(a(n,k),4).

Clear[A, a, b, n, m, k]; A[n_, 1] := 1; A[n_, n_] := 1; A[n_, k_] := (3*n - 3*k + 1)A[n - 1, k - 1] + (3*k - 2)A[n - 1, k]; a = Table[A[n, k], {n, 20}, {k, n}]; b = Table[a[[n + m - 1]][[n]], {m, 1, 10}, {n, 1, 10}]; Table[b[[4]][[n]], {n, 1, 10}]

Cf. A142458.

Sequence in context: A105987 A051963 A033675 this_sequence A011815 A144929 A163398

Adjacent sequences: A144377 A144378 A144379 this_sequence A144381 A144382 A144383

nonn,uned

Roger L. Bagula and Gary W. Adamson (rlbagulatftn(AT)yahoo.com), Oct 01 2008