The distinct advantage of the proposed method is its ability to
calculate the RF of any order and propagation mode in a
straightforward and fast-converging way, taking fully into account
the frequency dependence of the TL configuration, the effect of
TL transpositions or cross-bondings, the lumped elements
connected to the TL terminals and the influence of the imperfect
earth. In cases where only one RF calculation is required, the QR
diagonalisation algorithm is selected, enhancing the efficiency of
the numerical technique and reducing the total computational time
for the pul parameters. On the other hand, in studies where the
calculation of multiple RFs is needed, the LM diagonalisation
algorithm is implemented, leading to a lower total number of
required frequency points. Both numerical approaches guarantee a
stable solution in few iterations and present smaller execution
times compared to other traditional RF calculation methods,
especially in cases involving complex lossy earth representations.