diff --git a/sumpy/expansion/pde_utils.py b/sumpy/expansion/pde_utils.py
index bd856a622fd281006e2a5eeda5c2c42ddb1a883f..23152e9a24880730dd42f87900c30770e87c5692 100644
--- a/sumpy/expansion/pde_utils.py
+++ b/sumpy/expansion/pde_utils.py
@@ -96,7 +96,7 @@ class PDE(object):
 def laplacian(pde):
     eqs = []
     for j in range(len(pde.eqs)):
-        eqs.append(div(grad(pde[j])).eqs[0])
+        eqs.append(div(grad(pde[j])))
     return PDE(pde.dim, *eqs)
 
 
@@ -118,7 +118,7 @@ def grad(pde):
     for d in range(pde.dim):
         mi = [0]*pde.dim
         mi[d] += 1
-        eqs.append(diff(pde, mi).eqs[0])
+        eqs.append(diff(pde, mi))
     return PDE(pde.dim, *eqs)