function [R] = build3DYRotation(theta)
    R = [ cos(theta), 0, sin(theta);...
              0       1,      0    ;...
         -sin(theta), 0, cos(theta)];
 