Discrepancies in BCS theory results for finite Fermi systems with exact values of thermodynamic quantities are eliminated by means of averaging over pairing fluctuations. A form of the averaging function is proposed and its parameters are considered.