diff --git a/specs/_features/eip7594/polynomial-commitments-sampling.md b/specs/_features/eip7594/polynomial-commitments-sampling.md
index 9cd862589..0a90443e9 100644
--- a/specs/_features/eip7594/polynomial-commitments-sampling.md
+++ b/specs/_features/eip7594/polynomial-commitments-sampling.md
@@ -196,12 +196,15 @@ def multiply_polynomialcoeff(a: PolynomialCoeff, b: PolynomialCoeff) -> Polynomi
     """
     Multiplies the coefficient form polynomials ``a`` and ``b``
     """
+    assert len(a) + len(b) <= FIELD_ELEMENTS_PER_EXT_BLOB
+
     r = [0]
     for power, coef in enumerate(a):
         summand = [0] * power + [int(coef) * int(x) % BLS_MODULUS for x in b]
         r = add_polynomialcoeff(r, summand)
     return r
 ```
+
 #### `divide_polynomialcoeff`
 
 ```python
diff --git a/tests/core/pyspec/eth2spec/test/eip7594/unittests/polynomial_commitments/test_polynomial_commitments.py b/tests/core/pyspec/eth2spec/test/eip7594/unittests/polynomial_commitments/test_polynomial_commitments.py
index 0b60cfd28..f8a879e46 100644
--- a/tests/core/pyspec/eth2spec/test/eip7594/unittests/polynomial_commitments/test_polynomial_commitments.py
+++ b/tests/core/pyspec/eth2spec/test/eip7594/unittests/polynomial_commitments/test_polynomial_commitments.py
@@ -2,6 +2,7 @@ import random
 from eth2spec.test.context import (
     spec_test,
     single_phase,
+    expect_assertion_error,
     with_eip7594_and_later,
 )
 from eth2spec.test.helpers.sharding import (
@@ -105,3 +106,19 @@ def test_recover_polynomial(spec):
     # Now flatten the cells and check that they match the entirety of the recovered data
     flattened_cells = [x for xs in cells for x in xs]
     assert flattened_cells == recovered_data
+
+
+@with_eip7594_and_later
+@spec_test
+@single_phase
+def test_multiply_polynomial_degree_overflow(spec):
+    rng = random.Random(5566)
+
+    # Perform a legitimate-but-maxed-out polynomial multiplication
+    poly1_coeff = [rng.randint(0, BLS_MODULUS - 1) for _ in range(spec.FIELD_ELEMENTS_PER_BLOB)]
+    poly2_coeff = [rng.randint(0, BLS_MODULUS - 1) for _ in range(spec.FIELD_ELEMENTS_PER_BLOB)]
+    _ = spec.multiply_polynomialcoeff(poly1_coeff, poly2_coeff)
+
+    # Now overflow the degree by pumping the degree of one of the inputs by one
+    poly2_coeff = [rng.randint(0, BLS_MODULUS - 1) for _ in range(spec.FIELD_ELEMENTS_PER_BLOB + 1)]
+    expect_assertion_error(lambda: spec.multiply_polynomialcoeff(poly1_coeff, poly2_coeff))