[en] A code VMOMS is described which finds approximate solutions to the Grad-Shafranov equation describing scalar pressure-balance equilibria for axisymmetric tokamak plasmas. A Fourier series expansion of the flux surface coordinates (R,Z) is made in terms of two new coordinates (rho, theta), and the resulting equation is conveniently reduced to a system of ordinary differential equations (ODE's) using a variational principle. The solution of these simple equations with pressure and current as driving functions, yields, in principle, a complete description of the equilibrium. Complete axisymmetry is assumed, as well as up-down symmetry about the toroidal midplane