The membrane bending stiffness of nearly spherical lipid vesicles can be deduced from the analysis of their thermal shape fluctuations. The theoretical basis of this analysis [Milner and Safran, Phys. Rev. A: At., Mol., Opt. Phys., 1987, 36, 4371-4379] uses the mean field approximation. In this work we apply Monte Carlo simulations and estimate the error in the determination of the bending stiffness due to the approximations applied in the theory. It is less than 10%. The method presented in this work can also be used to determine the changes of the bending stiffness of biological membranes due to their chemical and/or structural modifications.