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

##</code><code>##

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

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

##</code><code>##

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]
]

##</code><code>##

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

[
 [32 64  1]
 [35 64  1]
 [38 67  1]
 [41 67  1]
 [44 67  1]
 [47 70  1]
 [50 70  1]
 [53 70  1]
 [56 73  1]
 [59 73  1]
]