function [ c ] = polyadd( a, b )
%POLYADD adds two polynomials a and b
%   The two polynomials can be of different degree
%   and both are represented by row vectors

Na=length(a);
Nb=length(b);

if Na>=Nb
    c = a + [zeros(1,Na-Nb),b];
else
    c = [zeros(1,Nb-Na),a] + b;
end

end

