w1=[0.4, 0.2, 0.2, 0.1, 0.1];
w2=[0.6, 0.2, 0.1];

C= [ 3, 5, 2;
     0, 2, 5;
     1, 1, 3;
     8, 4, 3;
     7, 6, 5 ];

[e,Flow]=emd_mex(w1,w2,C)