double summ2 = 0.0; for (int i = 1; i <= n2; ++i) { a[i] = ppnd((i - TH) / an25); summ2 += a[i] * a[i]; double ssumm2 = Math.sqrt(summ2); double rsn = 1.0 / Math.sqrt(an); double a1 = poly(C1, 6, rsn) - a[1] / ssumm2; if (n > 5) { i1 = 3; double a2 = -a[2] / ssumm2 + poly(C2, 6, rsn); fac = Math.sqrt((summ2 - 2.0 * a[1] * a[1] - 2.0 * a[2] * a[2]) / (1.0 - 2.0 * a1 * a1 - 2.0 * a2 * a2)); sa += sign(1, i - j) * a[Math.min(i, j)]; for (int i = 1; i <= n1; ++i) { if (i != j) { asa = sign(1, i - j) * a[Math.min(i, j)] - sa; } else { asa = -sa; double s = 1.0; if (n <= 11) { double gamma = poly(G, 2, an); if (y >= gamma) { pw[0] = SMALL;