En este artículo consideramos el problema de valores propios complementarios, de gran interés para muchos investigadores por sus numerosas aplicaciones en Ingeniería y Física, y abordamos su solución como un problema de complementariedad no lineal usando un método cuasi-Newton, un tipo de método que, hasta donde conocemos, no ha sido utilizado para dicho fin. Verificamos que el problema satisface ciertas hipótesis que permiten utilizar un algoritmo global cuasi-Newton y realizamos pruebas numéricas que muestran la eficacia del algoritmo utilizado y lo hacen una buena alternativa de solución para problemas de valores propios complementarios.