We introduce our main results on solving a nonlinear steady-state heat conduction problem in anisotropic and nonhomogeneous media and its finite element approximation. In particular, we concentrate on temperature distribution profile typical e.g. for large transformers and rotary machines. Several theorems on the existence and uniqueness of weak and discrete solutions and their convergence are presented.