function T=translate(tx,ty)
% translate

T=[1 0 tx; 0 1 ty; 0 0 1];
