euler2rotmat function

Make a rotation (or direction cosine) matrix