En el problema de autovalores de respuesta lineal que surge de la química cuántica y la física, es necesario calcular varios autovalores positivos junto con los correspondientes autovectores. Para tal tarea, en este artículo presentamos un algoritmo FEAST basado en integración de contornos complejos para el problema de autovalores de respuesta lineal. Al dividir simplemente el espectro en una colección de regiones disjuntas, el algoritmo es capaz de paralelizar el proceso de resolución del problema de autovalores de respuesta lineal. Se establecen los resultados de convergencia asociados para revelar la precisión del espacio propio aproximado. Se presentan ejemplos numéricos para demostrar la efectividad de nuestro algoritmo propuesto.