function R = rot_matrix(theta)

    R = [ cos(theta) sin(theta); 
        -sin(theta) cos(theta)];

end
