def voigt(x, a, b, c, d): z = (x-b+1j*d)/c/numpy.sqrt(2) result = a*numpy.real(wofz(z))/c/numpy.sqrt(2*numpy.pi) return result