function T=skewx(a);
% Transformation matrix for a skew in the x direction

T=[1 a 0; 0 1 0; 0 0 1];