* Add an R6 class PVfactory that encapsulated Present Value-calculation in arbitrary dimensions (the older calculatePV* functions were hardcoded to a one-dimensional vector!)
* Use the PVfactory instead of calculatePV* functions
* Convert cost PV to use the multi-dimensional PVfactory code
* Implement costs relative to benefits

Implements #80