function B = least_squares(X,Y) 

B = Y*pinv(X);