divisors.lib
LIB "divisors.lib"; ring r=31991,(x,y,z),dp; ideal I = y^2*z - x*(x-z)*(x+3*z); qring Q = std(I); divisor A = makeDivisor(ideal(x,z),ideal(1)); divisor B = makeDivisor(ideal(x,y),ideal(1)); intmat M[4][4]= 1,4,0,0, 1,0,3,0, 0,0,0,2, 1,1,1,1; polytope PP = polytopeViaPoints(M); makePDivisor(list(list(PP,A),list(PP,B)));