docs: def _separable_approx2(h, N=1):
      """ returns the N first approximations to the 2d function h
      whose sum should be h
      """
      return np.cumsum([np.outer(fy, fx) for fy, fx in _separable_series2(h, N)], 0)