function cost = J(theta,X,y)
  yhat = X*theta';
  cost = (y-yhat)'*(y-yhat);


