Este documento propone algoritmos para modelar comportamientos fraccionarios (dinámicos) utilizando núcleos racionales no singulares cuyo interés se demuestra primero en una función de ley de potencias pura. Luego se proponen dos algoritmos para encontrar un núcleo racional no singular que permita ajustar los datos de entrada-salida. El primero deriva la respuesta al impulso del sistema modelado a partir de los datos. El segundo encuentra los polos y ceros entrelazados de la función racional que ajusta la respuesta al impulso encontrada utilizando el primer algoritmo. Varias aplicaciones muestran la eficiencia del trabajo propuesto.