tmp(m,n)={
  my(M=matrix(m,n,i,j,0));
  M[1,1]=1;
  M[1,1]
};