The problem of reconstruction of the Cauchy data for the wave equation in ℝ¹ from the measurements of its solution on the boundary of a finite interval is considered. This is a one-dimensional model for the multidimensional problem of photoacoustics, which was studied previously. The method was adapted and simplified for the one-dimensional situation. The results of numerical testing to see the rate of convergence and the stability of the procedure are given. Some hints are also given on how the procedure of reconstruction can be simplified in the 2d and 3d cases.