We propose the use of fundamentally distinctive techniques to solve the problem of curve fitting a Mössbauer spectrum. The techniques we investigated are: evolutionary algorithm, basin hopping, and hill climbing. These techniques were applied in isolation and combined to fit different shapes of Mössbauer spectra. The results indicate that complex Mössbauer spectra can be automatically curve fitted using minimum user input, and combination of these techniques achieved the best performance (lowest statistical error). The software and sample of Mössbauer spectra have been made available through a link at the reference.