KKKR Ex. 3.16 polynomial: 38 terms, degree 6
Disc_A(phi):              38 terms, degree 6
Disc_A(phi) == KKKR polynomial : True
Disc_A(phi) == -KKKR polynomial: False
  Disc at p=q=r=373/2700: 51895117/72900000000
  Disc at p=q=r=4/27: 0
  Disc at p=q=r=427/2700: -77854483/72900000000
  n=2: exact Disc(lam X)/(lam^(n(n-1)) Vandermonde^2), lam=1e3,1e5,1e7: ['0.940036', '0.999400', '0.999994']
  n=3: exact Disc(lam X)/(lam^(n(n-1)) Vandermonde^2), lam=1e3,1e5,1e7: ['0.874761', '0.998918', '0.999989']
  n=4: exact Disc(lam X)/(lam^(n(n-1)) Vandermonde^2), lam=1e3,1e5,1e7: ['-2.388753', '0.979203', '0.999793']
  n=5: exact Disc(lam X)/(lam^(n(n-1)) Vandermonde^2), lam=1e3,1e5,1e7: ['421.406547', '0.669950', '0.997495']

KKKR boundary cross-check PASSED
