R<x> := PolynomialRing(RationalField());
X037 := HyperellipticCurve(x^6+14*x^5+35*x^4+48*x^3+35*x^2+14*x+1);
X037;
J037 := Jacobian(X037);
J037;
P := X037![0,1];
Q := X037![0,-1];
x := J037![P,Q];
x;
2*x;
3*x;
Order(x);   // point of infinite order.

