"Source: Text/Geometry/CrvHyp.text";
"Line: 4910";
"Date: Thu Aug 18 16:02:10 2022";
"Main: Thu May 25 10:02:38 2023";
// original file: Text/Geometry/CrvHyp.text, line: 4910
// Example: H134E35 ()
print "Example: H134E35";
ei := GetEchoInput();
SetEchoInput(true);
_<x> := PolynomialRing(Rationals());
C := HyperellipticCurve(x^6+x^2+2);
ptsC := Points(C : Bound:=1000); ptsC;
J := Jacobian(C);
RankBound(J);
assert $1 eq 1;
PJ1 := J! [ ptsC[3], ptsC[1] ]; 
Order(PJ1);
assert $1 eq 8;
PJ2 := J! [ ptsC[5], ptsC[1] ]; 
Order(PJ2);
assert $1 eq 0;
all_pts := Chabauty(PJ2); all_pts;
okay := all_pts eq Chabauty(12*PJ2);
assert okay; // Check that it is okay
BadPrimes(C);
Chabauty(PJ2, 3);
SetEchoInput(ei);
