The projections based on harmonic polynomials offer many advantages over traditional ones, among them the minimization of the scale alteration while allowing to represent a territory with topological continuity. Algorithms are usually implemented as “closed” functions without passing the polynomial coefficients as parameters, which forces to create a new function in case of small changes, for example when using different ellipsoids. This discourages the popularization of these advantageous projections. However, it is very simple to create algorithms that receive the coefficients as input parameters. In addition, the ellipsoid change does not significantly affect the scale distribution of the projection, so these coefficients are valid for any terrestrial ellipsoid currently used. On the other hand, using only two types of rectangular variables, Mercator and Gauss-Schreiber, any region of the world can be represented, be it polar or equatorial. This paper proposes a way to implement the algorithms to be used in any geographic software; also provides an algorithm for the reversion of polynomials with micro-metric accuracy.