On the calculation of the poles of multivariate meromorphic functions using the symbolic-numeric two-point qd-algorithm (Q780406)

The Rutishauser qd (quotient-difference) algorithm is first recalled in the univariate case and it is explained how it can be applied also in two-point (or Laurent expansion) case. The bivariate extension follows by setting \((x,y)=(\lambda_1 z,\lambda_2 z)\) so that the problem becomes again univariate in \(z\) but with coefficients depending on a parametric vector \(\vec{\lambda}=(\lambda_1,\lambda_2)\). After convergence one obtains parametric curves of pole locations. This property is illustrated with several numerical examples.