1. Introduction The stationary behaviour of semiconductor devices is governed by a set of nonlinear elliptic partial differential equations. This includes a nonlinear Poisson equation and two nonlinear continuity equations. Using Gummel's method [2] we can decouple the nonlinear elliptic system so that at each step we solve an equation of the form
