% This is the root Min-Norm function program.
% the input is the number of source and the othegonal subspace
function Theta=rmnorm(P,Uo)
ij=sqrt(-1);
[L,X]=size(Uo);
c=Uo(1,:)';
c=c/(c'*c);
d=Uo*c;
%
rmno=sort(roots(d)); rmno=rmno(L-1:-1:1);
Theta=sort(asin(angle(rmno(1:P))/pi))';

