return np.prod(s / (s + 1) ** 2, axis=0)