load(input('Ange datafil: ','s'))
tic;
[m,n] = size(A);

% Skapa nix, dvs indexvektorn för ickebasvariabler
nix = setdiff([1:n],bix);

% Skapa initial partition
B = A(:,bix);
N = A(:,nix);
cB = c(bix,:);
cN = c(nix,:);

opt=0;
iter = 0;
while opt==0,
    iter = iter + 1;

    % Beräkna reducerad kostnad och högerled
    % -------

    
    
    % Beräkna mest negativ reducerad kostnad, rc_min,
    %  och index för inkommande variabel, inkix
    % -------

    

    if rc_min >= -1.0E-10
        opt=1;
        disp('Optimum');
    else
        % Beräkna inkommande kolumn, a
        % --------

        
        if max(a) <= 0 
            disp('Obegränsad lösning');
        else
            % Bestäm utgående variabel, utgix
            % -------



            disp(sprintf('Iter: %d, z: %f, rc_min: %f, ink: %d, utg: %d ',iter,z,rc_min,nix(inkix),bix(utgix)));

            % Konstruera ny partitionering mha ink och utg
            % --------
            
            
            
        end
    end
end

toc
disp(sprintf('z: %f',z));
x = zeros(n,1);
x(bix) = xb;
disp(sprintf('sum(x-xcheat): %f',sum(x-xcheat)));
disp(sprintf('z-zcheat: %f',z-zcheat));

