model of heat conduction through 1d rod, using partial differential equation approximation