Division polynomials of Y^2 = X^3 + 1 computed up to n = 9.

n = 3: |P^1[3]| = degree = 2 (expected 3*(n/3)^2 - 1 = 2)
  factor (deg 2): u**2 - u + 1
  P^1[3] = {-eps, -eps^2}: True

n = 6: |P^1[6]| = degree = 11 (expected 3*(n/3)^2 - 1 = 11)
  factor (deg 1): u - 2
  factor (deg 1): u + 1
  factor (deg 1): 2*u - 1
  factor (deg 2): u**2 - u + 1
  factor (deg 6): u**6 - 3*u**5 + 60*u**4 - 115*u**3 + 60*u**2 - 3*u + 1
  linear factors: [u - 2, u + 1, 2*u - 1]
  the sextic (monic): u**6 - 3*u**5 + 60*u**4 - 115*u**3 + 60*u**2 - 3*u + 1
  equals (u^2-u+1)^3 + 54u^2(u-1)^2 : True
  BHRS's printed P(u) = u^6-3u^5+5u^3-3u+1 equals (u^2-u-1)^3+2 : True ; equals (u^2-u+1)^3 - 6u^2(u-1)^2 : True
  BHRS's P equals the computed sextic: False
  minimal polynomial of tau = -eps(b eps^2-1)/(b eps-1), b^3 = 2: x**6 - 3*x**5 + 60*x**4 - 115*x**3 + 60*x**2 - 3*x + 1
  equals the computed sextic: True
  J on the computed sextic's roots = -54; J on BHRS's P roots = 6 (identities above)
  J(tau) numerically: -54.0 + 1.92592994438723585305597794258e-34*I
  BHRS C_{v tau} proportional to -tau x^3 + (tau-1) y^3 + z^3 (ratio checks = 0): 0 0
  numeric P^1[6]: [-1.00000000000, 0.500000000000, 2.00000000000, 0.00877999689001 - 0.132222938383*I, 0.00877999689001 + 0.132222938383*I, 0.5 - 0.866025403784*I, 0.5 + 0.866025403784*I, 0.5 - 7.52978275729*I, 0.5 + 7.52978275729*I, 0.99122000311 - 0.132222938383*I, 0.99122000311 + 0.132222938383*I]

n = 9: |P^1[9]| = degree = 26 (expected 3*(n/3)^2 - 1 = 26)
  factor (deg 2): u**2 - u + 1
  factor (deg 3): u**3 - 6*u**2 + 3*u + 1
  factor (deg 3): u**3 + 3*u**2 - 6*u + 1
  factor of degree 18

Theorem 1.5 class: sum over 4 pencils of k(3L - E_t) = 12k L - k * (sum_t E_t) = 12k L - 3k E;
  12k - printed coefficient of L: 0 ; correct E-coefficient 3k = 9*m**2 - 3 ; printed 9k = 27*m**2 - 9
