  We consider an iterative procedure to solve quasilinear elliptic systems with
$p$-growth. The scheme was first considered by Koshelev in the quadratic case
$p=2$. We present numerical applications as well as applications to higher
regularity properties.
