Se considera un conjunto de sistemas singularmente perturbados que comprenden ecuaciones integro-diferenciales de Fredholm asociadas con problemas de reacción-difusión. Para aproximar soluciones a estos sistemas, se utiliza un esquema de segundo orden para las derivadas y la regla del trapecio para los términos integrales. La discretización se realiza en mallas no estándar conocidas como mallas de tipo Shishkin. El método numérico demuestra una tasa de convergencia de segundo orden con respecto a los parámetros pequeños en las ecuaciones, y se derivan estimaciones de error en la norma máxima discreta. Se realizan experimentos numéricos para verificar los resultados teóricos.