In many engineering applications variation in thermal properties is significant especially when there are large temperature gradients or the material is not homogeneous. To account for these variations in thermal properties, develop a finite difference formulation for an internal node in case of a three dimensional steady state heat conduction equation with variable thermal conductivity.