docs: def _separable_approx3(h, N=1):
      """ returns the N first approximations to the 3d function h
      """
      return np.cumsum([np.einsum("i,j,k", fz, fy, fx) for fz, fy, fx in _separable_series3(h, N)], 0)