riemann_sum_1d_plot: 1D Riemann sums with optional plot
Description
Computes lower, upper, and midpoint Riemann sums for a scalar function
f(x) on an interval [xmin, xmax]. Optionally draws a
2D plot with rectangles and, if requested, the true curve.