In this paper, we have presented the optimized BCS formalism using the isothermal probability distribution. The effect of statistical fluctuations on thermodynamical properties of nuclei has been investigated. The average gap parameter is calculated and then the energy, the entropy and the heat capacity are evaluated. The resulting values are compared with results of the standard BCS and the Static Path Approximation plus Random Phase Approximation, (SPA+RPA). We have shown that the resulting heat capacity versus temperature using the optimized BCS model is very similar in shape to the experiments in the other models. The peak in the heat capacity is interpreted as the transition from the paired to the normal phase.