OFFSET
3,1
FORMULA
Numerators of sequence a[ 3, n ] in (a[ i, j ])^2 where a[ i, j ] = 1/i if j<=i, 0 if j>i
MATHEMATICA
rows = 30;
M = MatrixPower[Table[If[j <= i, 1/i, 0], {i, 1, rows}, {j, 1, rows}], 2];
T = Table[M[[n]]*LCM @@ Denominator[M[[n]]], {n, 1, rows}];
a[n_] := T[[n, 3]];
Table[a[n], {n, 3, rows}] (* Jean-François Alcover, May 06 2022 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from Sean A. Irvine, Nov 04 2019
STATUS
approved