y = (1 + np.linspace(0, 5, np.prod(ys))).astype(dtype).reshape(ys)