Newton-Raphson solution of Poisson's equation in a pn diode

This paper presents a numerical solution, using MATLAB, of the electrostatic potential in a pn junction, which obeys Poisson's equation. This numerical method is based on the Newton-Raphson technique and is useful for educational purposes. It can be incorporated in an undergraduate course on semiconductor devices to demonstrate the applications of numerical analysis techniques in device physics. It may also be appropriate for a section on numerical analysis techniques in an engineering mathematics course.