In this paper, the system of the semiconductor device equations with heat effect is considered. An approximation to the system that makes use of a mixed finite element method
for the electrostatic potential equation combined with a finite element method for the densities
equations and the temperature equation is proposed. Existence and uniqueness of the approximate
solution are proved. A convergence analysis is also given.