def gd(X,y):
    theta = np.zeros(2)
    lr = 0.01
    theta1 = theta
    theta = gd_update_epoch(X, y, theta, lr)
    while abs_diff(theta, theta1) > 0.01:
        theta_hold = theta
        theta = gd_update_epoch(X, y, theta, lr)
        theta1 = theta_hold

    return theta