Add example for solve

Author: termoshtt

examples/solve.rs

@@ -0,0 +1,21 @@
+
+extern crate ndarray;
+extern crate ndarray_linalg;
+
+use ndarray::*;
+use ndarray_linalg::*;
+
+// Solve `Ax=b` for many b with fixed A
+fn factorize() -> Result<(), error::LinalgError> {
+    let a: Array2<f64> = random((3, 3));
+    let f = a.factorize_into()?; // LU factorize A (A is consumed)
+    for _ in 0..10 {
+        let b: Array1<f64> = random(3);
+        let x = f.solve(Transpose::No, b)?; // solve Ax=b using factorized L, U
+    }
+    Ok(())
+}
+
+fn main() {
+    factorize().unwrap();
+}