%Diagrama P-x-y

X1=0:0.05:1; %fracción molar del compuesto 1 en el líquido (el más volátil)

X2=1-X1 %fracción molar del compuesto 2 en el líquido

A1=input('Introduzca el valor de la constante A de la sustancia 1')

B1=input('Introduzca el valor de la constante B de la sustancia 1')

C1=input('Introduzca el valor de la constante C de la sustancia 1')

A2=input('Introduzca el valor de la constante A de la sustancia 2')

B2=input('Introduzca el valor de la constante B de la sustancia 2')

C2=input('Introduzca el valor de la constante C de la sustancia 2')

T=input('Introduzca la temperatura en °C')

%Determinación de la presión de vapor del compuesto 1 puro

Po1=10.^(A1-(B1/(T+C1)));

%Determinación de la presión de vapor del compuesto 2 puro

Po2=10.^(A2-(B2/(T+C2)));

%Determinación de la presión de los componentes en el sistema

P1=Po1*X1

P2=Po2*X2

%Determinación de la presión total de la mezcla

P=P1+P2

%Determinación de la fracción molar del compuesto 1 en el vapor

Y1=P1./P

plot(X1,P,Y1,P)