function deriv = dJ(theta,X,y)
  deriv = (y-X*theta')'*X;



