Least-squares fitting applied to nuclear mass formulas. Solution by the Gauss-Seidel method

In this paper, 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.