Fractrace.Formulas.H2 C# (CSharp) Method

        long H2(double ar, double ai, double aj, double ak, double br, double bi, double bj, double bk, long zkl)
        {
            double aar, aai, aaj, aak;
            long tw;
            int n;

            aar = ar * ar; aai = ai * ai; aaj = aj * aj; aak = ak * ak; tw = 0L;
            for (n = 1; n < zkl; n++)
            {
                ar = aar - aai - aaj - aak - br;
                ai = 2 * ai * ar - bi;
                aj = 2 * (ar * aj) - bj;
                if (aj > 1) aj = 1;
                if (aj < -1) aj = -1;
                ak = 2 * (ar * ak) - bk;
                aar = ar * ar; aai = ai * ai; aaj = aj * aj; aak = ak * ak;
                if ((aar + aai + aaj + aak) > gr) { tw = n; break; }
            }
            return (tw);
        }