Mehrfachintegrale

restart:
f:={x,y}->         ;# (* Funktion eingeben *)
fy:=x->            ;# (* Obere Integrationsgrenze für y *)
xu:=               ;# (* Untere Integrationsgrenze für x *)
xo:=               ;# (* Obere Integrationsgrenze für x *)
yu:=               ;# (* Untere Integrationsgrenze für y *)
int(int(f(x,y),y=yu..fy(x)),x=xu..xo);

Beispiel:

restart:
f:={x,y}->x^2+y^2;
fy:=x->x^2;
xu:=1;
xo:=2;
yu:=1;
int(int(f(x,y),y=yu..fy(x)),x=xu..xo);
#(* Darstellung der Integrationsgrenzen *)
g1:=plot([[2,0],[2,5]]):
g2:=plot({fy(x),1},x=0..2.2):
with(plots):display({g1,g2});