diff -up src/sage/sandpiles/sandpile.py.orig src/sage/sandpiles/sandpile.py
--- src/sage/sandpiles/sandpile.py.orig	2018-08-03 05:10:09.000000000 -0600
+++ src/sage/sandpiles/sandpile.py	2018-08-29 08:47:29.728055776 -0600
@@ -342,9 +342,6 @@ from sage.arith.all import falling_facto
 from sage.rings.all import Integer, PolynomialRing, QQ, ZZ
 from sage.symbolic.all import I, pi
 
-# TODO: remove the following line once 4ti2 functions are removed
-path_to_zsolve = os.path.join(SAGE_LOCAL,'bin','zsolve')
-
 class Sandpile(DiGraph):
     """
     Class for Dhar's abelian sandpile model.
@@ -5197,7 +5194,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: