%
% Display Confidence Intervals (raw scale)
%
head=[];
if Cov_Mat==1
  head=[head 'PCA on Covariance Matrix, '];
else
  head=[head 'PCA on Correlation Matrix, '];
end
if Normal_Dist==1
  head = [head 'Multivariate Normality Assumed'];
else
  head = [head 'Multivariate Normality Not Assumed'];
end
headb=['  psi_hat   SE(psi_hat)   Lower    Upper '];
form=[' %6.2f     %6.2f       %6.2f   %6.2f'];
disp(head)
disp([' '])
disp(['lambda_hat'])
disp(lam')
disp(['Gamma_hat'])
disp(Gam)
disp(['Confidence Intervals for Eigenvalues'])
disp(headb)
for i=1:size(C,2)
  [Str, errmsg]=sprintf(form, CI_eigvalue(i,:));
  disp(Str)
end  
disp([' '])
disp(['Confidence Intervals for Partial Sum of Eigenvalues/trace(Sigma)'])
disp(headb)
for i=1:size(C,2)
  [Str, errmsg]=sprintf(form, CI_ratio(i,:));
  disp(Str)
end  

%
% Display Confidence Intervals (log scale)
%
disp([' '])
disp(['Confidence Intervals for Eigenvalues (Using Log Transformation)'])
disp(headb)
for i=1:size(C,2)
  [Str, errmsg]=sprintf(form, CI_eigvalue_log(i,:));
  disp(Str)
end  
disp([' '])
disp(['Confidence Intervals for Partial Sum of Eigenvalues/trace(Sigma)'])
disp(['(Using Log Transformation)'])
disp(headb)
for i=1:size(C,2)
  [Str, errmsg]=sprintf(form, CI_ratio_log(i,:));
  disp(Str)
end