for ( j = r + 1 ; j < m ; ++j ) {

                Aj = A[j];

                f = div( Aj[c], pivot );