i = 1;
j = 1;

m = input('The X size of the matrix?');
n = input('The Y size of the matrix?');

while (i <= m && j <= n)
    maxi = i;
    for k = i+1:1:m
        if abs(JARED(k,j)) > abs(JARED(maxi,j))
            maxi = k;
        end
    end
    
    if JARED(maxi, j) ~= 0
        for u = i+1:1:m
            %SUBTRACTION HERE
        end
        i = i + 1;
    end
    j = i + 1;
end