function exact = f(x);
exact = exp(2) / (exp(4) -1) * (exp(2*x) - exp(-2*x)) + x;
