An algorithm based on the choice of a single monotone function (subject to boundary conditions) is presented which generates all regular static spherically symmetric perfect fluid solutions of Einstein's equations. For physically relevant solutions the generating functions must be restricted by non-trivial integral-differential inequalities. Nonetheless, the algorithm is demonstrated here by the construction of an infinite number of previously unknown physically interesting exact solutions.

Final form to appear in Phys Rev D. Includes a number of clarifications