For the first lime in electro-thermal device and circuit simulation thermal nonlinearity due to temperature dependent diffusivity is fully treated. It is shown that transformation of the time variable must be employed, in combination with thr well-known Kirchhoff transformation for temperature-dependent conductivity to solve accurately the time-dependent heat diffusion equation.