Sub RED(IZ1, IZ2, JZ1, JZ2, JM1, JM2, JMF, IC1, JC1, JCF, KC, C(), NCI, NCJ, NCK, S(), NSI, NSJ)
    LOFF = JC1 - JM1
    IC = IC1
    For J = JZ1 To JZ2
        For L = JM1 To JM2
            VX = C(IC, L + LOFF, KC)
            For I = IZ1 To IZ2
                S(I, L) = S(I, L) - S(I, J) * VX
            Next I
        Next L
        VX = C(IC, JCF, KC)
        For I = IZ1 To IZ2
            S(I, JMF) = S(I, JMF) - S(I, J) * VX
        Next I
        IC = IC + 1
    Next J
End Sub