Riemann Sums
You can work with partitions.
> gRiemann(x^2+3,x=[-1,-0.5,0.5,1.5,2],[-0.75,0.0,1.1,1.5]);
You can select from standard methods.
> gTrap(x^2+3,x=-1..3,4);
You can even work with "n" rectangles.
> Ln:=Lsum(x^2+3,x=-1..3,n);
Find the limit of the last result.
> limit(Ln,n=infinity);
Find definite or indefinite integrals.
> int(x^2+3,x=-1..3);
Estimate definite integrals given data.
> xd:=[1,2,3,4,5,6];
> yd:=[2.2,3.4,2.7,1.5,2.7,3.2];
> gdiscreteLsum(xd,yd);
Visualize Riemann sums for functions of two variables.
> gMsum3d(x^2+y^2,x=-2..2,y=-2..2,4,4);
>