ydot=MaxwellStefan(t,y)

function ydot=MaxwellStefan(t,y)

global N1 N2 N3

N3=0;
F12 = 1.286e-3; 
F13 = 2.081e-3; 
F23 = 3.019e-3;

ydot(1)=(y(1)*N2-y(2)*N1)/F12+(y(1)*N3-(1-y(1)-y(2))*N1)/F13;
ydot(2)=(y(2)*N1-y(1)*N2)/F12+(y(2)*N3-(1-y(1)-y(2))*N2)/F23;

ydot=ydot';