Initial data for the exact solution subroutine d_const/u_exa.f: 3. Initial temperature 'u0' in the center of the bar. Copyright 1996 Leon van Dommelen