function res=itera_www(N_est)

%Fracciones molares en la superficie del l´iquido
x1_0=0.319;
x2_0=0.528;

x0=[x1_0, x2_0, 1-x1_0-x2_0];


% Longitud del tubo
zL=0.238; %m

z=[0 zL];

[z,x_cal]=ode15s(@ec_maxstef_2,z,x0,[],N_est);


res=x_cal(end,1:2)-[0,0]