цикл при помощи индексов

Здесь написан обратный ход для слау 4 на 5 Помогите пожалуйста написать его в общем случае, т.е при помощи индексов P.S это пример для обычной матрицы , но для слау сделать никак не получается

for (int j = size-1; j > 0; j--)//обратный ход
    {
        int iDop = j - 1;
        for (int i = j-1; i >= 0; i--, iDop--)
        {
            //R[i][j] = R[i][j] - R[j][j] * (double(A[iDop][j] / A[j][j]));
            A[i][j] = A[i][j] - A[j][j] * (double(A[iDop][j] / A[j][j]));
        }
    }

tmp=a[2][3]/a[3][3];
a[2][3]=a[2][3]-a[3][3]*tmp;
a[2][3]=a[2][4]-a[3][4]*tmp;

tmp=a[1][3]/a[3][3];
a[1][3]=a[1][3]-a[3][3]*tmp;
a[1][4]=a[1][4]=a[1][4]-a[3][4]*tmp;

tmp=a[0][3]/a[3][3];
a[0][3]=a[0][3]-a[3][3]*tmp;
a[0][4]=a[0][4]-a[3][4]*tmp;
___________________________
tmp=a[1][2]/a[2][2];
a[1][2]]=a[1][2]-a[2][2]*tmp;
a[1][3]=a[1][3]-a[2][3]*tmp;
a[1][4]=a[1][4]-a[2][4]*tmp;

tmp=a[0][2]/[2][2];
a[0][2]=a[0][2]-a[2][2]*tmp;
a[0][3]=a[0][3]-a[2][3]*tmp;
a[0][4]=a[0][4]-a[2][4]*tmp;
_____________________________
tmp=a[0][1]/a[1][1];
a[0][1]=[0][1]-a[1][1]*tmp;
a[0][2]=a[0][2]-a[1][2]*tmp;
a[0][3]=a[0][3]-a[1][3]*tmp;
a[0][4]=a[0][4]-a[1][4]*tmp;