import numpy as np 
x0=np.ones((10))
#print(x0)
x1=[64.3,99.6,145.45,63.75,135.46,92.85,86.97,144.76,59.3,116.03]
x2=[2,3,4,2,3,4,2,4,1,3]
y=np.array([62.55,82.42,132.62,73.31,131.05,86.57,85.49,127.44,55.25,104.84])
#y=[62.55,82.42,132.62,73.31,131.05,86.57,85.49,127.44,55.25,104.84]
X=np.stack((x0,x1,x2),axis=1)
X1=np.mat(X)
#print(type(X1))
#print(X)
Y=y.reshape(10,1)
Y1=np.mat(Y)
#print(Y)
a=(X1.T)*X1
#print(a)
W=(a.I)*(X1.T)*Y1
#print(W)
print("X=",X)
print("Y=",Y)
print("W=",W)
print("W的shape=",W.shape)