models_X0Nplus_g2g3.m

Qx<x> := PolynomialRing(RationalField());
Qxy<y> := PolynomialRing(Qx);

f67  := x^6 + 2*x^5 +   x^4 - 2*x^3 + 2*x^2 - 4*x + 1; 
f73  := x^6 + 2*x^5 +   x^4 + 6*x^3 + 2*x^2 - 4*x + 1; 
f103 := x^6 + 6*x^5 + 5*x^4 + 2*x^3 + 2*x^2       + 1;
f107 := x^6 + 2*x^5 + 5*x^4 + 2*x^3 - 2*x^2 - 4*x - 3;
f167 := x^6 - 4*x^5 + 2*x^4 - 2*x^3 - 3*x^2 + 2*x - 3;
f191 := x^6         + 2*x^4 + 2*x^3 + 5*x^2 - 6*x + 1;
Q_97  := y^3 + (-x^2+x)*y^2 +(-x^3-x^2-1)*y + (x^3 + x^2); 
Q_109 := x*y^3 + y^4 - x^3 - 4*x^2*y - 4*x*y^2 - x^2 - 2*x*y + 2*y^2 + x + 3*y; 
Q_113 := x^3*y - 3*x^2*y^2 - 4*x^2*y - 2*x^2 + x*y^3 + 3*x*y^2 + 3*x*y + x + y^4 + 2*y^3 + 2*y^2 + y; 
Q_139 := x^3 - x^2*y^2 + x^2*y - x*y^3 - 2*x*y^2 - 3*x*y - x + y^4 + y^3 + y^2 + y; 
Q_149 := -(x^3 - x^2*y^2 + x*y^3 + x*y^2 - 2*x*y - x - y^4 + y^3 + y^2 - y); 
Q_151 := 20*x^3*y + 24*x^3 + 22*x^2*y^2 + 20*x^2*y - 40*x^2 + 8*x*y^3 + 6*x*y^2 - 32*x*y + 16*x + y^4 + y^3 - 6*y^2 + 4*y; 
Q_179 := 2*x^3*y + x^3 + x^2*y^2 + x^2*y + 3*x*y^3 + 3*x*y^2 - x*y - x + y^4 + 4*y^3 + 4*y^2 + y;  
Q_239 := x^3 + (-y^2 + y + 1)*x^2 + (-y^3 - y^2 - y)*x + (y^4 + 3*y^3 + 2*y^2 + y); 
Q_127 := 20*x^3*y + 24*x^3 + 4*x^2*y^2 - 20*x^2*y - 32*x^2 + 9*x*y^3 + 18*x*y^2 + 20*x*y + 16*x + y^4 + 4*y^3 + 6*y^2 + 4*y; 


x0plus_quartics_list := [Q_97, Q_109, Q_113, Q_127, Q_139, Q_149, Q_151, Q_179, Q_239 ];
x0plus_hyp_list := [f67, f73, f103, f107, f167, f191];
levels_quartics_list := [97, 109, 113, 127, 139, 149, 151, 179, 239 ];