function [omegabar] = threshold(Y, sigma)
Y = norminv(0.0075);
logombegabar = Y*sigma - 0.5*sigma^2;
omegabar = exp(logomegabar);

% Calibrate steady state probability of default at 3% (see BGG), then fix
% value of sigma and compute corresponding threshold productivity level.