p = 3
m = 2
F = GF(p)
Rx.<x> = PolynomialRing(F)
f = x^3 - x
C = superelliptic(f, m)
a = superelliptic_drw_form(C.one, 0*C.dx, 0*C.x)
b = a+a+a+a+a+a+a+a+a
print(b)