matplotlib
diff --git a/‎lib/matplotlib/axes/_axes.py
Lines changed: 6 additions & 28 deletions b/‎lib/matplotlib/axes/_axes.py
Lines changed: 6 additions & 28 deletions
diff --git a/‎lib/matplotlib/tests/baseline_images/test_axes/fill_between_interpolate.pdf
44 Bytes b/‎lib/matplotlib/tests/baseline_images/test_axes/fill_between_interpolate.pdf
44 Bytes
diff --git a/‎lib/matplotlib/tests/baseline_images/test_axes/fill_between_interpolate.png
147 Bytes b/‎lib/matplotlib/tests/baseline_images/test_axes/fill_between_interpolate.png
147 Bytes
@@ -4667,22 +4667,11 @@ def fill_between(self, x, y1, y2=0, where=None, interpolate=False,
                 raise ValueError('Input passed into argument "%r"' % name +
                                  'is not 1-dimensional.')
 
-        if y1.ndim == 0:
-            y1 = np.ones_like(x) * y1
-        if y2.ndim == 0:
-            y2 = np.ones_like(x) * y2
-
         if where is None:
-            where = np.ones(len(x), np.bool)
-        else:
-            where = np.asarray(where, np.bool)
-
-        if not (x.shape == y1.shape == y2.shape == where.shape):
-            raise ValueError("Argument dimensions are incompatible")
+            where = True
+        where = where & ~reduce(np.logical_or, map(np.ma.getmask, [x, y1, y2]))
 
-        mask = reduce(ma.mask_or, [ma.getmask(a) for a in (x, y1, y2)])
-        if mask is not ma.nomask:
-            where &= ~mask
+        x, y1, y2 = np.broadcast_arrays(np.atleast_1d(x), y1, y2)
 
         polys = []
         for ind0, ind1 in mlab.contiguous_regions(where):
@@ -4828,22 +4817,11 @@ def fill_betweenx(self, y, x1, x2=0, where=None,
                 raise ValueError('Input passed into argument "%r"' % name +
                                  'is not 1-dimensional.')
 
-        if x1.ndim == 0:
-            x1 = np.ones_like(y) * x1
-        if x2.ndim == 0:
-            x2 = np.ones_like(y) * x2
-
         if where is None:
-            where = np.ones(len(y), np.bool)
-        else:
-            where = np.asarray(where, np.bool)
-
-        if not (y.shape == x1.shape == x2.shape == where.shape):
-            raise ValueError("Argument dimensions are incompatible")
+            where = True
+        where = where & ~reduce(np.logical_or, map(np.ma.getmask, [y, x1, x2]))
 
-        mask = reduce(ma.mask_or, [ma.getmask(a) for a in (y, x1, x2)])
-        if mask is not ma.nomask:
-            where &= ~mask
+        y, x1, x2 = np.broadcast_arrays(np.atleast_1d(y), x1, x2)
 
         polys = []
         for ind0, ind1 in mlab.contiguous_regions(where):