function s=sumapoly(p1,p2)
n=length(p1);
m=length(p2);
if n>m
    p2=[zeros(1,n-m) p2];
else
    p1=[zeros(1,m-n) p1];
end
s=p1+p2;