// polyval.h  simple Horners evaluation y = p[0] + p[1]*x + ...p[pwr]x^pwr
double polyval(int pwr, double p[], double x);