Tengo un conjunto de funcionesf_t con varias raíces (dos en realidad). Quiero encontrar la "primera" raíz y hacer esto confsolve funciona bien la mayor parte del tiempo. El problema es que las dos raíces convergen, ya que t va al infinito. (un simple ejemplo de mis funciones seríaf_t(x) = x^2 - 1/t) Entonces cuanto más grandet obtiene, más erroresfsolve hace. Hay unpredefinido función, similar afsolve a lo que puedo decir que solo debe buscar en un rango determinado (por ejemplo, encontrar siempre la raíz en[0, inf)).

La pregunta es esencialmente la misma quehttps://mathematica.stackexchange.com/questions/91784/how-to-find-numerically-all-roots-of-a-function-in-a-given-range?noredirect=1&lq=1Sin embargo, las respuestas que hay para Mathematica las quiero en Python.

PD: Ahora sé cómo puedo escribir mi propio algoritmo, pero como estos tienden a ser más lentos como los que esperaba, esperaba encontrar uno que hiciera lo mismo. Especialmente he leído esta publicaciónEncuentra la raíz de una función en un intervalo dado