Implement (not completely, not tested) the rest of Kuzmin's slope limiting for quadratic elements
	- use it in the example (works).
Make some Solution - methods available outside.
Add centroids of elements to global.h.