Examples of the use of the package "Riemann.m".

<<riemann.m

e=2+(1/3)(x-3)^2+Sin[3 x];

PlotRiemann[e,{x,1,5},{AspectRatio -> Automatic},
{Partition -> Uniform}]

PlotRiemann[e,{x,1,5},{AspectRatio -> Automatic},
{Partition -> Random,Height -> Random}]

PlotRiemann[e,{x,1,5},{AspectRatio -> Automatic},

{Partition -> Uniform,Height -> Midpoint,Subdivisions -> 25}]