function [mr] = makerotate2d ( theta )

   mr = [ cos(theta) -sin(theta) 0; sin(theta) cos(theta) 0; 0 0 1 ];