diff -up src/sage/sandpiles/sandpile.py.orig src/sage/sandpiles/sandpile.py
--- src/sage/sandpiles/sandpile.py.orig	2015-10-11 17:41:27.067777358 -0300
+++ src/sage/sandpiles/sandpile.py	2015-10-11 17:41:57.909778539 -0300
@@ -347,9 +347,6 @@ from sage.rings.arith import falling_fac
 from sage.rings.all import Integer, PolynomialRing, QQ, ZZ, lcm
 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.
@@ -5230,7 +5227,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: