To find the coefficients for bi- and tri-exponential fitting functions to represent the radial dose functions of 16 commercially available brachytherapy sources.The search for the coefficients was done using a genetic algorithm. Coefficients were encoded into chromosomes, which were subjected to crossover and mutation. After each operation, chromosomes were evaluated according to their fitness and the better ones were chosen with higher probability for the next generation. The best chromosomes obtained after 2000 operations were used for the coefficients.For all brachytherapy sources, tri-exponential dose functions agreed with the respective input data within 1.4%. The mean deviation, obtained by averaging absolute deviations of all sources and input data, was <1.0%. For 8 of the 16 sources, the fit offered by bi-exponential functions was virtually identical to that of tri-exponential ones.Tri-exponential functions can accurately represent the radial dose functions of all commercially available brachytherapy sources. For the eight sources where bi-exponential functions provide nearly equally accurate fits, their continued usage is recommended.