\\ Get and compile the branch "nicolas-crv", and cd to the subdirectory "Algcurves" of the pari directory.

read("Algcurves.gp");

B=Branches0(y^3+2*x^3*y-x^7,t,a)[2][1];
BranchExpand(B,10)
BranchExpand(B,100)

C=CrvInit(y^3+2*x^3*y-x^7);
CrvPrint(C);

C1=CrvInit(-256*x^56 + 6144*x^55 - 62464*x^54 + 333824*x^53 - 859648*x^52 - 120832*x^51 + 7252992*x^50 - 16046080*x^49 - 9891072*x^48 + 90136576*x^47 - 73076736*x^46 - 237805568*x^45 + 420485120*x^44 + 341843968*x^43 - 1165840384*x^42 - 192667648*x^41 + 2178936320*x^40 - 238563328*x^39 - 3063240704*x^38 + 639488000*x^37 + 3412593664*x^36 - 639488000*x^35 - 3063240704*x^34 + 238563328*x^33 + 2178936320*x^32 + 192667648*x^31 - 1165840384*x^30 - 341843968*x^29 + (-288*y^4 + 420485120)*x^28 + (3456*y^4 + 237805568)*x^27 + (-14400*y^4 - 73076736)*x^26 + (14976*y^4 - 90136576)*x^25 + (56160*y^4 - 9891072)*x^24 + (-142848*y^4 + 16046080)*x^23 + (-52992*y^4 + 7252992)*x^22 + (400896*y^4 + 120832)*x^21 + (-55872*y^4 - 859648)*x^20 + (-624384*y^4 - 333824)*x^19 + (134784*y^4 - 62464)*x^18 + (624384*y^4 - 6144)*x^17 + (-55872*y^4 - 256)*x^16 + (16*y^6 - 400896*y^4)*x^15 + (-96*y^6 - 52992*y^4)*x^14 + (-384*y^6 + 142848*y^4)*x^13 + (3232*y^6 + 56160*y^4)*x^12 + (-5424*y^6 - 14976*y^4)*x^11 + (960*y^6 - 14400*y^4)*x^10 - 3456*y^4*x^9 + (960*y^6 - 288*y^4)*x^8 + 5424*y^6*x^7 + 3232*y^6*x^6 + 384*y^6*x^5 - 96*y^6*x^4 - 16*y^6*x^3 + 27*y^8);

C[7] \\ yx/(2x^3+3y^2) dx, x^3/(2x^3+3y^2) dx
u = C[7][1][1]/C[7][1][2]
w = x
factor(MorImg(y^3+2*x^3*y-x^7,u,w))
poldisc(%[2,1],y)
DivPrint(FnDiv(C,u-2/3))

CrvPrint(C)
RiemannRoch(C,[2,5])
L=RiemannRoch(C,[[-1,1],3;3,1;1,-2])
DivPrint(FnDiv(C,L[1]))

C1 = CrvInit((x+y+1/x+1/y+1)*x*y);
CrvPrint(C1)
CrvEll(C1,[1,0,0])

f = x^5+y^4+x^2*y^3;
C0 = CrvInit(f);
CrvPrint(C0)
[X,Y] = CrvRat(C0,1)
substvec(f,[x,y],[X,Y])

C=CrvInit(x^5 + y^5 - 6*x^3 + 6*x^2 + x*y - 3*y^2);
CrvPrint(C)
CrvPicTorsGalRep(C,2,13,700)