[en] A self-consistent kinetic treatment is presented here, where the Boltzmann equation is solved for a particle conserving Krook collision operator. The resulting equations have been implemented numerically. The treatment solves for the entire quasineutral column, making no assumptions about λ_mfp/L where λ_mfp is the ion-neutral collision mean free path and L is the size of the device. Coulomb collisions are neglected in favour of collisions with neutrals, and the particle source is modeled as a uniform Maxwellian. Electrons are treated as an inertialess but collisional fluid. The ion distribution function for the trapped and the transiting orbits is obtained. Interesting findings include the anomalous heating of ions as they approach the presheath, the development of strongly non-Maxwellian features near the last λ_mfp, and strong modifications of the sheath criterion. (author)