The int(f, dom) command computes the integral of the function f over the region specified by dom.

If f is an expression, dom must have the form list(name)=domain where the left-hand side is a list of the variables of integration and domain is an unevaluated function call describing the region. If f is a Vector, the top-level int command is mapped over the components of f with dom passed to int as well as the extra arguments.

The regions of integration are represented by unevaluated function calls. The possible regions of integration are Circle(cen, rad), Ellipse(cen, a, b, phi), Ellipse(eqn), Parallelepiped(r1, ..., rk), Rectangle(r1, r2), Region(r1, ..., rk), Sector(obj, start, finish), Sphere(cen, rad), Tetrahedron(v1, v2, v3, v4), and Triangle(v1, v2, v3).

Circle(cen, rad) The first parameter of Circle, cen, must have type Vector(2, algebraic) and rad must have type algebraic. These represent the center and the radius of the circle, respectively. If a coordinate system attribute is specified on cen, the center is interpreted in that coordinate system.

Ellipse(cen, a, b, phi)

The parameter cen is the center of the ellipse and must have type 'Vector'(algebraic). If a coordinate system attribute is specified on cen, it is interpreted in that coordinate system. The parameters a and b are the lengths of the semimajor and semiminor axes, respectively. The resulting ellipse is constructed via the following process: Start with an ellipse centered at the origin having the specified axes lengths, with its major axis initially on the x-axis. It is rotated through an angle of phi in the counterclockwise direction and translated to cen.

Ellipse(eqn)

The parameter eqn is either an equation specifying the ellipse or an algebraic expression such that the equation eqn = 0 specifies the ellipse. A Cartesian equation for a general conic section is of the form $A{x}^2+Bxy+C{y}^2+Ex+Fy+G=0$, and specifies the locus of all points $⟨x\,y⟩$ that satisfy the equation; this locus is a non-degenerate, real ellipse if and only if three conditions hold:

The two variable names that appear in eqn must be the ones specified on the left-hand side of dom.

Parallelepiped(r1, ..., rk) This specifies a k-fold region of integration. Each ri must have type algebraic..algebraic and the expressions in the ranges cannot depend on the variables of integration. These are the ranges that are used for each variable during integration.

Rectangle(r1, r2) This is a specialization of Parallelepiped to R^2. In the xy-plane, the integration has the order dx dy, and range r1 can be of the form x1(y)..x2(y), in contrast to the restrictions on ranges in Parallelepiped.

Region(r1, ..., rk) This specifies a k-fold region of integration. Each ri must have type algebraic..algebraic and the ith range can only depend on the first i-1 variables. Thus, the region can easily be transformed into an iterated integral.

Sector(obj, start, finish) The first parameter of Sector, obj, must be a Circle or Ellipse structure. The Sector structure defines a sector of the circle or ellipse with boundaries specified by the start and finish angles. The sector is defined by the region inside the circle or ellipse swept out by a ray emanating from its center as it rotates from the start angle to the finish angle.

To define precisely how angles are measured, it suffices to discuss only circles and ellipses centered at the origin. For a circle or ellipse centered elsewhere, the start and finish angles are determined as if the circle or ellipse were first translated to the origin.

For a Circle centered at the origin, angle is measured counterclockwise from the positive x-axis. Therefore, the angle $\frac{3\mathrm{\pi }}{2}$ specifies the negative y-axis.

To define how angle is measured for an Ellipse centered at the origin, we first define the right semimajor axis of the ellipse to be the semimajor axis in the right half-plane (the first and fourth quadrants of the plane). If the major axis of the ellipse is coincident with the y-axis, then its right semimajor axis is defined to be the one on the negative y-axis. Thus, for an ellipse centered at the origin with its major axis sitting on the line y = x, its right semimajor axis is the one inside the first quadrant.

For an Ellipse centered at the origin, angle is measured counterclockwise from its right semimajor axis. Therefore, in the example ellipse given in the previous paragraph, the angle $\frac{\mathrm{\pi }}{4}$ specifies the positive y-axis.

If the start and finish angles are more than $2\mathrm{\pi }$ apart, int does not discard the repeated area. That is, int counts the multiplicity of the region over which the integration is being done. For example, if start = 0, finish = 3Pi, and the region is a Circle structure, then the result is equal to twice the integral of f over the top half of the circle plus its integral over the bottom half of the circle.

If the finish angle is less than the start angle, the sign of the integral is flipped. In general, the integral of f over the sector from start to finish is always the negative of the same integral with the start and finish angles swapped.

Sphere(cen, rad) The first parameter of Sphere, cen, must have type 'Vector'(3, algebraic) and rad must have type algebraic. These represent the center and the radius of the sphere, respectively. If a coordinate system attribute is specified on cen, the center is interpreted in that coordinate system.

Tetrahedron(v1, v2, v3, v4) Each vi must have type 'Vector'(3, algebraic). These represent the vertices of the tetrahedron. If coordinate system attributes are specified on these Vectors, each vertex is interpreted in its own coordinate system.

Triangle(v1, v2, v3) Each vi must have type 'Vector'(2, algebraic). These represent the vertices of the triangle. If coordinate system attributes are specified, each vertex is interpreted in its own coordinate system.

If the specified region contains an additional argument of type list(name), the integration variables will be substituted with the variables in that list. Note that this substitution will only be visible if the integral is returned unevaluated.

The int(f, dom, inert) command returns the unevaluated (inert) form of the integral of f over dom.

$\mathrm{with}\left(\mathrm{VectorCalculus}\right)\:$

$\mathrm{int}\left(x^2+y^2\,\left[x\,y\right]=\mathrm{Circle}\left(⟨0\,0⟩\,r\right)\right)$

$\frac{\mathrm{\pi }{r}^4}{2}$

$\mathrm{int}\left(\exp \left(-x^2-y^2-z^2\right)\,\left[x\,y\,z\right]=\mathrm{Sphere}\left(⟨0\,0\,0⟩\,r\right)\right)$

$\mathrm{\pi }\left(\sqrt{\mathrm{\pi }}\mathrm{erf}\left(r\right){\ⅇ}^{r^2}-2r\right){\ⅇ}^{-r^2}$

$\mathrm{int}\left(\exp \left(-x^2-y^2-z^2\right)\,\left[x\,y\,z\right]=\mathrm{Sphere}\left(⟨0\,0\,0⟩\,1\,\left[r\,\mathrm{\phi }\,\mathrm{\theta }\right]\right)\,'\mathrm{inert}'\right)$

${\int }_0^1{\int }_0^{\mathrm{\pi }}{\int }_0^{2\mathrm{\pi }}{r}^2{\ⅇ}^{-r^2}\sin \left(\mathrm{\phi }\right)\ⅆ\mathrm{\theta }\ⅆ\mathrm{\phi }\ⅆr$

$\mathrm{int}\left(xy\,\left[x\,y\right]=\mathrm{Triangle}\left(⟨0\,0⟩\,⟨1\,0⟩\,⟨0\,1⟩\right)\right)$

$\frac{1}{24}$

$\mathrm{int}\left(\sin \left(x\right)\cos \left(y\right)\tan \left(z\right)\,\left[x\,y\,z\right]=\mathrm{Parallelepiped}\left(0..\mathrm{\pi }\,0..\frac{\mathrm{\pi }}{2}\,0..\frac{\mathrm{\pi }}{4}\right)\right)$

$\ln \left(2\right)$

$\mathrm{int}\left(\frac{\sin \left(x\right)}{x}\,\left[x\,y\right]=\mathrm{Rectangle}\left(0..\frac{\mathrm{\pi }}{2}\,0..\frac{\mathrm{\pi }}{2}\right)\right)$

$\frac{\mathrm{Si}\left(\frac{\mathrm{\pi }}{2}\right)\mathrm{\pi }}{2}$

$\mathrm{int}\left(1\,\left[x\,y\right]=\mathrm{Region}\left(0..1\,x^2..x\right)\right)$

$\frac{1}{6}$

$\mathrm{int}\left(x+y+z\,\left[x\,y\,z\right]=\mathrm{Tetrahedron}\left(⟨0\,0\,0⟩\,⟨1\,0\,0⟩\,⟨0\,1\,0⟩\,⟨0\,0\,1⟩\right)\right)$

$\frac{1}{8}$

$\mathrm{int}\left(x+y+z\,\left[x\,y\,z\right]=\mathrm{Tetrahedron}\left(⟨0\,0\,0⟩\,⟨1\,0\,0⟩\,⟨0\,1\,0⟩\,⟨0\,0\,1⟩\right)\,'\mathrm{inert}'\right)$

${\int }_0^1{\int }_0^{1-t}{\int }_0^{\frac{\left(-1+t\right)\left(x-1+t\right)}{1-t}}\left(x+y+t\right)\ⅆy\ⅆx\ⅆt$

$\mathrm{int}\left(1\,\left[x\,y\right]=\mathrm{Ellipse}\left(\frac{x^2}{4}+\frac{y^2}{9}-1\right)\right)$

$6\mathrm{\pi }$

$\mathrm{int}\left(y^2\,\left[x\,y\right]=\mathrm{Sector}\left(\mathrm{Circle}\left(⟨0\,0⟩\,r\right)\,0\,\mathrm{\pi }\right)\right)$

$\frac{\mathrm{\pi }{r}^4}{8}$

$\mathrm{int}\left(x\,\left[x\,y\right]=\mathrm{Sector}\left(\mathrm{Ellipse}\left(\frac{x^2}{4}+\frac{y^2}{9}-1\right)\,0\,\frac{\mathrm{\pi }}{2}\right)\right)$

$4$

$\mathrm{int}\left(⟨t\,t^2\,t^3⟩\,t\right)$