function q=polyad(y1,y2)                  %  lowest power on right
% { ************************************************************ }
% { *  Adds polynomial y1 and y2        ( q=y1+y2 )            * }
% { ************************************************************ }

n1=max(size(y1));  n2=max(size(y2));  n1new=n1;  n2new=n2;

if n1>n2 while n1>n2new, y2=[0 y2]; n2new=n2new+1; end,  end;
if n2>n1 while n2>n1new, y1=[0 y1]; n1new=n1new+1; end,  end;

q=y1+y2;  
