“…As for the integration of the system of matrix equations, the left-hand sides of which contain a matrix of an nth order, it should be noted that the method of reduction of the considered system to a system of ordinary differential equations of an order of n 2 is the universal method. If the coefficients of the matrix on the right-hand sides are independent of time, the system can be integrated using the Laplace transformation or by passing from matrix (i, j = a, b) to matrix ρ = U U † [16], where matrix U is the matrix of eigenvectors of Hamilton operator (11). We integrated system (13) using the series expansion method the implementation of which can be described as follows.…”