bivariateDiagEntries -- computes diagonal entries and eigenvalues for a determinantal representation of a bivariate polynomial

This method computes the eigenvalues and diagonal entries of a monic symmetric determinantal representation of a real bivariate polynomial $f$, or gives an error if certain necessary conditions for existence of such a representation are not met. For a symmetric determinantal representation $f = det(I + x_1A_1 + x_2A_2)$, this method computes diagonal entries and eigenvalues of $A_1$ and $A_2$. The output is a 4-tuple of column vectors: (eigenvalues of A_1, eigenvalues of $A_2$, diagonal entries of $A_1$, diagonal entries of $A_2$).

When working over an InexactFieldFamily like RR or CC, the option Tolerance can be used to specify the internal threshold for checking equality (any floating point number below the tolerance is treated as numerically zero).