Function basis_function 

fn basis_function(n: usize, i: usize, u: f64) -> f64