#include <cassert>
#include <cmath>

// Approximate the following integral using Simpson's rule
//  - this finds the ell-two error between the two functions
//
//                ---------------------
//               /  b
//              /  /              2
//             /   | (f(x) - g(x))  dx
//         \  /    / 
//          \/    a


template <typename T>
T l2error( T f( T x ), T g( T x ), T a, T b,
           unsigned int n, T eps_abs, unsigned int max_iterations );

#include "l2error.tpp"