pub fn rotation_matrix(roll: f64, pitch: f64, yaw: f64) -> Matrix3<f64>
Compute the ZYX rotation matrix from Euler angles (roll, pitch, yaw).