5 5];
    bs = [1 0.1 0.05 0.1];
    
    as = [0.01 0.75 1]; % 2];
    bs = ones(1,4);
    
    for i=1:length(as)
        [temp outNEG{i}(k)]=min(.5*(z(k)-x).^2+normalExpGammaNeglogpdf(x, as(i), bs(i)));
    end
    %}
end
%% Plot
fs = 14;

[styles, colors, symbols] =  plotColors();
figure();
plot(z, x(outLap), [styles{1}, colors(1)], 'linewidth', 3);
hold on
plot(z, z, [styles{2}, colors(2)], 'linewidth', 3)
title('Lasso', 'fontsize', fs)
xlabel('w^{MLE}', 'fontsize', fs)
ylabel('w^{MAP}', 'fontsize', fs