x_ff3 := FiniteField(3); x_ff3x := PolynomialAlgebra(x_ff3); x_ff3xy := PolynomialAlgebra(x_ff3x); x_f := x_ff3xy.1^2+1; x_F := FunctionField(x_f); x_O := MaximalOrderFinite(x_F); x_I := Element(x_O,[3,1])*x_O; Modinv(PrimitiveElement(x_O), x_I);