function J = costFunction(X, y, theta)
m = size(X,1);
disp("No. of training examples");
disp(m);
predictions = X * theta;
disp("Predictions");
disp(predictions);
sqrErrors = (predictions-y).^2;
disp("Squared Errors");
disp(sqrErrors);
J = 1/(2*m) * sum(sqrErrors);