  In the framework of BCS model, we have applied the isothermal probability
distribution to take into account the statistical fluctuations in calculation
of thermodynamical properties of nuclei. The energy and the heat capacity are
calculated in $^{94}Mo$ nucleus using the mean gap parameter. The results are
compared with the values obtained based on the most probable values,
experimental data as well as some other theoretical models. We have shown that
heat capacity versus temperature behaves smoothly instead of singular behaviour
predicted by the standard BCS model. Also a smooth peak in heat capacity is
observed which is a signature of transition from normal to super fluid phase.
