For each q∈N0, we construct positive linear polynomial approximation operators Mn that simultaneously preserve k-monotonicity for all 0≤k≤q and yield the estimate
for x∈[0,1] and λ∈[0,2), where and is the second Ditzian–Totik modulus of smoothness corresponding to the “step-weight function” ψ. In particular, this implies that the rate of best uniform q-monotone polynomial approximation can be estimated in terms of .