UDC 537.87We use the well-known method of analysis of oversized electrodynamic systems, which is based on numerical solution of the integral equation for the surface current and is employed to calculate diffraction of an electromagnetic field by arbitrary conducting surfaces, for calculation of the characteristics of waveguide structures. In contrast to the finite-difference algorithms for straightforward solution of the Maxwell equations, which are dealt with the field components in the calculation space divided into cells with small sizes (compared with the wavelength), the desired quantity in this method is the current on the conducting surface bounding the structure. This reduces the required RAM capacity and speeds up calculations by several orders of magnitude in the case of oversized systems. In this work, the integral equation method is realized in the form of a computer code and is used to calculate three-dimensional multimode waveguide systems. The calculation results are compared with the results of perturbation theory and the experimental data.