pdl> p $points = rint transpose cat sequence(10),sequence(10)/3

[
 [0 0]
 [1 0]
 [2 1]
 [3 1]
 [4 1]
 [5 2]
 [6 2]
 [7 2]
 [8 3]
 [9 3]
]