ENH: support extended axis argument in np.median

makes it more compatible with the np.mean ufunc