long H1(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++)
{
ai = 2.0 * ar * ai + bi;
aj = 2.0 * ar * aj + bj;
ak = 2.0 * ar * ak + bk;
ar = aar - aai - aaj - aak + br;
aar = ar * ar; aai = ai * ai; aaj = aj * aj; aak = ak * ak;
if ((aar + aai + aaj + aak) > gr) { tw = n; break; }
}
return (tw);
}