def feq(a, b, max_relative_error=1e-12, max_absolute_error=1e-12):
    a = float(a)
    b = float(b)
    # if the numbers are close enough (absolutely), then they are equal
    if abs(a - b) < max_absolute_error:
        return True
    # if not, they can still be equal if their relative error is small
    if abs(b) > abs(a):
        relative_error = abs((b - a)/b)
    else:
        relative_error = abs((a - b)/a)
    return relative_error <= max_relative_error
dataset/ETHPy150Open sympy/sympy/sympy/physics/tests/test_hydrogen.py/feq/swapped operands of `-`