|  1 |    initial version    |  
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
 
                
                Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.