Add support for setting an initial solution in BVP

Introduce the ability to initialize BVP problems with an optional initial solution. This includes changes to manage `_solution` internally, a new `set_initial_solution` method, and getter for the `solution` property.
 Unused debug print statements were also removed from plasticity.py.