function compute_arc_length(f, a, b, steps = 1000, epsilon = 0.01)
{
let delta = (b-a)/steps;
let sum = 0.0;
let integral = new RiemannSum(steps);
let arc_length_function = (x) => Math.sqrt(1 + Math.pow( (f(x + epsilon) - f(x)) / epsilon , 2) )
return integral.compute(arc_length_function, a, b);
}