matrixInverseVector

function A_inv_b = matrixInverseVector(A, b, x_init, alpha)
    while norm(A*x_init - b)^2 > 1e-6,
        x_init -= alpha * 2*A*(A*x_init - b);
    end;
    A_inv_b = x_init;
end;