num_iters = 1500
alpha = 0.01
theta = c(0,0)
X = pop$V1
y = pop$V2
X = cbind(1,X)
# Call GradientDescent()
result <- GradientDescent(X,y,theta,alpha,num_iters )
theta <- result$theta
j_hist <- result$J_history
summary(j_hist)