In many applications, such as the heat conduction and hydrology, there is a need to recover the (possibly discontinuous) diffusion coefficient a from boundary measurements of solutions of a parabolic equation. The complete inverse problem is ill posed and nonlinear, so numerical solution is quite difficult, and we linearize the problem around constant a. We study and solve numerically the linea...