pub fn forward_substitution(l: &[f64], b: &[f64]) -> Vec<f64>
Solve a matrix equation of the form Lx=b, where L is a lower triangular matrix. See the Wikipedia page.