The exact analytical solution of inhomogeneous differential equations, resulting from the improvement of the known methods due to Post (1956) and Hall (1967), is presented. As an illustration, the energy lower bounds of N-fermion systems with oscillator and gravitational interaction are calculated and the results compared with the earlier data.