basis_matrix function

Generate 2-d basis in directions i, j in n dimensions (i,j <= n)