Systems of integral equations with a degenerate kernel and an algorithm for their solution using the Maple program

In the mathematical literature, a scalar integral equation with a degenerate kernel is well described (see below (1)), where all the written functions are scalar quantities). The authors are not aware of publications where systems of integral equations of (1) type with kernels in the form of a product of matrices would be considered in detail. It is usually said that the technique for solving such systems is easily transferred from the scalar case to the vector one (for example, in the monograph A.L. Kalashnikov "Methods for the approximate solution of integral equations of the second kind" (Nizhny Novgorod: Nizhny Novgorod State University, 2017), a brief description of systems of equations with degenerate kernels is given, where the role of degenerate kernels is played by products of scalar rather than matrix functions). However, as the simplest examples show, the generalization of the ideas of the scalar case to the case of integral systems with kernels in the form of a sum of products of matrix functions is rather unclear, although in this case the idea of reducing an integral equation to an algebraic system is also used. At the same time, the process of obtaining the conditions for the solvability of an integral system in the form of orthogonality conditions, based on the conditions for the solvability of the corresponding algebraic system, as it seems to us, has not been previously described. Bearing in mind the wide applications of the theory of integral equations in applied problems, the authors considered it necessary to give a detailed scheme for solving integral systems with degenerate kernels in the multidimensional case and to implement this scheme in the Maple program. Note that only scalar integral equations are solved in Maple using the intsolve procedure. The authors did not find a similar procedure for solving systems of integral equations, so they developed their own procedure.