Multiconfiguration approximation nowadays is one of the main and most developed methods of investigation of atoms and ions. One of the basic stages in establishing this method was the work by A. Jucys [1], where the multiconfiguration Hartree–Fock equations were correctly and consistently derived for the first time. Taking into account this consideration and the following tradition of many years, in this paper these equations will be referred to as the Hartree–Fock– Jucys (HFJ) equations. The publication of the paper [1] is closely related to the fact that A. Jucys, still being a young scientist before the World War II, has been trained by R.D. Hartree, and in the beginning of fifties – by V. Fock. Right after the above-mentioned publication, a whole series of works were performed where these equations for particular atoms were solved both numerically and analytically [2–7]. This enabled determination of the main properties of the solutions of HFJ equations and formulation of simpler to solve and easier to perceive simplified (two-configuration) equations. In all the above-mentioned publications, as well as in many others, the equations were solved “by hand”, i. e. using only mechanical calculating machines. By the initiative of Prof. A. Jucys, the first electronic computer in Lithuania (BESM-2M) was obtained and a Computation Centre was established. This allowed development of a program for solving the ordinary Hartree–Fock (HF) equations numerically [8], i. e. considerable automation of calculations of the radial orbitals. Later on, such programs were developed in series for more advanced computers as well, first in the machine codes [9] and then in the FORTRAN algorithmic language [10]. The development of one’s own programs was justified mostly by the significant differences between the locally available computers and those from the West both in the structure and in the peculiarities of the algorithmic languages. The programs for numerical solution of the HF equations, developed at the Department of Theory of Atom (previously the Department of Quantum-Mechanical Computations), and many other programs, once were widely accepted in the scientific centres of the former Soviet Union. This was enhanced by the publication of 20 issues of “Collection of Programs for the Automation of Atomic Calculations” in our Institute from 1977 to 1988, where 48 FORTRAN programs were published. It must be noted that it was a personal initiative of Prof. A. Jucys to use FORTRAN for program development as early as in the beginning of the seventies of the previous century. This fact is now undoubted. It is the language that the overwhelming majority of the programs used for the atomic calculations all over the world are written in. The author of this paper developed the program for numerical solution of the simplified HFJ equations in the course of his thesis. It was communicated at the