En este trabajo, consideramos un sistema de ecuaciones diferenciales ordinarias con difusión discreta no local y retardo finito y con un número finito o infinito de ecuaciones. Probamos varias propiedades de las soluciones como comparación, estabilidad y simetría. Creamos una simulación numérica que muestra que este modelo puede ser apropiado para modelar tablas de vida dinámicas en ciencias actuariales o demográficas. De esta manera, algunos indicadores de bondad y suavidad se mejoran al comparar con técnicas clásicas.