def ghiggs(t): if t<=1: return csqrt(1-1/t)/2. * ( clog((1 + csqrt(1-1/t))/(1 - csqrt(1-1/t)))-pi*1j) else: return csqrt(1/t-1)*casin(csqrt(t))
def fhiggs(t): if t <= 1.: return casin(sqrt(t))**2. else: return -(log((sqrt(t) + sqrt(t - 1.)) / (sqrt(t) - sqrt(t - 1.))) - pi * 1j)**2. / 4.
def fhiggs(t): if t<=1.: return casin(sqrt(t))**2. else: return -(log((sqrt(t) + sqrt(t-1.))/(sqrt(t) - sqrt(t-1.))) - pi*1j )**2./4.
def ghiggs(t): if t <= 1: return csqrt(1 - 1 / t) / 2. * (clog( (1 + csqrt(1 - 1 / t)) / (1 - csqrt(1 - 1 / t))) - pi * 1j) else: return csqrt(1 / t - 1) * casin(csqrt(t))