function m=rotationMatrix(theta) m = [cos(theta) -sin(theta); sin(theta) cos(theta)];