Function osculating_circle 

pub fn osculating_circle(
    pos: Vec3,
    tangent: Vec3,
    normal: Vec3,
    curvature: f32,
) -> (Vec3, f32)

Compute the osculating circle at a point: returns (center, radius)