First time here? Check out the FAQ!
answered 2022-07-28 19:36:04 +0100
This is not shiny but works
var('n p') n = 2 p = 5 vars1 = [var('a_%d%d' %(i,j)) for i in [1..n] for j in [1..p]] matrix(SR, n, p, vars1)
adapt for n, p > 9