Select Code to Download

                          | a b 0 |
[ x' y' 1 ] = [ x y 1 ] x | c d 0 |
                          | e f 1 |

pdl> p $m = pdl '3 0 0 ; 0 3 0 ; 32 64 1'

...
 [ 0  3  0]
 [32 64  1]
]

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

...
 [8 3]
 [9 3]
]

pdl> p $points->append(1) x $m

...
 [56 73  1]
 [59 73  1]
]