tNG{i}(k)] = ...
</span><span class="comment">            min(.5*(z(k)-x).^2+normalGammaNeglogpdf(x, delta, 1)');
</span><span class="comment">    end
</span><span class="comment">    %}
</span>

<span class="comment">    %{
</span><span class="comment">        %% normal inverse Gaussian
</span><span class="comment">        deltas = [0.01 0.75 1 2];
</span><span class="comment">        for i=1:length(deltas)
</span><span class="comment">            delta = deltas(i);
</span><span class="comment">           