function q = polyint(p,a,b)
  %POLYINT 多項式積分快速算法
  %polyint(p)返回多項式p的不定積分,常數(shù)項用0;
  %polyint(p,a,b) 返回多項式p在[a,b]的定積分.
  
  %by L.J. Hu 1999
  
  p=p(:)';
  n=length(p);
  q=p./[n:-1:1];
  q=[q,0];
  if nargin>2,
     q=polyval(q,b)-polyval(q,a);
  end