Se presenta el modelado de los equilibrios líquidos multifásicos de sistemas cuaternarios utilizando los modelos de coeficiente de actividad NRTL y UNIQUAC. El procedimiento numérico utilizado para calcular los equilibrios líquidos multifásicos de estos sistemas, está basado en la minimización de la energía de Gibbs del sistema con pruebas de estabilidad termodinámica para encontrar el estado más estable del sistema. Los algoritmos que se presentan para la estimación de los parámetros de interacción binaria de los modelos NRTL y UNIQUAC fueron aplicados en la correlación de datos experimentales para dos sistemas cuaternarios con dos y tres fases líquidas en equilibrio. La concordancia entre los datos experimentales y los calculados con ambos modelos termodinámicos fue satisfactoria para los dos sistemas estudiados.