%Time test for the different flash calculations.
filename='SSresults.mat';
load (filename)


%For PH Flash
%Definition of initial points
NC=5;
x0=[0.0080 0.0860 0.4590 0.3450];
x=x0;
y0=[0.1620 0.4570 0.3280 0.1840]; 
y=y0;
x(NC)=1-sum(x(1:NC-1));        
y(NC)=1-sum(y(1:NC-1));
VF0=0.65;
[Zl0,phil,Hl,Vl]= srks(z(16,:),T0,P(15),1);
[Zg0,phig,Hg,Vg]= srks(z(16,:),T0,P(15),2);
Hz= 

%Explicit compressibility
tic
[xeo, yeo, VFeo, Teo] = flashPHeo(  P(15), Hz,  x0, y0, z(16,:), T0, VF0, Zl0, Zg0 );
toc
Teo=toc

%Implicit compressibility
tic
[xs, ys, vfs, Tss] = flashPH( P(15), Hz, x0, y0, z(16,:), T0, VF0 );
toc
Ts=toc
