Voigt profile fit

Hi,
Would anybody help me to create a program that can find the best-fit profile for a data based on the Voigt function (combined Gaussian or Lorentian)? Is anybody who has created such program already?

Attached is a sample data that I want to fit.

Thanks.