In this paper we investigate the problem of recovering the source f in the elliptic system −∇ · α∇u + βu = f in Ω, α∇u · n + σu = j on ∂Ω from an observation z of the state u on a part Γ of the boundary ∂Ω, where the functionals α, β, σ and j are given. For the particular interest in reconstructing probably discontinuous sources, we use the standard least squares method with the total variation regularization, i.e. we consider a minimizer of the minimization problem 2 = O 1 n is satisfied. Finally, a numerical experiment is presented to illustrate our theoretical findings.