Sep-19-2019, 02:09 AM
This problem can be solved in many ways, e.g. consider the following:
from scipy.optimize import root_scalar root_scalar(lambda x: func(x, *popt) - your_y_value, bracket=[0, 10000], method='bisect')