El proceso iterativo en el programa 'fractrace' se muestra a continuación.
Fue necesario añadir un número pequeño pero finito en los denominadores para evitar divergencias en los cálculos numéricos. Ésta no es una solución satisfactoria pues introduce una pequeña deformación que no debe estar allí. Recordemos que el producto de escatores es cerrado en el conjunto $\mathbb{E}_{-}^{1+n}$, de manera que se debe poder resolver este problema numérico de una forma mas satisfactoria.
/* ix - fractrace - tomotrace */
aar = ar * ar; aai = ai * ai; aaj = aj * aj;
double r2 = aar + aai + aaj + aai * aaj /(aar+.001);
double r2 = aar + aai + aaj + aai * aaj /(aar+.001);
for (n = 1; n < zkl; n++) {
r2 = aar + aai + aaj + aai * aaj /(aar+.001);
tar = aar - aai - aaj + aai * aaj /(aar+.001) +br;
tai = 2 * ar * ai - 2 * ai * aaj / (ar+.001) +bi;
taj = 2 * ar * aj - 2 * aj * aai / (ar+.001) +bj;
ar=tar; ai=tai; aj=taj;
aar = ar * ar; aai = ai * ai; aaj = aj * aj;
r2 = aar + aai + aaj + aai * aaj /(aar+.001);
if (r2 > 400) { tw = n; break; }
}
r2 = aar + aai + aaj + aai * aaj /(aar+.001);
tar = aar - aai - aaj + aai * aaj /(aar+.001) +br;
tai = 2 * ar * ai - 2 * ai * aaj / (ar+.001) +bi;
taj = 2 * ar * aj - 2 * aj * aai / (ar+.001) +bj;
ar=tar; ai=tai; aj=taj;
aar = ar * ar; aai = ai * ai; aaj = aj * aj;
r2 = aar + aai + aaj + aai * aaj /(aar+.001);
if (r2 > 400) { tw = n; break; }
}
Ejemplo de archivo .xml para bajar y procesar en 'fractrace': fractrace-formula-example.zip