Function fixed_point 

pub fn fixed_point(
    g: impl Fn(f64) -> f64,
    x: f64,
    tol: f64,
    max_iter: usize,
) -> Option<f64>

Fixed-point iteration: x_{n+1} = g(x_n).