function [l,li,r,ri,f,fi,x,xi]=genpai(gen)
gen=sort(gen);genn=length(gen);l=0;li=0;r=0;ri=0;f=0;fi=0;x=0;xi=0;
for i=1:genn
    if abs(real(gen(i)))<0.00001
       gen(i)=gen(i)-real(gen(i));
    end
    if abs(imag(gen(i)))<0.00001
       uuuu=real(gen(i));
       gen(i)=uuuu;
   end
end
%gen
for i=1:genn
   if abs(real(gen(i)))<0.00001&abs(imag(gen(i)))<0.00001           %判num的零根，sumzi0：零根在z中的下標；%lingzi0i：零根的個數；
      l(li+1)=i;li=li+1;
   else if abs(real(gen(i)))>=0.00001&abs(imag(gen(i)))<0.00001  %判num的實數根，sumzi0：實數根在z中的下標；lingzi0i：實數根的個數；
         r(ri+1)=i;ri=ri+1;
   else if abs(real(gen(i)))>=0.00001&abs(imag(gen(i)))>=0.00001  %判num的復數根，sumzi0：復數根在z中的下標；lingzi0i：復數根的個數；
         f(fi+1)=i;fi=fi+1;
      else
         x(xi+1)=i;xi=xi+1;                  %虛根
      end
   end
end 
end

if li>=1
   for i=1:li
      l(i)=0;                                           %l：零根矩陣
   end
else
   l=[];
end
if ri>=1
   for i=1:ri
      r(i)=real(gen(r(i)));                                           %l：實數根矩陣
   end
   else
   r=[];
end
if fi>=1
   for i=1:fi
      f(i)=gen(f(i));                                           %l：復數根矩陣
   end
   else
   f=[];
end
if xi>=1
   for i=1:xi
      xu(i)=imag(gen(x(i)));                                           %l：虛根矩陣
   end
   x=j*xu;
else
   x=[];
end
l=sort(l);r=sort(r);