pub fn is_on_circumference(
    center: &Point2d,
    point: &Point2d,
    radius: f64
) -> bool