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);