* cosmo.el (cosmo-trapzd, cosmo-qsimp): Adapt Simpson's rule to
automatically check convergence. At similar accuracy, the code is
about one order of magnitude slower than the previous
implementation. Try to improve efficiency by replacing loops with
more functional solutions. For the time being the newly added code
is commented out, the old implementation is used.