El objetivo de este trabajo es obtener una fórmula de representación integral para las soluciones de problemas de valor inicial para sistemas neutrales fraccionarios lineales autónomos con derivadas de tipo Caputo y retardos distribuidos. Los resultados obtenidos mejoran y extienden los resultados correspondientes en el caso particular de sistemas fraccionarios con retardos constantes y serán una herramienta útil para estudiar diferentes tipos de propiedades de estabilidad. Los resultados propuestos coinciden con los correspondientes para sistemas diferenciales lineales neutros de primer orden con derivadas de orden entero.