Least square regression

function p = least_square_regression(x, y, degree)
    rows = length(x);
    A = zeros(rows, degree + 1);

    for i = 1:rows
        for j = 1: degree + 1
            A(i, j) = x(i).^(j-1);
        end
    end

    % Note: A'*A ~= A*A'
    p = (A'*A)\(A'*y');
    p = flipr(p');

end