function T=skewy(a);
% Transformation matrix for a skew in the y direction

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