Numerical determination of pseudorotation constants