added function to evaluate uh in any point

* we can evaluate the linear combination of base functions (read: the
calculated solution) in any point now with Mesh::eval
* added ePT_GNUPLOT_SURF visualization type that creates a uniform grid
and then plots uh over it using Mesh::eval