clc
clear
NC=2;
PROPS=zeros(NC,4);%N2-40% CH4-60%
PROPS(:,1)=[126.2;190.6];%critical temperature [K]
PROPS(:,2)=[34;45.99];%critical pressure [bar]
BIC=zeros(NC,NC);%Binary interaction coefficients
COMP=zeros(NC,1);%
COMP=[0.4;0.6];%mole fraction
ABCD=[NC,4];%Zarayeb

PRESS=30;%bar
TEMP=200;%K
phasetype = 0;
[ q,BETA,z,phi ] = EOS_CAL(NC,PRESS,TEMP,phasetype,COMP,PROPS,BIC);