Added a python routine that analyses neat's output to derive prodability density functions, modes and confidence intervals for all the quantities without resorting to binning. The output needs some work, and the spline fits fail if the same values are repeated many times, but otherwise it seems to work. It can also provide some quick and dirty plots so the PDFs can be inspected. Still needs a bit of polishing, however.