NUMERICAL CALCULATION OF THE COMPLEX BERRY PHASE IN NON-HERMITIAN SYSTEMS