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.
PalmW. J.III, Introduction to Matlab 6 for Engineers (McGraw-Hill, Boston, 2001).
2.
GummelH. K., ‘A self-consistent iterative scheme for one-dimensional steady state transistor calculations’, IEEE Trans. Electron Devices, 11 (1964), 455–465.
3.
SzeS. M., Semiconductor Devices Physics and Technology, 2nd edn (John Wiley, New York, 2002).
