Ex 10.1
Code
import numpy as np
import matplotlib.pyplot as plt
import scipy.linalg as lin
m=3
n=3
A=np.random.random((m,n))
b=np.random.random((m,1))
x,res,rnk,s=lin.lstsq(A,b)
print(lin.norm(b-A.dot(x)))
Result
Ex 10.2
Ex 10.1
Code
import numpy as np
import matplotlib.pyplot as plt
import scipy.linalg as lin
m=3
n=3
A=np.random.random((m,n))
b=np.random.random((m,1))
x,res,rnk,s=lin.lstsq(A,b)
print(lin.norm(b-A.dot(x)))
Result
Ex 10.2