Add option to save kernel density for posterior objects

This allows for example getting the posterior density of forecasts.

It simply requires adjusting pm3.m to save the objects that can already be computed via posterior_moments. I would propose to call the option posterior_kernel_density and save it to options_.estimation.posterior_kernel_density.indicator=1