The solution of systems of equilibrium is an important though commonplace operation, rendered difficult by non-linearity. Previous methods of numerical solution have suffered from slow convergence, unreliability, and inefficient structure. An algorithm, RAMESES, is presented in matrix algebra terms, which is both simple in structure and efficient, involving only one matrix inversion per system and giving an exact solution of the set of dependent equilibrium equations at each iteration.