// Solve `Ax=b` for tridiagonal matrix
Solve `Ax=b` for tridiagonal matrix
Solve ` Ax = b ` for tridiagonal matrix
fn solve() -> Result<(), error::LinalgError> {
    let mut a: Array2<f64> = random((3, 3));
    let b: Array1<f64> = random(3);
    a[[0, 2]] = 0.0;
    a[[2, 0]] = 0.0;
    let _x = a.solve_tridiagonal(&b)?;
    Ok(())
}
fn solve ( ) -> Result < ( ) , error :: LinalgError > { let mut a : Array2 < f64 > = random ( ( 3 , 3 ) ) ; let b : Array1 < f64 > = random ( 3 ) ; a [ [ 0 , 2 ] ] = 0.0 ; a [ [ 2 , 0 ] ] = 0.0 ; let _x = a . solve_tridiagonal ( & b ) ? ; Ok ( ( ) ) }
Solve `Ax=b` for tridiagonal matrix
Solve ` Ax = b ` for tridiagonal matrix