A numerical method for monitoring temperature distribution in which boundary flux and initial state are unknown is presented. Regularizations based on Tikhonov's and Beck's method are employed. And then, regularization parameters are evaluated by L-curve. The method is applied to an actual piping problem in a steam power plant and compared with measured data, and it is also applied to a two-dimensional thermal shock problem.