from math import sqrt as s

funcs =(
    lambda a,b: (s(a**2 + a * s(a**2 - b**2)) + s(a**2 - a * s(a**2 - b**2))) / (4 * a**3 * b) / (s(a / b) + s(b / a) - 2),
    lambda a,b: (s(a**2 + a * s(a**2 - b**2)) - s(a**2 - a * s(a**2 - b**2))) / (4 * a**3 * b) / (s(a / b) + s(b / a) - 2),
    lambda a,b: (s(a) + s(b)) / (s(a) - s(b)),
)

print([f(4.0, 1.0) for f in funcs])