The integro-differential equation (IDE) of the reactor point kinetics was numerically solved as an ODE. The integral term of the IDE was taken as an independent term of the ODE and it was calculated with an estimate of the neutron density (ND) using the Lagrange interpolation polynomials, then the exact solution of the resultant ODE (IQLn algorithms) is obtained. Four applications were used to assess the accuracy of the algorithms: a mild ramp, a fast ramp, a sinusoidal and a feedback type of reactivity insertion. The 2nd degree polynomial algorithm (IQL 2 ) showed that for about 3-to-4-digit accuracy and Δt⩽10-3, it was as accurate as other highly accurate codes, and that it is very efficient.