f(x):=

Breite(a,b,n):=(b-a)/n

Rechteck(xa,ya,xc,yc):=[[xa,ya],[xc,ya],[xc,yc],[xa,yc],[xa,ya]]

U_Hoehe(a,b,n,num):=f(a+num*(b-a)/n)

U_A_Rechteck(a,b,n,num):=U_Hoehe(a,b,n,num)*Breite(a,b,n)

US(a,b,n):=SUM(U_A_Rechteck(a,b,n,num),num,0,n-1)

US_Plot(a,b,n):=VECTOR(Rechteck(a+i*(b-a)/n,0,a+(i+1)*(b-a)/n,f(a+i*(b-a)/n))~
,i,0,n-1)

O_Hoehe(a,b,n,num):=f(a+(num+1)*(b-a)/n)

O_A_Rechteck(a,b,n,num):=O_Hoehe(a,b,n,num)*Breite(a,b,n)

OS(a,b,n):=SUM(O_A_Rechteck(a,b,n,num),num,0,n-1)

OS_Plot(a,b,n):=VECTOR(Rechteck(a+i*(b-a)/n,0,a+(i+1)*(b-a)/n,f(a+(i+1)*(b-a)~
/n)),i,0,n-1)