diff --git a/sumpy/expansion/__init__.py b/sumpy/expansion/__init__.py
index e1da810e632822e67e167470e1f02e87538b5b67..00849b56d1746e02e1f9aa84dd0753223dd4e933 100644
--- a/sumpy/expansion/__init__.py
+++ b/sumpy/expansion/__init__.py
@@ -378,7 +378,7 @@ class LinearRecurrenceBasedExpansionTermsWrangler(ExpansionTermsWrangler):
                         c = tuple(pde_mi2 + diff)
                         if c not in coeff_ident_enumerate_dict:
                             break
-                        eq[coeff_ident_enumerate_dict[c]] = 1
+                        eq[coeff_ident_enumerate_dict[c]] = coeff
                     else:
                         pde_mat.append(eq)
         if len(pde_mat) > 0: