// Problema de blending e producao de uma refinaria
// Edgar & Himmelblau, 1988
// 
//  max c'x = min -c'x
//      A x <= b
//        x >= 0

c=[8.1; 10.8];
A=[0.8   0.44
    0.05  0.1
    0.1   0.36];
b=[24000; 2000; 6000];
lb=zeros(2,1);
ub=[];
x0=[0; 0];

// active set: simplex
//exec lipsol.sci;
//[x,S,eflag,out,lambda]=lipsol(-c,A,b,[],[],lb,ub)

// interior point: S. Mehrotra, 1992
[xi,Si,eflagi,iter,outi]=karmarkar([],[],-c,[],[],[],[],[],A,b,lb,ub)
