// polyval.c  simple Horners evaluation y = p[0] + p[1]*x + ...p[pwr]x^pwr
double polyval(int pwr, double p[], double x)
{
  int i;
  double y = p[pwr]*x;
  for(i=pwr-1; i>0; i--) y = (p[i]+y)*x;
  return y + p[0];
}