  We apply the derivative expansion of the effective action in the exact
renormalization group equation up to fourth order to the $Z_2$ and $O(N)$
symmetric scalar models in $d=3$ Euclidean dimensions. We compute the critical
exponents $\nu$, $\eta$ and $\omega$ using polynomial expansion in the field.
We obtain our predictions for the exponents employing two regulators widely
used in ERG computations. We apply Wynn's epsilon algorithm to improve the
predictions for the critical exponents, extrapolating beyond the
next-to-next-to-leading order prediction of the derivative expansion.
