Implement the solver class: - [x] `interface` method: keep the option open to use more solvers - [x] modify the structure to be more "C++ like" (less pointers and list, more std::) - [ ] prototype how to handle the centrifugal term for spherical and cilindrical, 3D and 2D