An algorithm for solving the pulsar equation

We present an algorithm of finding numerical solutions of pulsar equation. The problem of finding the solutions was reduced to finding expansion coefficients of the source term of the equation in a base of orthogo- nal functions defined on the unit interval by minimizing a multi-variable mismatch function defined on the light cylinder. We applied the algorithm to Scharlemann & Wagoner boundary conditions by which a smooth solu- tion is reconstructed that by construction passes success- fully the Gruzinov's test of the source function exponent.