In this paper, the authors propose a simple and numerically efficient interpolation algorithm capable to approximate the First Order Reversal Curves (FORC) distribution in any point. Using this algorithm as input for a classical Preisach Model one can obtain the magnetization curves of the system used in FORC measurements. In order to improve the precision of the calculus one can adapt different strategies form different variants of Preisach model.