En este trabajo se propone el esquema de Runge-Kutta implícito de orden 1.5 para solucionar de forma numérica las ecuaciones de la cinética puntual estocástica, considerando diferentes valores de reactividad. Se obtiene la matriz de covarianza usando la teoría de difusión de una velocidad para un reactor nuclear. Para esto, se asume que en el tiempo inicial un pulso de neutrones es inyectado en un ensamble no multiplicativo. Los resultados obtenidos con este método logran una buena precisión en la desviación estándar y el valor esperado en comparación con otros métodos reportados en la literatura. Se concluye que el método propuesto es eficiente para el cálculo de la densidad de población de neutrones y la concentración de precursores de neutrones atrasados.