for a = 1:ml-1
    for c = 1:a
        k = m(ml-a,mc-c);
        for b = 1:mc
            m(ml-a,b)=m(ml-a,b)-k*m(ml+1-c,b);
        endfor
    endfor
endfor