diff --git a/sklearn/utils/plotting.py b/sklearn/utils/plotting.py
new file mode 100644
index 0000000000000..bcad03b72828a
--- /dev/null
+++ b/sklearn/utils/plotting.py
@@ -0,0 +1,135 @@
+"""
+Plotting utilities.
+"""
+# Authors: Sebastian Raschka <se.raschka@gmail.com>
+#
+# License: BSD 3 clause
+
+
+from itertools import cycle
+import matplotlib
+import matplotlib.pyplot as plt
+import numpy as np
+
+
+def plot_decision_regions(X, y, classifier, X_highlight=None, res=0.02,
+                          cycle_marker=True, legend=1, cmap=None):
+    """
+    Plots decision regions of a classifier.
+
+    Parameters
+    ----------
+    X : array-like, shape = [n_samples, n_features]
+        Feature Matrix.
+    y : array-like, shape = [n_samples]
+        True class labels.
+    classifier : Estimator for classification.
+    X_highlight : array-like, shape = [n_samples, n_features] (default: None)
+        An array with data points that are used to highlight samples in `X`.
+    res : float (default: 0.02)
+        Grid width. Lower values increase the resolution but
+        slow down the plotting.
+    cycle_marker : bool
+        Use a different marker for each class.
+    legend : int
+        Integer to specify the legend location.
+        No legend if legend is 0.
+    cmap : Custom matplotlib colormap object.
+
+    Returns
+    ---------
+    None
+
+    Examples
+    --------
+    from sklearn import datasets
+    from sklearn.svm import SVC
+    iris = datasets.load_iris()
+    X = iris.data[:, [0,2]]
+    y = iris.target
+    svm = SVC(C=1.0, kernel='linear')
+    svm.fit(X,y)
+    plot_decision_region(X, y, classifier=svm, res=0.02, cycle_marker=True, legend=1)
+    plt.xlabel('sepal length [cm]')
+    plt.ylabel('petal length [cm]')
+    plt.title('SVM on Iris')
+    plt.show()
+
+    """
+    # check if data is numpy array
+    for a in (X, y):
+        if not isinstance(a, np.ndarray):
+            raise ValueError('%s must be a NumPy array.' % a.__name__)
+
+    # check if test data is provided
+    plot_testdata = True
+    if not isinstance(X_highlight, np.ndarray):
+        if X_highlight is not None:
+            raise ValueError('X_test must be a NumPy array or None')
+        else:
+            plot_testdata = False
+
+    if len(X.shape) == 2 and X.shape[1] > 1:
+        dim = '2d'
+    else:
+        dim = '1d'
+
+    marker_gen = cycle('sxo^v')
+
+    # make color map
+    colors = ['red', 'blue', 'lightgreen', 'gray', 'cyan']
+    n_classes = len(np.unique(y))
+    if n_classes > 5:
+        raise NotImplementedError('Does not support more than 5 classes.')
+
+    if not cmap:
+        cmap = matplotlib.colors.ListedColormap(colors[:n_classes])
+
+    # plot the decision surface
+
+    if dim == '2d':
+        y_min, y_max = X[:, 1].min() - 1, X[:, 1].max() + 1
+    else:
+        y_min, y_max = -1, 1
+
+    x_min, x_max = X[:, 0].min() - 1, X[:, 0].max() + 1
+    xx, yy = np.meshgrid(np.arange(x_min, x_max, res),
+                         np.arange(y_min, y_max, res))
+
+    if dim == '2d':
+        y_min, y_max = X[:, 1].min() - 1, X[:, 1].max() + 1
+        Z = classifier.predict(np.array([xx.ravel(), yy.ravel()]).T)
+    else:
+        y_min, y_max = -1, 1
+        Z = classifier.predict(np.array([xx.ravel()]).T)
+
+    Z = Z.reshape(xx.shape)
+    plt.contourf(xx, yy, Z, alpha=0.4, cmap=cmap)
+    plt.xlim(xx.min(), xx.max())
+    plt.ylim(yy.min(), yy.max())
+
+    # plot class samples
+
+    for c in np.unique(y):
+        if dim == '2d':
+            y_data = X[y == c, 1]
+        else:
+            y_data = [0 for i in X[y == c]]
+
+        plt.scatter(x=X[y == c, 0],
+                    y=y_data,
+                    alpha=0.8,
+                    c=cmap(c),
+                    marker=next(marker_gen),
+                    label=c)
+
+    if legend:
+        plt.legend(loc=legend, fancybox=True, framealpha=0.5)
+
+    if plot_testdata:
+        if dim == '2d':
+            plt.scatter(X_highlight[:, 0], X_highlight[:, 1],
+                        c='', alpha=1.0, linewidth=1, marker='o', s=80)
+        else:
+            plt.scatter(X_highlight, [0 for i in X_highlight],
+                        c='', alpha=1.0, linewidth=1, marker='o', s=80)