A two-dimensional non-linear equilibrium code, able to simulate reversed field pinch (RFP) discharges, has been developed. It is based on the two-dimensional MHD equilibrium code MAXFEA, which has been properly modified by introducing RFP profiles of the poloidal current density flux function in accordance with the μ&p model. In particular, equation parameters consistent with available data from the RFX machine have been used. A complete mesh of the RFX poloidal field system has been prepared and the whole poloidal electrical circuit has been implemented into the code. The RFP equilibrium code has been successfully validated against some RFX experiments, with and without plasma.