function M = rotmat(theta)

M = [cos(theta) -1*sin(theta); sin(theta) cos(theta)];
