This algorithm solves the otimization problem min_x \| y - A x \|_2^2 + lambda Phi(x) and has applications in compressed sensing, image restoration and reconstruction, sparse regression, and several other problems.