diff --git a/sympy/polys/modulargcd.py b/sympy/polys/modulargcd.py
index 6fa835b98a1b..9294bdba9a03 100644
--- a/sympy/polys/modulargcd.py
+++ b/sympy/polys/modulargcd.py
@@ -1009,7 +1009,7 @@ def _modgcd_multivariate_p(f, g, p, degbound, contbound):
     d = 0
     evalpoints = []
     heval = []
-    points = set(range(p))
+    points = list(range(p))
 
     while points:
         a = random.sample(points, 1)[0]
@@ -1605,7 +1605,7 @@ def _func_field_modgcd_p(f, g, minpoly, p):
     evalpoints = []
     heval = []
     LMlist = []
-    points = set(range(p))
+    points = list(range(p))
 
     while points:
         a = random.sample(points, 1)[0]