The paper concerns numerical analysis of superparamagnetic clusters by means of Langevin paramagnetism theory and simulated annealing (SA) procedure supplemented with some additional constraints. In order to determine a proper calculation algorithm and its parameters the three methods, i.e., simple SA, SA with maximum entropy, and SA with local maximum entropy were tested. Within the analyzed SA procedures (carried out for the computer generated semi-empirical curves) the best results were obtained by making use of the local entropy constraint. It was also shown that this algorithm can be successfully applied to determination of superparamagnetic cluster distribution. The tested SA procedure was applied to analysis of measured magnetic isotherms for Al-RE-Ni amorphous alloys and thin SiC-Mn layers deposited on a Si substrate.