function M=transgeom(x,y,z);

M=[cos(z) -sin(z) x
    sin(z) cos(z) y
    0       0     1];
end