A discrete sampling inversion scheme for the heat equation (Q1104062)

A simple and extremely accurate procedure is presented for approximating the initial temperature for the heat equation on the line using a discrete time and spatial sampling. The procedure is based on the ``sinc expansion'' which for functions in a particular class yields a uniform exponential error bound with exponent depending on the number of spatial sample locations chosen. Further the temperature needs only be sampled at one and the same temporal value for each of the spatial sampling points. For N spatial sample points, the approximation is reduced to solving a linear system with \(a(2N+1)\times (2N+1)\) coefficient matrix. This matrix is a symmetric centrosymmetric Toeplitz matrix and hence can be determined by computing only \(2N+1\) values using quadratures.