f:read("Gib die Funktion f(x):= ein");

a:read("Gib die Untergrenze ein");

b:read("Gib die Obergrenze ein");

n:read("Gib die Streifenzahl ein");

dx:(b-a)/n;

g(k):=x[k]:a+k*dx /* Berechnung der Fußpunkte */;

x_wert:makelist(g(k),k,0,n);

US:sum(f(x[k])*dx,k,0,n-1);

OS:sum(f(x[k])*dx,k,1,n);

F:(US+OS)/2.0;

print("Das bestimmte Integral ist näherungsweise",F);