with(LinearAlgebra); SumRowT:=proc(n) f := (r, s) -> (s - 1)^(r - 1); M := Matrix(n + 1, f); g := m -> n^(m - 1)/(2*m - 1); K := Matrix(n + 1, 1, g); A := MatrixMatrixMultiply(MatrixInverse(M), K); S := ilcm(seq(denom(A[m, 1]), m = 1 .. n + 1)); end; seq(SumRowT(n),n=0..40);