function [sol,flag]=simplex(A,b,c)
%% La función [sol,flag]=simplex(A,b,c)
%% resuelve el problema
%%  min c'x,     S.A.  A.x <= b, x>=0
%%
%% CON b>=0 y indBase=indices de las vars. de los vectores base
%%              (vars. holgura o artificiales)
%%
%% y devuelve
%% flag= 0, si no existe solución,
%%     = 1, si la solución es única,
%%     = 2, si existe un segmento de soluciones,
%%     = inf, si la solución es no acotada
%% sol = solucion, en los casos (flag=1 o =2)
%%     = NaN, en otro caso (flag=0 o =inf).
%%

m=length(b);
n=length(c);

Aholgurada=[A eye(m)];
cholgurada=[ c ; zeros(m,1) ];
indBase=n+1:n+m;

[sol,flag]=simplex2fases(Aholgurada,b,cholgurada,indBase);
sol=sol(1:n);
