% Emacs mode header: -*-Octave-*- source. %

function hlen=show_bc1(omega,q3ar);

%>> q2=logspace(-40,-20,20);
 [x,P,Q,rh,eh,data_ph] = tw_nonlin_4M2(omega,q3ar);
len=length(q3ar);
hlen=floor(len/2);

subplot(1,2,1)
semilogy(1:length(q3ar),q3ar,'-o',1:length(q3ar),data_ph','-o')
grid minor

subplot(1,2,2)
plot(-x,rh(hlen,:),-x,eh(hlen,:))
grid

end
