A numerical method for a singular perturbation problem arising in the modelling of plasma sheaths. (8th January 2008)