% omega_ratio=input(' Enter (omega22/omega11): ');
%
omega_ratio=1;
%
nf=301;
df=0.8/(nf-1);
%
arg=zeros(nf,1);
tf=zeros(nf,1);
%
for i=1:nf
%    
     mu(i)=(i-1)*df;
%
     a=1;
     b=-2-mu(i);
     c=1;
%     
     alpha21=(-b+sqrt(b^2-4*a*c))/(2*a);
     alpha22=(-b-sqrt(b^2-4*a*c))/(2*a);
%
     a1(i)=+sqrt(alpha21);
     a3(i)=+sqrt(alpha22);
%
     out1=sprintf(' %8.4g %8.4g %8.4g ',mu(i),a1(i),a3(i));
     disp(out1);
%
end
%
Q1=[a1',mu'];
Q2=[a3',mu'];