Some stuff

src/funcs.cpp

@@ -4,6 +4,10 @@ double f(double x) {
     return 100*std::exp(-10*x);
 }
 
+double u(double x) {
+    return 1. - (1. - std::exp(-10.))*x - std::exp(-10.*x);
+}
+
 void build_g_vec(int n_steps, arma::vec& g_vec) {
     g_vec.resize(n_steps-1);