Column numbers of elements in a triangular representation of a symmetric matrix
jtriag(blocksize, istart = 1, diag = "lower")
Useful links