El documento considera el modelo de un centro de llamadas en forma de un sistema de colas de varios servidores con llegadas de Poisson y un área de espera ilimitada. En el modelo en consideración, las llamadas entrantes no difieren en términos de condiciones de servicio, servicio solicitado y períodos entre llegadas. Se asume que una llamada entrante puede utilizar cualquier servidor libre y todos son idénticos en cuanto a capacidades y calidad. El problema objetivo es encontrar la distribución estacionaria del número de llamadas en el sistema para un servicio recurrente arbitrario. Esto nos permitirá evaluar las medidas de rendimiento de tales sistemas y resolver varios problemas de optimización para ellos. Considerar modelos con tiempos de servicio no exponenciales proporciona soluciones para una amplia clase de modelos matemáticos, lo que hace que los resultados sean más adecuados para centros de llamadas reales. La solución se basa en la aproximación de la función de distribución dada del tiempo de servicio por la función de distribución hiperexponencial. Por lo tanto, primero, se resuelve el problema de estudiar un sistema con servicio hiperexponencial utilizando el método matricial-geométrico. Posteriormente, sobre la base de este resultado, se construye una aproximación de la distribución estacionaria del número de llamadas en un sistema de varios servidores con una función de distribución arbitraria del tiempo de servicio. Se consideran varios problemas en la aplicación de esta aproximación, y se analiza su precisión en comparación con el resultado analítico conocido para un caso particular, así como con los resultados de la simulación.