function ng = ij2ng (N,i,j)
    ng = i + (j-1)*N;
endfunction