diff --git a/diofant/tests/integrals/test_integrals.py b/diofant/tests/integrals/test_integrals.py
index 98b2e33535..a9c268fcd0 100644
--- a/diofant/tests/integrals/test_integrals.py
+++ b/diofant/tests/integrals/test_integrals.py
@@ -1334,3 +1334,13 @@ def test_sympyissue_20360():
 
 def test_sympyissue_20941():
     assert integrate(x**2*sqrt(1 - x**2), (x, 0, 1)) == pi/16
+
+
+@pytest.mark.slow
+def test_sympyissue_21034():
+    f1 = x*(-x**4/asin(5)**4 - x*sinh(x + log(asin(5))) + 5)
+    f2 = (x + cosh(cos(4)))/(x*(x + 1/(12*x)))
+
+    assert (f1.integrate(x).diff(x) - f1).simplify() == 0
+    assert (f2.integrate(x).diff(x) -
+            f2).simplify().rewrite(exp).simplify() == 0
diff --git a/docs/release/notes-0.13.rst b/docs/release/notes-0.13.rst
index d535661da2..f58d3f2539 100644
--- a/docs/release/notes-0.13.rst
+++ b/docs/release/notes-0.13.rst
@@ -43,3 +43,4 @@ These Sympy issues also were addressed:
 * :sympyissue:`20973`: cancel raises PolynomialError for exp(1+O(x))
 * :sympyissue:`20985`: TypeErrors appearing for simple plynomial manipulations (did not happen in v1.6.1)
 * :sympyissue:`21031`: Limit of "limit (((1+x)**(1/x)-(1+2*x)**(1/(2*x)))/asin (x),x,0)" is wrong with v1.7.1
+* :sympyissue:`21034`: (Integration) regressions?