function beta = leastSquaresGD(y, tX, alpha)
%Computes beta using gradient descent

  % initialize
  beta = zeros([size(tX, 2) 1]);
  maxIters = 2000;

  % iterate
  for k = 1:maxIters
    % INSERT YOUR FUNCTION FOR COMPUTING GRADIENT
    g = computeGradient(y, tX, beta);

    % INSERT YOUR FUNCTION FOR COMPUTING COST FUNCTION
    %L = computeCost(y, tX, beta);

    % INSERT GRADIENT DESCENT UPDATE TO FIND BETA
    beta = beta - alpha.*g;

    % INSERT CODE FOR CONVERGENCE
    if (norm(g) <  0.00000001)
       break
    end
  end
end