 [temp outNIG{i}(k)]=min(.5*(z(k)-x).^2+pen_NIG(x, delta, 1));
</span><span class="comment">        end
</span><span class="comment">    %}
</span>
<span class="comment">    %{
</span><span class="comment">    %% normal exponential Gaussian
</span><span class="comment">    as = [0.5 0.5 0.5 5];
</span><span class="comment">    bs = [1 0.1 0.05 0.1];
</span><span class="comment">
</span><span class="comment">    as = [0.01 0.75 1]; % 2];
</span><span class="comment">    bs = ones(1,4);
</span><span class="co