The task of this master thesis was to review a frequency-dependent model of metal oxide surge arrester developed by the IEEE working group 3.4.11. and its implementation within the EMTP – RV software package. ABB's EXLIM T surge arresters have been successfully implemented within the software package. The theoretical basis of the IEEE frequency-dependent surge arrester model is described in the introductory part of the paper and in the second chapter. The stages of development and the final version of the model are described, along with the equations that lead to the values of the parameters of individual model elements. The second chapter also describes the individual elements of the model within the EMTP – RV software package, and among other things, their capabilities and impact on simulations are listed. The third chapter contains catalog data of 28 metal-oxide surge arresters EXLIM T. The tables list the surge arresters according to the maximum system voltage and the surge arrester rated voltage. The maximum residual voltage on the surge arrester in case of current surges with different duration and intensities is stated. The third chapter also includes catalog data on each surge arrester dimension as well as illustrations of their shape. The process of modeling a metal oxide surge arrester within the EMTP – RV software package is described in chapter four. A detailed procedure is described in the example of only one surge arrester with a nominal voltage of 444 kV. The procedure starts from the selection of the necessary model elements, followed by the calculation of parameter values and their entry into individual elements. By performing the first iteration of the model, the simulation can be started, and the most accurate model of surge arrester is pursued by adjusting the parameter values of the model elements. For other surge arresters, no detailed procedure is given. Numerical results of the simulations for all surge arresters are shown in a table, as well as the changed parameters of the model. In the fourth chapter, the Heidler function with its parameters is presented. Finally, in the fifth chapter, an analysis of the obtained parameters values of the surge arrester model is performed to develop a method by which the same solution would be achieved using an accelerated procedure in future modeling of surge arresters. The adjustment time of the parameter Vref for the slow wave and the adjustment time of the parameter L1 for the fast wave are successfully shortened by developing equations that connect the specified parameters with the residual voltage on the surge arrester.