restart: with(linalg): fx:= : fy:= : fz:= : v:=([fx,fy,fz]); xu:= : xo:= : yu:= : yo:= : zu:= : zo:= : int(int(int(diverge(v,[x,y,z]),x=xu..xo),y=yu..yo),z=zu..zo): simplify(");