Create lin_alg_solve.py

Author: lucylow

lin_alg_solve.py

@@ -0,0 +1,16 @@
+from numpy import array
+from scipy.linalg import solve
+
+M = array([[1/5.,1,-1,0,1,1/2.],
+           [1,0,-3,5,-2,-2],
+           [1,1,2,-1,3,-2],
+           [1,-1,-1,1,1./2,-1],
+           [-1,5,-1,2,0,1],
+           [-1,1,-1,0,0,-6]], float)
+
+answer = array([24.1312, 46.2798,-61.8372,31.1466,51.2106,-5.7008])
+
+print solve(M, answer)
+
+
+# [-225.76371598  -88.06783231   -3.09392272  119.13851089  142.05030225 24.41510107]