Kartesische Koordinaten:
Remove["Global`*"] Needs["Calculus`VectorAnalysis`"] SetCoordinates[Cartesian]; f={fx,fy,fz} fx= ; (* Funktion eingeben *) fy= ; (* Funktion eingeben *) fz= ; (* Funktion eingeben *) xu= ; (* untere Integrationsgrenze für x *) xo= ; (* obere Integrationsgrenze für x *) yu= ; (* untere Integrationsgrenze für y *) yo= ; (* obere Integrationsgrenze für y *) zu= ; (* untere Integrationsgrenze für z *) zo= ; (* obere Integrationsgrenze für z *) Integrate[Div[f],{x,xu,xo},{y,yu,yo},{z,zu,zo}] //Simplify