x(outNEG{i}),[styles{i+2}, colors(i+2)], 'linewidth', 3);
    str{i} = sprintf('a = %5.3f, b=%5.3f', as(i), bs(i));
end
legend(str, 'location', 'southeast')
plot(z, z, [styles{2}, colors(2)], 'linewidth', 3)
title('normalExpGamma', 'fontsize', fs)
xlabel('w^{MLE}', 'fontsize', fs)
ylabel('w^{MAP}', 'fontsize', fs)
printPmtkFigure('NEGthreshold')
%}
%%

end

%%
function out = pen_NIG(w, delta, gamma)
tmp = sqrt(delta^2 + w.^2);
out = log(tmp)-log(besselk(1,gamma*tmp));
end

##### SOURCE END #####
--></body><