Una observación preliminar: La ecuación $f(x)=0$ equivale a $\tan(x)=\frac1x$ . Ahora, $\tan$ es $\pi$ -periódica y continua, por lo que para grandes $x$ las soluciones serán arbitrariamente cercanas a las soluciones de $\tan(x)=0$ que es sólo $x=\pi n$ para $n\in\mathbb Z$ . Por eso (casi) todas sus soluciones se acercan a los múltiplos de $\pi$ .
Utilizando el método de Newton podemos dar las siguientes aproximaciones para la siguiente raíz:
Si $x_0$ es una solución a $f(x_0)=0$ , entonces la siguiente solución $x_1$ es aproximadamente $$x_1=h_n(x_0+\pi),$$
donde $h_n$ denota $n$ aplicaciones del método de Newton. Por ejemplo (perdón por lo feo) $h_3$ expresión), $$h_1(x)=x-\frac{f(x)}{f'(x)}=x-\frac{x \tan (x)-1}{\tan (x)+x \sec ^2(x)}$$ y
$$h_2(x)=x-\frac{x \tan (x)-1}{\tan (x)+x \sec ^2(x)}-\frac{\left(x-\frac{x \tan (x)-1}{\tan (x)+x \sec ^2(x)}\right) \tan \left(x-\frac{x \tan (x)-1}{\tan (x)+x \sec ^2(x)}\right)-1}{\left(x-\frac{x \tan (x)-1}{\tan (x)+x \sec ^2(x)}\right) \sec ^2\left(x-\frac{x \tan (x)-1}{\tan (x)+x \sec ^2(x)}\right)+\tan \left(x-\frac{x \tan (x)-1}{\tan (x)+x \sec ^2(x)}\right)}$$ y $$h_3(x)=x+\frac{1-x \tan (x)}{x \sec ^2(x)+\tan (x)}-\frac{\frac{\left(x^2 \sec ^2(x)+1\right) \tan \left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}{x \sec ^2(x)+\tan (x)}-1}{\frac{\left(x^2 \sec ^2(x)+1\right) \sec ^2\left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}{x \sec ^2(x)+\tan (x)}+\tan \left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}-\frac{\left(x+\frac{1-x \tan (x)}{x \sec ^2(x)+\tan (x)}-\frac{\frac{\left(x^2 \sec ^2(x)+1\right) \tan \left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}{x \sec ^2(x)+\tan (x)}-1}{\frac{\left(x^2 \sec ^2(x)+1\right) \sec ^2\left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}{x \sec ^2(x)+\tan (x)}+\tan \left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}\right) \tan \left(x+\frac{1-x \tan (x)}{x \sec ^2(x)+\tan (x)}-\frac{\frac{\left(x^2 \sec ^2(x)+1\right) \tan \left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}{x \sec ^2(x)+\tan (x)}-1}{\frac{\left(x^2 \sec ^2(x)+1\right) \sec ^2\left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}{x \sec ^2(x)+\tan (x)}+\tan \left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}\right)-1}{\left(x+\frac{1-x \tan (x)}{x \sec ^2(x)+\tan (x)}-\frac{\frac{\left(x^2 \sec ^2(x)+1\right) \tan \left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}{x \sec ^2(x)+\tan (x)}-1}{\frac{\left(x^2 \sec ^2(x)+1\right) \sec ^2\left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}{x \sec ^2(x)+\tan (x)}+\tan \left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}\right) \sec ^2\left(x+\frac{1-x \tan (x)}{x \sec ^2(x)+\tan (x)}-\frac{\frac{\left(x^2 \sec ^2(x)+1\right) \tan \left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}{x \sec ^2(x)+\tan (x)}-1}{\frac{\left(x^2 \sec ^2(x)+1\right) \sec ^2\left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}{x \sec ^2(x)+\tan (x)}+\tan \left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}\right)+\tan \left(x+\frac{1-x \tan (x)}{x \sec ^2(x)+\tan (x)}-\frac{\frac{\left(x^2 \sec ^2(x)+1\right) \tan \left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}{x \sec ^2(x)+\tan (x)}-1}{\frac{\left(x^2 \sec ^2(x)+1\right) \sec ^2\left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}{x \sec ^2(x)+\tan (x)}+\tan \left(\frac{2 x^2+\cos (2 x)+1}{2 x+\sin (2 x)}\right)}\right)}$$
2 votos
¿Qué hay de malo en utilizar Método de Newton ¿y sólo ajustar los puntos de inicio?
0 votos
¿No es demasiado complicado aplicar el método de Newton una y otra vez, cuando las raíces siguen un patrón predecible? La distancia entre las raíces aumenta, y la velocidad a la que aumenta esa distancia disminuye. Parece que la distancia entre raíces converge a $2\pi$ . ¿Podría explotarse esto de algún modo para obtener una fórmula de la distancia a la siguiente raíz, o de la posición de la siguiente raíz?
0 votos
Parece que esto es más que nada un duplicado de su pregunta anterior aquí , sólo aplicando algunas reducciones señaladas en los comentarios; no es tan bueno tener estas dos preguntas, ya que podría suceder fácilmente que una recibiera respuesta, y luego la otra recibiera una respuesta casi idéntica de alguien que no conocía la solución anterior.
0 votos
@MiloBrandt Tienes razón, me he adelantado y he borrado la pregunta anterior.
2 votos
@Kos A pesar de tu preocupación voy a publicar una respuesta a tu pregunta utilizando el método de Newton. Es un poco feo pero funciona bien
0 votos
Gracias, lo estoy deseando.
2 votos
Pero el comentario sobre el método de Newton podría ser un buen comienzo, que eligiendo $x_0+\pi$ como conjetura inicial, entonces una aproximación podría ser $$x_1 = x_0 + \pi - \frac{f(x_0+\pi)}{f'(x_0+\pi)} = x_0 +\pi+g(x_0)$$ donde $g(x) = - \dfrac{f(x+\pi)}{f'(x+\pi)}$ .