Add Poisson 3D solver
Update the solver to use 3D
Implement function integral
Fix a comment
Implement get_rhs
use get_rhs
Simplify a bit
Fix a bug
Polish printing