That interpolating scheme guarantees the accuracy \(O(h(k)^4)\)  ([1]}). We notice, that the order of the interpolating scheme is lower than the numerical integration order \( p=5 \)  and affects the integration order after the discrete transition. Nevertheless, the estimated orders for reduced gradients calculation presented in Table REF  remains valid if we assume integration order \( p=4 \) .
