Algorithme de Durand-Kerner degré 4

On peut par itération trouver les racines d'un polynôme. Ici de degré 4, défini par ses racines [math]a_j[/math]. On part de quatre valeurs complexes [math]z_j[/math] et on itère le système dynamique défini par[br][math]z'_j=z_j-\frac{P(z_j)}{\prod_{k\not =j}(z_j-z_k)}[/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 de l'autre point de départ [math]z_k[/math].

Information: Algorithme de Durand-Kerner degré 4