Accurate numerical steady-state solutions for a diffused one-dimensional junction diode