function[poly]=polyadd1(poly1,poly2)
% Copyright 1996 Justin Shriver
% polyadd(poly1,poly2) adds two polynominals possibly of uneven length

if length(poly1)&lt;length(poly2)
  short=poly1;
  long=poly2;
else
  short=poly2;
  long=poly1;
end
mz=length(long)-length(short);
if mz>0
  poly=[zeros(1,mz),short]+long;
else
  poly=long+short;
end
