Function coplanar_3d_epsilon 

pub fn coplanar_3d_epsilon<S>(
    a: Point3<S>,
    b: Point3<S>,
    c: Point3<S>,
    d: Point3<S>,
    epsilon: S,
) -> bool
where
    S: OrderedField + AbsDiffEq<Epsilon = S>,