LIB "brnoeth.lib";
kill X, X2,R,G,LG;
ring R=11,(x,y),lp;
list X = Adj_div(x^7 + y^7 - 1);
def X2 = NSplaces(1,X);
def X3 = extcurve(1,X2);
def RR =X3[1][5];
setring RR;
print("POINTS");
print(POINTS);
/* PROBLEM -- this G defined a different divisor every time the
this code is run!!! Need a way to compute G from a list of points */
intvec G=(10,-1,0,0,9,0,0,0,0,0,0,0,0,0);
def R = X2[1][2];
setring R;
list LG = BrillNoether(G,X2);
print(LG);
