A method for orthonormalizing atomic functions is proposed and substantiated. The efficiency of using a basic system to solve a number of boundary value electrodynamics problems is shown. Numerical experiment and analysis of the physical results showed that the R-functions method combined with basis functions is universal for waveguides and resonators of arbitrary shape.