Feature/scalar tests

Added a test for the scalar example, that verifies that it reproduces (approximately) the theoretically expected convergence order.
Also modified the scalar example, so that it fits into the testing framework.
So far, only up to six Lobatto nodes are tested.

Also, the MATCHER used in EXPECT_THAT is a pretty ugly construction and one can probably using something that is much more clever.
Unfortunately, I missed a little point: `Pointwise` expects a *tuple* matcher and `DoubleNear` is not. A workaround would be to get rid of `expected_cr` and put the `EXPECT_THAT` in the same loop where `convrate` is computed:

    EXPECT_THAT(convrate[i], DoubleNear(double(2*nnodes-2), 0.99)) << "Convergence rate at node " << i << " not within expected range";

where `0.99` is the safety factor (somehow it fails with `0.9`).