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