The Ising model for finite systems, e.g. for clusters of different sizes and crystal lattices, was solved analytically by the mean field approach. The magnetization was calculated from the number of accessible microstates, using the gamma function and its derivatives, unlike the usual solution in the microcanonical which uses the Stirling approximation. We determined a scaling exponent of ̃1/3, which shows how the Curie temperature decreases with decreasing nanoparticle size. Moreover, the model predicts the behaviour of surface and core regions and it explains in simple terms several effects previously observed in experiments and Monte Carlo simulations of small magnetic systems.