Merge pull request matplotlib#19534 from ianthomas23/18986_fill_between_nan

BUG: fill_between with interpolate=True and NaN.

Author: jklymak

lib/matplotlib/axes/_axes.py

@@ -5286,7 +5286,8 @@ def _fill_between_x_or_y(
         where = where & ~functools.reduce(
             np.logical_or, map(np.ma.getmask, [ind, dep1, dep2]))
 
-        ind, dep1, dep2 = np.broadcast_arrays(np.atleast_1d(ind), dep1, dep2)
+        ind, dep1, dep2 = np.broadcast_arrays(
+            np.atleast_1d(ind), dep1, dep2, subok=True)
 
         polys = []
         for idx0, idx1 in cbook.contiguous_regions(where):