for ( j = 0 ; j < r ; ++j ) {

                Aj = A[j];

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