function m=reg(y,x)

m=inv(x'*x)*(x'*y);