An estimate of error in $L_p$-approximation in terms of higher order integral modulus of smoothness is obtained using the device of Steklov means for an iterative combination, due to Micchelli, of Phillips operators.