Os autores apresentam uma abordagem semi-analítica para obter soluções para as equações da cinética pontual de nêutrons. A reatividade dependente do tempo é aproximada por aproximação constante por partes em tempo discretizado, e o sistema de equações diferenciais é resolvido por um fator de integração. O termo de fonte externa na solução possui uma integral, que é resolvida numericamente pela interpolação de polinômios de Lagrange e uma fórmula simples de quadratura. Mesmo que a fonte externa apresente não linearidade nas equações, a metodologia é apropriada. A metodologia proposta é utilizada em casos de referência padrão e seus resultados são comparados com os da literatura.