A program to solve rotating plasma problems

It is shown that the solution of a rotating plasma problem minimizes a quitably chosen funtional. This variational problem is solved by the Ritz-Galerkin methud using piecewise bilinear functions and applying some Newton-Côtes-like quadrature. The resulting linear system with a sparse nonegative definite matrix is dolved by an adapted conjugate gradient method.