function y = cohesivepoly(x,p)
  y = polyval(p,x);
  