% Este algoritmo resuelve problemas del tipo:
%
%==================================================================
%             Min d'X+X'H X    H (matríz positiva definida de mxn).
%             s.a. AX<=b
%                   X>=0
%==================================================================

% Aproximación inicial Xo;
     X0=(-.5)*((inv(H))*d;
     
% Los multiplicadores de lagrange estan dados por:
    lamda=inv(.5*Aa*inv(H)*(Aa)')*(-ba-.5*(Aa')*(inv(H))*d); 

% La distancia perpendicular del punto Xo a cada una de las 
% restricciones está dada por D:
    D=(aj*(X_)-bj)/(sqrt(aj*(inv(H))*aj')

    % De todas las distancias del punto inicial a cada una de las 
    % restricciones tomamos la máxima. Y de aquí obtenemos nuestro
    % siguiente aproximación X1.
    
    % Calculo el valor de lambda correspondiente a esta restricción
    % si lambda < 0 entonces esta restricción se tiene que quitar
    % si lambda > 0 entonces vuelvo a aplicar el método pero con la 
    % restricción fija.
    
    
    % Continúo así hasta obtener la solución.
    
    














