  A numerical method optimizing the coefficients of the semi empirical mass
formula or those of similar mass formulas is presented. The optimization is
based on the least-squares adjustments method and leads to the resolution of a
linear system which is solved by iterations according to the Gauss-Seidel
scheme. The steps of the algorithm are given in detail. In practice the method
is very simple to implement and is able to treat large data in a very fast way.
In fact, although this method has been illustrated here by specific examples,
it can be applied without difficulty to any experimental or statistical data of
the same type, i.e. those leading to linear system characterized by symmetric
and positive-definite matrices.
