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)); double s = 1.0; if (n <= 11) { double gamma = poly(G, 2, an); if (y >= gamma) { pw[0] = SMALL; m = poly(C3, 4, an); s = Math.exp(poly(C4, 4, an)); } else { m = poly(C5, 4, xx); s = Math.exp(poly(C6, 3, xx)); double z90f = Z90 + bf * Math.pow(poly(C7, 2, Math.pow(XX90, xx)), ld); double z95f = Z95 + bf * Math.pow(poly(C8, 2, Math.pow(XX95, xx)), ld); double z99f = Z99 + bf * Math.pow(poly(C9, 2, xx), ld);