riemann_sum_1d_plot function

1D Riemann sums with optional plot