Algorithme de Durand-Kerner degré 3

On peut par itération trouver les racines d'un polynôme. Ici de degré 3, défini par ses racines [math]a_1[/math], [math]a_2[/math] et [math]a_3[/math]. On part de trois valeurs complexes [math]z_1[/math], [math]z_2[/math] et [math]z_3[/math] et on itère le système dynamique défini par[br][math]z'_j=z_j-\frac{P(z_j)}{(z_j-z_k)(z_j-z_\ell)}[/math][br][br]Cette méthode est dûe à [url=https://en.wikipedia.org/wiki/Durand–Kerner_method]Durand-Kerner[/url].
Modifiez la position du point [math]z_1[/math] et observez comme son itération converge vers l'une ou l'autre des racines, partitionnant le plan en l'ensemble de Julia du système dynamique.[br][br]Vous pouvez également modifier la position des racines [math]a_j[/math] , ainsi que la position des autres points de départ [math]z_2[/math] et [math]z_3[/math].

Information: Algorithme de Durand-Kerner degré 3