private static double polevl(double x, double[] coef, int N) { double ans; ans = coef[0]; for (int i = 1; i <= N; i++) { ans = ans*x + coef[i]; } return ans; }