Determine transformation matrix for planar rotation
function planarRotation input Real n[3] "normalized axis of rotation"; input SI.Angle angle "rotation angle"; output Real S[3, 3] "transformation matrix to rotate a frame_along axis n around angle"; end planarRotation;