In this paper, the Cauchy problem for the Helmholtz equation, also known as the continuation problem, is considered. The continuation problem is reduced to a boundary inverse problem for a well-posed direct problem. A generalized solution to the direct problem is obtained and an estimate of its stability is given. The inverse problem is reduced to an optimization problem solved using the gradient method. The convergence of the Landweber method with respect to the functionals is compared with the convergence of the Nesterov method. The calculation of the gradient in discrete form, which is often used in the numerical solutions of the inverse problem, is described. The formulation of the conjugate problem in discrete form is presented. After calculating the gradient, an algorithm for solving the inverse problem using the Nesterov method is constructed. A computational experiment for the boundary inverse problem is carried out, and the results of the comparative analysis of the Landweber and Nesterov methods in a graphical form are presented.

This research presents a comparative analysis of numerical methods for solving first-kind Fredholm integral equations using the Bubnov–Galerkin method with various wavelet and orthogonal polynomial bases. The bases considered are constructed from Legendre, Laguerre, Chebyshev, and Hermite wavelets, as well as Alpert multiwavelets and CAS wavelets. The effectiveness of these bases is evaluated by measuring errors relative to known analytical solutions at different discretization levels. Results show that global orthogonal systems—particularly the Chebyshev and Hermite—achieve the lowest error norms for smooth target functions. CAS wavelets, due to their localized and oscillatory nature, produce higher errors, though their accuracy improves with finer discretization. The analysis has been extended to incorporate perturbations in the form of additive noise, enabling a rigorous assessment of the method’s stability with respect to different wavelet bases. This approach provides insight into the robustness of the numerical scheme under data uncertainty and highlights the sensitivity of each basis to noise-induced errors.

In this research, we consider the inverse problem for the wave equation under an unknown initial condition. A generalized solution to the direct problem was formulated, its correctness was established, and the stability assessment was obtained. The inverse problem was reduced to an optimization problem, where the objective function was minimized using gradient methods, including the accelerated Nesterov algorithm. The conjugate problem was constructed, and the functional gradient was computed, while the existence of the Frechet derivative was proved. For the first time, the quaternion Fourier transform (QFT) was applied to the numerical solution of a direct problem, making it possible to analyze multidimensional wave processes more efficiently. A computational experiment was carried out, which demonstrated that if there is insufficient additional information, the restoration of the initial condition is incomplete. The introduction of the second boundary condition makes it possible to significantly improve the accuracy and stability of the solution. The results confirm the importance of an integrated approach and the availability of sufficient a priori information when solving inverse problems.

This article addresses a two-dimensional inverse problem of magnetotelluric sounding under the assumption of E-polarized electromagnetic fields. The main focus is on the construction of a forward numerical model based on the Helmholtz equation with a complex coefficient, and the recovery of electrical conductivity from boundary measurements. The second-order finite difference method is employed for numerical simulation, providing stable approximations of both the direct and the conjugate problems. The inverse problem is formulated as a minimization of a data misfit functional, and solved using Nesterov’s accelerated gradient descent method, which ensures fast convergence and robustness to noise. Numerical experiments are presented for a synthetic model featuring a smooth background conductivity and a localized anomaly. Comparison between the exact and reconstructed solutions demonstrates the high accuracy and efficiency of the proposed algorithm. The developed approach can serve as a foundation for constructing practical inversion schemes in geophysical exploration problems.