import numpy as np
from scipy.optimize import root_scalar

def f(x, *my_arg, **my_kwargs): # note dummy variable needed by root_scalar
    return np.tanh(2.8 - x)

