function K = polyK(x1,x2,p,b)
  if nargin<4
      b=1;
  end
  if nargin<3
      p=1;
  end
  if nargin<2
      x2=x1;
  end
  K = (b+x1*x2').^p;
  