Harmonic interpolating wavelets in Neumann boundary value problem in a circle

The Neumann boundary value problem (BVP) in a unit circle is discussed. For the solution of the Neumann BVP, we built a method employing series representation of given $2 \pi$-periodic continuous boundary function by interpolating wavelets consisting of trigonometric polynomials. It is convenient to use the method due to the fact that such series is easy to extend to harmonic polynomials inside a circle. Moreover, coefficients of the series have an easy-to-calculate form. The representation by the interpolating wavelets is constructed by using an interpolation projection to subspaces of a multiresolution analysis with basis $2 \pi$-periodic scaling functions (more exactly, their binary rational compressions and shifts). That functions were developed by Subbotin and Chernykh on the basis of Meyer-type wavelets. We will use three kinds of such functions, where two out of the three generates systems, which are orthogonal and simultaneous interpolating on uniform grids of the corresponding scale and the last one generates only interpolating on the same uniform grids system. As a result, using the interpolation property of wavelets mentioned above, we obtain the exact representation of the solution for the Neumann BVP by series of that wavelets and numerical bound of the approximation of solution by partial sum of such series.