   %{
    deltas = [0.01 0.75 1];
    for i=1:length(deltas)
        delta = deltas(i);
        [temp outNG{i}(k)] = ...
            min(.5*(z(k)-x).^2+normalGammaNeglogpdf(x, delta, 1)');
    end
    %}
    
    
    %{
        %% normal inverse Gaussian
        deltas = [0.01 0.75 1 2];
        for i=1:length(deltas)
            delta = deltas(i);
            [temp outNIG{i}(k)]=min(.5*(z(k)-x).^2+pen_NIG(x, delta, 1));
        end
    %}
    
    %{
    %% normal exponential Gaussian
    as = [0.5 0.5 0.