El objetivo de este artículo es demostrar la viabilidad de utilizar modelos termodinámicos de coeficientes de actividad (Wilson, NRTL, UNIQUAC) programados con matrices, en lugar de estar programados con ciclos. Se determina la relación de equilibrio líquido-vapor de las mezclas Metanol-Etanol-Benceno y Acetona-Cloroformo-Metanol representados en mapas de curvas de residuo y en mapas de líneas de destilación. Para obtener resultados más confiables y conclusiones objetivas, el estudio fue apoyado con el uso del diseño estadístico de experimentos. La programación con matrices mostró ser mucho más rápida que el modelo cíclico convencional, independiente del tipo de modelo de coeficiente de actividad usado.