diff --git a/examples/images_contours_and_fields/interpolation_methods.py b/examples/images_contours_and_fields/interpolation_methods.py
new file mode 100644
index 000000000000..eda8cfa7efee
--- /dev/null
+++ b/examples/images_contours_and_fields/interpolation_methods.py
@@ -0,0 +1,32 @@
+'''
+Show all different interpolation methods for imshow
+'''
+
+import matplotlib.pyplot as plt
+import numpy as np
+
+# from the docs:
+
+# If interpolation is None, default to rc image.interpolation. See also
+# the filternorm and filterrad parameters. If interpolation is 'none', then
+# no interpolation is performed on the Agg, ps and pdf backends. Other
+# backends will fall back to 'nearest'.
+#
+# http://matplotlib.org/api/pyplot_api.html#matplotlib.pyplot.imshow
+
+methods = [None, 'none', 'nearest', 'bilinear', 'bicubic', 'spline16',
+           'spline36', 'hanning', 'hamming', 'hermite', 'kaiser', 'quadric',
+           'catrom', 'gaussian', 'bessel', 'mitchell', 'sinc', 'lanczos']
+
+grid = np.random.rand(4, 4)
+
+fig, axes = plt.subplots(3, 6, figsize=(12, 6),
+                         subplot_kw={'xticks': [], 'yticks': []})
+
+fig.subplots_adjust(hspace=0.3, wspace=0.05)
+
+for ax, interp_method in zip(axes.flat, methods):
+    ax.imshow(grid, interpolation=interp_method)
+    ax.set_title(interp_method)
+
+plt.show()