The aim of this work considers a second order point reactor kinetics model based on the P1 approximation of transport theory, called in this work as P1 point reactor model. The P1 point reactor model implicitly considers the time derivative of the neutron source which has not been thus considered previously. The inverse method to calculate the reactivity in nuclear reactors -chosen because its high accuracy- Matrix Formulation. The numerical results shown that the Matrix Formulation for the reactivity estimation constitutes a method with insignificant calculation errors.