diff -up src/sage/sandpiles/sandpile.py.orig src/sage/sandpiles/sandpile.py
--- src/sage/sandpiles/sandpile.py.orig	2020-10-24 11:37:59.000000000 -0600
+++ src/sage/sandpiles/sandpile.py	2020-10-30 14:00:00.744208307 -0600
@@ -347,9 +347,6 @@ from sage.arith.all import falling_facto
 from sage.rings.all import Integer, PolynomialRing, QQ, ZZ
 from sage.symbolic.all import I, pi, SR
 
-# TODO: remove the following line once 4ti2 functions are removed
-path_to_zsolve = os.path.join(SAGE_LOCAL, 'bin', 'zsolve')
-
 
 
 def _sandpile_help(cls, usage, verbose=True):
@@ -5174,7 +5171,7 @@ class SandpileDivisor(dict):
         sign_file.close()
         # compute
         try:
-            os.system(path_to_zsolve+' -q ' + lin_sys + ' > ' + lin_sys_log)
+            os.system('zsolve -q ' + lin_sys + ' > ' + lin_sys_log)
             # process the results
             zhom_file = open(lin_sys_zhom,'r')
         except IOError: