As of 1 April 2026, the PSJD database will become an archive and will no longer accept new data. Current publications from Polish scientific journals are available through the Library of Science: https://bibliotekanauki.pl
The numerical stability of the solution of recursion relation for mean occupation numbers derived by Schönhammer for ideal Fermi gas trapped in 1D harmonic potential is studied. In low temperature region there exists a solution of this recursion relation. In high temperature region the iteration becomes unstable. In low and high temperature regions with growing number of particles the region of numerical instability diminishes.