fn cross_product(a: &[f64; 3], b: &[f64; 3]) -> [f64; 3]
Helper function to compute cross product of two 3D vectors.