Let $A$ be a real symmetric invertible matrix and $b$ a real non-zero vector. Consider the problem of finding a real number non-zero $\lambda$ and a real valued vector $x$ such that $$Ax=\lambda x + b.$$
How can I numerically and efficiently solve this problem?