diff --git a/test_separability.py b/test_separability.py
new file mode 100644
index 000000000..8c52e09ff
--- /dev/null
+++ b/test_separability.py
@@ -0,0 +1,13 @@
+from astropy.modeling import models as m
+from astropy.modeling.separable import separability_matrix
+
+cm = m.Linear1D(10) & m.Linear1D(5)
+
+print("separability_matrix(cm)")
+print(separability_matrix(cm))
+
+print("separability_matrix(m.Pix2Sky_TAN() & m.Linear1D(10) & m.Linear1D(5))")
+print(separability_matrix(m.Pix2Sky_TAN() & m.Linear1D(10) & m.Linear1D(5)))
+
+print("separability_matrix(m.Pix2Sky_TAN() & cm)")
+print(separability_matrix(m.Pix2Sky_TAN() & cm))