Nuclear reactivity and neutron population density are critical parameters for quantifying and describing the behavior of a nuclear reactor. Accurate knowledge of reactivity is essential for the safe control and operation of a reactor. Various approaches and methods for calculating reactivity are reported in the literature. In this paper, we propose using the Euler-Maclaurin formula to numerically solve reactivity through the inverse point kinetics equation, a key model in the design of digital reactivity meters. By approximating an integral as a sum in the expansion of a continuous integral into a discrete version, and considering only the first three Bernoulli numbers while reducing higher-order derivatives, we derive a model to calculate nuclear reactivity as a function of neutron population density and physical constants in thermal reactors using uranium as nuclear fuel.