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)