function [theta, J_history] = gradientDescent(X, y, theta, alpha, num_iters)
m = length(y);
J_history = zeros(num_iters, 1);
for iter = 1:num_iters
theta=theta-alpha*sum((X*theta-y).*X)'/m; %注意!对θ的更新要同步进行
J_history(iter) = computeCost(X, y, theta);
end
end