function q=polyint(y,y0)      %  lowest power on right
% { ************************************************************ }
% { *  Integrates polynomial y,                                * }
% { *  y0 represents the integration constant                  * }
% { ************************************************************ }

n=max(size(y));                        
poly=fliplr(y);
for i=1:n, poly(i)=poly(i)/i; end
poly=[y0 poly];
q=fliplr(poly);

    


