You need to provide an irreducible bivariate polynomial
R.<x,y> = ZZ[]
f1 = 127*x*y - 1207*x - 1461*y + 21
as well as bounds for the solutions and a parameter k
X, Y = 30, 20
k = 2
Then we can call the solver
bivariate_solver(f1, X, Y, k)
which returns the solutions
[(21, 21), (23, 19)]
This project is licensed under the MIT license.