One method for the calculation of current density distribution in a finite number of long parallel conductors, having rectangular cross section, is proposed in this paper. Numerical results aim to highlight the importance of the skin effect, which can be combined with the proximity effect. The method of superposition of these two effects was applied to the calculation of the electromagnetic field in electric power busbars systems. It has been shown that the skin effect has a much greater impact, especially when the conductors are thin and strong electric currents flow through them, so special attention is paid to its calculation. For numerical solution the integral equations are used. The function of current density is approximated by the finite functional series. This way leads to a very accurate solution with only two terms. Differential evolution method is applied for minimization of error function. To demonstrate the application of the proposed approach, numerical values for busbars are presented and compared with values obtained by using the finite elements method.