Change of variables formula (${\int }_{T(S)}f(x)dx={\int }_{S}f(T(u))|\det DT(u)|du$)

The change of variables formula expresses an integral in coordinates that make its integrand or its region simpler. Three ingredients change together: the region, the function being integrated, and the element of area or volume. The last change is measured by the absolute Jacobian determinant; forgetting it amounts to counting transformed cells as though their sizes had stayed the same.

In one dimension, substitution is the integrated chain rule. If $g$ is continuously differentiable on $[a,b]$ and $f$ is continuous on its image, thenThese are directed definite integrals: if $g$ decreases, the transformed endpoints reverse. This identity does not require $g$ to be injective, because an antiderivative composed with $g$ already accounts for any backtracking.^[ref]

For a concrete computation, set $u=1+x^2$, so $du=2xdx$. The endpoints $x=0,1$ become $u=1,2$:The numerator is precisely the derivative of the expression in the denominator. Substitution reveals the familiar natural logarithm hidden inside this rational function. Keeping the original endpoints after changing the integrand would describe a different integral.

For the elementary two-dimensional theorem, let $T:U\to V$ be a $C^1$ diffeomorphism between open sets in the plane, let $S\subset U$ be a bounded region with piecewise smooth boundary, and let $f$ be continuous on $T(S)$, with the displayed integrals finite. ThenHere $DT$ is the Jacobian matrix. The diffeomorphism assumption means a bijection with a continuously differentiable inverse, and ensures a nonzero Jacobian determinant. This excludes overlapping sheets; the more general area formula handles multiplicity.^[ref]

Take $S=[0,1{]}^2$ and the linear mapIts inverse is $(u,v)=((x-y)/2,y)$. The square becomes a parallelogram with vertices $(0,0),(2,0),(3,1),(1,1)$: base two, height one, area two. Consequently ${\iint }_{T(S)}1dxdy={\iint }_{S}2dudv=2$. A nonconstant integrand changes too: ${\iint }_{T(S)}xdxdy={\int }_0^1{\int }_0^1(2u+v)2dudv=3$.

The direction of the map matters when selecting the factor. Here $T$ goes from the new coordinates to the original coordinates, so the factor is two. Its inverse has Jacobian determinant $1/2$, appropriate for converting area in the opposite direction. Writing down the map before computing its Jacobian matrix prevents a reciprocal-factor error.

Now reflect the image vertically: $\tilde{T}(u,v)=(2u+v,-v)$ has Jacobian determinant $-2$, but still produces area two. The sign records reversal of the orientation of a vector space; ordinary area uses its absolute value. For oriented differential forms, the pullback of a differential form instead retains the signed determinant.^[ref]

The geometric beauty is that a nonlinear differentiable map has a local linear approximation: tiny rectangles behave approximately like these parallelograms. Summing their scaled contributions gives the change of variables formula. In polar coordinates this factor is $r$, explaining $dxdy=rdrd\theta$; equal angular steps cover more area farther from the origin. Choosing suitable coordinates can simplify both the boundary and the integrand.^[ref]

A change of dependent variable expresses the unknown function of a differential equation in terms of a new unknown, while keeping the independent variable fixed. For example, substituting $y=e^{-ax/2}v$ in $y^{\prime \prime }+a{y}^{\prime }+by=f$ removes the first derivative and gives $v^{\prime \prime }+(b-a^2/4)v=e^{ax/2}f$ when $a,b$ are constant. Transforming the boundary conditions as well is essential, particularly for a Robin boundary condition.

Let $\varphi :I\to J$ be a $C^1$ diffeomorphism between open intervals, which may be unbounded. The change of variables formula holds for every nonnegative Lebesgue measurable function $f:J\to [0,\infty ]$, allowing infinite integrals. Applying it to $|f|$ also proves equivalence of Lebesgue integrability of $f$ and of the transformed integrand; for an integrable signed or complex $f$ the same identity then holds by splitting into positive and negative, or real and imaginary, parts. The absolute value of the derivative is essential for a decreasing substitution such as $\varphi (x)=1/x$. Endpoints can be omitted because they are null sets.

On $x>0$, $0<y<1$, $0<z<x$, this change of variables is a bijection onto $0<\gamma <\beta <\alpha$. Its inverse is $x=\sqrt{\alpha \beta }$, $y=\sqrt{\beta /\alpha }$, $z=\gamma \sqrt{\alpha /\beta }$, and its forward Jacobian determinant is $2x$. Thus $xdxdydz=\frac{1}{2}d\alpha d\beta d\gamma$. Passing to increments $(\alpha -\beta ,\beta -\gamma ,\gamma )$ then maps the ordered cone to the positive octant, turning an exponential with linear exponent into a product of three one-dimensional integrals.

On the positive quadrant, the inverse coordinate change is $x=\sqrt{uv}$, $y=\sqrt{v/u}$. Its Jacobian determinant is $1/(2u)>0$. Regions bounded by positive constant-ratio lines and constant-product hyperbolas become rectangles in these coordinates. An integrand $x/y=u$ then cancels the Jacobian factor, making weighted area integrals particularly simple. The positivity restriction selects a unique square-root branch and fixes orientation.

In the open first quadrant this is one-to-one, with $x=\sqrt{(u+v)/2}$, $y=\sqrt{(u-v)/2}$ and Jacobian determinant $-8xy$. Its image satisfies $u>v$ and $u+v>0$. In particular $1\le v\le u\le U$ is a triangular image for a circular-hyperbolic region. A factor $xy$ in the original integrand cancels the inverse-Jacobian factor $1/(8xy)$; moreover $x^4-y^4=uv$ and $4x^2{y}^2=u^2-v^2$. Boundary degeneracy on $y=0$ is confined to a measure-zero set and must not be confused with an interior failure of invertibility.

For a smooth local diffeomorphism $a:{\mathbb{R}}^d\to {\mathbb{R}}^d$, each inverse branch contributes an inverse Jacobian determinant in the area formula. Their sum is the weighted inverse multiplicity. A lower determinant bound controls each branch separately, not the number of branches. Bounded weighted inverse multiplicity is sufficient for dispersion with a nonlinear velocity map.

The displayed map is a diffeomorphism of the open unit square onto itself. For fixed $u$, the second coordinate decreases strictly from one to zero as $v$ increases from zero to one. Its inverse is $u=1-x$, $v=(1-y)/[1-(1-x)y]$, and its positive Jacobian determinant is $[1-(1-x)y{]}^2/x$. Composing it with $(t,w)\mapsto ((1-t)/(1-wt),1-w)$ gives determinant $[1-(1-x)y][1-(1-x^2)y{]}^2/[x(1-y)]$. The reciprocal determinant therefore integrates to one over the target square by the change of variables formula, despite boundary singularities.

A hyperbolic change of integration variable uses identities such as $1+{\sinh }^{2}u={\cosh }^{2}u$ to simplify square roots. For an open matter-dominated Friedmann solution, $a=b{\sinh }^2(\eta /2)$ makes the cosmic-time integral elementary.

For a nondecreasing real function $z$ on an interval $I$ and nonnegative measurable function $h$,The pushforward of $z^{\prime }(u)du$ is dominated by Lebesgue measure: on each interval, the integral of the derivative is at most the increase of $z$. This handles jumps and singular parts without a smoothness assumption.

For a locally Lipschitz function $T$ between Euclidean spaces of the same dimension, the area formula integrates its absolute Jacobian determinant against the multiplicity of its fibers. In particular, for an injective $C^1$ map $T:X\to {\mathbb{R}}^d$ and nonnegative measurable $a$,This injective version extends the change of variables formula to maps with a possibly singular derivative.

Polar coordinates use $x=r\cos \theta$, $y=r\sin \theta$ and area element $rdrd\theta$.

In planar polar coordinates $(r,\theta )$, the polar radius $r\ge 0$ is the Euclidean distance to the origin. It satisfies $x=r\cos \theta$ and $y=r\sin \theta$. For a differentiable plane curve with parameters $r(t),\theta (t)$, the orthogonality of radial and angular directions gives squared speed ${\dot{r}}^2+r^2{\dot{\theta }}^2$.

Suppose $0\le r_1(\theta )\le r_2(\theta )$ and the region is described without overlap by $\alpha \le \theta \le \beta$ and $r_1(\theta )\le r\le r_2(\theta )$. Integrating the polar area element gives its area:For a region bounded by a single radial graph and two rays, take $r_1=0$. If the two bounding graphs cross, split the angular interval at their intersections so that the larger radius is used consistently.

This change of variables replaces radius by its negative logarithm: $r=\sqrt{x^2+y^2}=e^{-s}$. The angular convention in the title measures $t$ from the positive $y$ axis toward the positive $x$ axis. Its Jacobian determinant is $e^{-2s}>0$, so it supplies local coordinates on the punctured plane. For a single-valued function $u$, its expression $v(s,t)=u(x,y)$ is periodic in $t$ with period $2\pi$.

The multivariable chain rule givesand, for a twice continuously differentiable function,The first-order correction comes from differentiating the coefficients $x,y$ in the angular operator. Consequently the partial differential equation with the right-hand expression equal to zero becomes the heat equation $v_s=v_{tt}$, with the logarithmic radial coordinate acting as the evolution variable.

For polar coordinates $x=r\cos \theta$, $y=r\sin \theta$ with $r\ge 0$, the Jacobian determinant is $r$. The change of variables formula therefore gives the area element $dA=rdrd\theta$. Integrating over the radial interval $r_1(\theta )\le r\le r_2(\theta )$ gives $[r_2(\theta {)}^2-r_1(\theta {)}^2]/2$ per unit angle. Bounds must describe the intended region without counting it more than once.

The cylindrical coordinate system extends planar polar coordinates by an axial coordinate: $(x,y,z)=(R\cos \varphi ,R\sin \varphi ,z)$.

For $x=r\cos \varphi$, $y=r\sin \varphi$ and unchanged $z$, the Jacobian determinant of Cartesian coordinates with respect to cylindrical coordinates is $r$. Hence a volume integral acquires the displayed factor for $r\ge 0$. The inverse determinant is $1/r$ where $r>0$; the axis is a coordinate degeneracy.

The azimuthal angle $\varphi$ in cylindrical coordinates is measured in the horizontal plane from a chosen reference axis. In the usual mathematical convention $x=s\cos \varphi$, $y=s\sin \varphi$, so positive $\varphi$ runs counterclockwise from positive $x$ when viewed from positive $z$. This same angular coordinate occurs in spherical coordinates. It is defined modulo $2\pi$ away from the polar axis and is nonunique on the axis.

The local orthonormal basis associated with cylindrical coordinates: ${\mathbf{e}}_r=(\cos \theta ,\sin \theta ,0)$, ${\mathbf{e}}_{\theta }=(-\sin \theta ,\cos \theta ,0)$, and ${\mathbf{e}}_z=(0,0,1)$. The first two basis vectors depend on the angular coordinate and are undefined on the axis without a choice of angle.

The azimuthal unit vector is tangent to increasing azimuthal angle at fixed cylindrical radius and height. For $r>0$ it equals $(-y/r,x/r,0)$. A counterclockwise circle has tangent displacement $d\mathbf{r}=r\hat{\varphi }d\varphi$. Unlike a fixed Cartesian basis vector, its direction depends on position.

For a differentiable scalar field and cylindrical coordinates with $r>0$, the gradient has the displayed components in the local orthonormal basis. The chain rule applied to $r=\sqrt{x^2+y^2}$ and a local angle $\theta$ gives $r_x=\cos \theta$, $r_y=\sin \theta$, ${\theta }_x=-\sin \theta /r$, and ${\theta }_y=\cos \theta /r$. Substitution in $\nabla \varphi ={\varphi }_x{\mathbf{e}}_x+{\varphi }_y{\mathbf{e}}_y+{\varphi }_z{\mathbf{e}}_z$ proves the formula. The angular factor $1/r$ appears because an angular change $d\theta$ sweeps arc length $rd\theta$.

Differentiate the position vector $\mathbf{r}=r{\mathbf{e}}_r+z{\mathbf{e}}_z$ using ${\dot{\mathbf{e}}}_r=\dot{\theta }{\mathbf{e}}_{\theta }$, ${\dot{\mathbf{e}}}_{\theta }=-\dot{\theta }{\mathbf{e}}_r$, and fixed ${\mathbf{e}}_z$. The first derivative is the velocity $\dot{r}{\mathbf{e}}_r+r\dot{\theta }{\mathbf{e}}_{\theta }+\dot{z}{\mathbf{e}}_z$; differentiating once more gives the displayed acceleration. In particular, $r$ times the azimuthal acceleration is $d(r^2\dot{\theta })/dt$.

If the cylindrical components of a vector field have dependence $e^{im\varphi }$, differentiation of the rotating basis adds the displayed cross product term. Scalar-component differentiation alone misses this contribution. For $\mathbf{B}=(J/2)s\hat{\mathbf{\varphi }}$, one also has $(\mathbf{v}\cdot \nabla )\mathbf{B}=(J/2)\hat{\mathbf{z}}\times \mathbf{v}$. These identities are useful in cylindrical magnetohydrodynamic waves.

The cylindrical radius $R=\sqrt{x^2+y^2}$ is the perpendicular distance from the symmetry axis.

A cardioid is an epicycloid with one cusp; a standard polar equation is $r=a(1+\cos \theta )$.

For $a>0$, the cardioid $r=a(1+\cos \theta )$ traces its boundary once over $0\le \theta \le 2\pi$. The polar area element givesThe polar arc length integrand is $2a|\cos (\theta /2)|$. Splitting at the cusp parameter $\theta =\pi$ gives perimeter $8a$. The powers of $a$ distinguish the scaling of area and length.

For $a>0$ and $0<c<2a$, consider the region $c\le r\le a(1+\cos \varphi )$ in polar coordinates. It exists precisely for $|\varphi |\le \beta$, where $\beta =\arccos (c/a-1)\in (0,\pi )$. The polar area element $rdrd\varphi$ givesExpanding the square and integrating its even terms gives the title formula. This cap lies inside the cardioid and outside the concentric circle, so reversing the radial bounds would select the wrong region or a signed area.

Orthogonal coordinates have pairwise orthogonal coordinate tangent directions. Their metric is diagonal, with scale factors of orthogonal coordinates $h_i=|\partial \mathbf{x}/\partial q_i|$. Spherical coordinates and cylindrical coordinates are examples.

These plane orthogonal coordinates use two families of confocal parabolas. On the region $y\ge 0$, take $u,v\ge 0$ and defineThen $x=(u^2-v^2)/2$ and $y=uv$. The coordinate tangent vectors $(u,v)$ and $(-v,u)$ have zero dot product, and both have length $\sqrt{u^2+v^2}$. Thus the Jacobian determinant is $u^2+v^2$, the area element is $(u^2+v^2)dudv$, and these orthogonal coordinates degenerate only at the origin. Fixed $u=a>0$ gives $x=(a^2-y^2/a^2)/2$; fixed $v=b>0$ gives $x=(y^2/b^2-b^2)/2$. This normalization is a rescaling of the plane section of parabolic cylindrical coordinates.

Without a restricted parameter domain, simultaneous sign reversal $(u,v)\mapsto (-u,-v)$ preserves the image, so these planar parabolic coordinates are globally two-to-one. The inverse function theorem supplies local inverse branches away from the origin. The nonnegative quadrant is one useful branch chart for the upper half-plane.

For planar parabolic coordinates and $g(u,v)=f(x(u,v),y(u,v))$, the multivariable chain rule gives $g_u=u{f}_x+v{f}_y$ and $g_v=-v{f}_x+u{f}_y$. Differentiating once more and adding cancels the first partial derivatives and the mixed partial derivatives:Consequently, on every local coordinate patch away from the origin, $f$ is a harmonic function if and only if $g$ is a harmonic function. The Laplacian changes by a positive scale factor, so solutions of the Laplace equation are preserved. The restriction away from the origin is needed to divide by the Jacobian determinant.

These orthogonal coordinates extend two confocal families of parabolas parallel to an axis. In the displayed scaling their scale factors of orthogonal coordinates are $h_u=h_v=2\sqrt{u^2+v^2}$ and $h_w=1$, so the Laplacian is $\Delta \varphi =({\varphi }_{uu}+{\varphi }_{vv})/[4(u^2+v^2)]+{\varphi }_{ww}$ away from $u=v=0$. With $u\ge 0$, signed $v$ distinguishes the upper and lower half-planes; the two signs at $u=0$ describe the same negative $x$ half-axis. Fixed positive $u^2$ or $v^2$ describes a parabolic cylinder. This differs from the rotational three-dimensional parabolic coordinates, whose coordinate surfaces are elliptic paraboloids.

For a right-handed regular system of orthogonal coordinates, the gradient in orthogonal coordinates gives $\nabla u_i={\mathbf{e}}_i/h_i$. Thus $\varphi {\mathbf{e}}_1=(\varphi h_1)\nabla u_1$, and the curl product rule givesAlso $\nabla u_2\times \nabla u_3={\mathbf{e}}_1/(h_2{h}_3)$ and the divergence of this cross product is zero, by the divergence-of-a-cross-product identity and the vanishing curl of each gradient. Applying the divergence product rule to $\varphi {\mathbf{e}}_1=(\varphi h_2{h}_3)(\nabla u_2\times \nabla u_3)$ givesCyclic permutation and addition give the full component formulas for any vector field. This derivation accounts for the varying coordinate basis without differentiating its unit vectors explicitly. The factors $h_i>0$ are scale factors of orthogonal coordinates; the chart must avoid coordinate singularities.

For a regular chart of orthogonal coordinates $\mathbf{r}(q)$, define $h_i=|{\partial }_{q_i}\mathbf{r}|>0$ and ${\mathbf{e}}_i=h_i^{-1}{\partial }_{q_i}\mathbf{r}$. The tangent unit vectors form an orthonormal basis. Matching $d\Phi =\nabla \Phi \cdot d\mathbf{r}$ with $d\mathbf{r}={\sum }_i{h}_i{\mathbf{e}}_{i}d{q}_i$ proves the formula. Dividing by a scale factor converts a derivative per coordinate increment into a derivative per physical distance in that direction.

For a ring torus of major radius $A>0$, circular torus coordinates use distance $r$ and angle $\alpha$ in the circular cross-section, and azimuth $\beta$ around the central axis. The coordinate tangent vectors are the three partial derivatives of the displayed position. They are mutually orthogonal, with scale factors of orthogonal coordinates $1,r,A+r\sin \alpha$. For $0<r<A$ their ordered Jacobian determinant is $r(A+r\sin \alpha )>0$, so $dV=r(A+r\sin \alpha )drd\alpha d\beta$. At $r=0$ the cross-section angle degenerates. These are tubular coordinates for a circular torus, rather than the confocal coordinate convention also called toroidal coordinates.

For a length scale $c>0$, planar bipolar coordinates may be defined byThe common scale factor is $h=c/(\cosh \tau -\cos \sigma )$, so Laplace's equation becomes ${\Phi }_{\sigma \sigma }+{\Phi }_{\tau \tau }=0$ away from coordinate singularities. For $\tau >0$, constant-$\tau$ curves are circlesTheir radii decrease as $\tau$ increases; for $0<a<b$ the $\tau =b$ circle lies inside the $\tau =a$ circle. The region $a<\tau <b$ is therefore an eccentric circular annulus. This ordering is useful when translating coordinate boundary data into geometric language.

For a harmonic function between $\tau =a$ and $\tau =b$, $0<a<b$, in bipolar coordinates, suppose $\Phi (\sigma ,a)=0$ and $\Phi (\sigma ,b)=F(\sigma )$. Write the periodic boundary values as a Fourier seriesSeparation of variables givesEach radial factor solves $M^{\prime \prime }-n^{2}M=0$ and has endpoint values zero and one; the mean mode solves $M^{\prime \prime }=0$. Sufficiently regular boundary data give the required Fourier series solution. Uniqueness follows from the maximum principle for a harmonic function with zero Dirichlet boundary conditions on both circles. The formula follows the coordinate values; $\tau =b$ is the geometrically smaller circle.

For $a>0$, $u\ge 0$, $0\le v<2\pi$ and real $z$, these orthogonal coordinates have coordinate tangents $r_u=(a\sinh u\cos v,a\cosh u\sin v,0)$, $r_v=(-a\cosh u\sin v,a\sinh u\cos v,0)$ and $r_z=(0,0,1)$. Their scale factors of orthogonal coordinates are $h_u=h_v=a\sqrt{{\sinh }^{2}u+{\sin }^{2}v}$ and $h_z=1$; the ordered Jacobian determinant is positive away from the two focal lines. Positive constant $u$ gives confocal elliptical cylinders, while generic constant $v$ gives a hyperbolic-cylinder half-branch selected by the signs of $x$ and $y$. The limiting $u=0$ surface is the strip $y=0$, $|x|\le a$, with arbitrary $z$. Constant $z$ gives horizontal planes.

In elliptic cylindrical coordinates, write a vector field as $A=A_u{e}_u+A_v{e}_v+A_z{e}_z$ using the local orthonormal unit vectors. The scale factors of orthogonal coordinates are $(h,h,1)$, so the divergence and curl areApplying the divergence to the gradient of a scalar function gives $\Delta \varphi =({\varphi }_{uu}+{\varphi }_{vv})/h^2+{\varphi }_{zz}$. These formulas hold where $h>0$. They use physical components along unit vectors, rather than components along the unnormalized coordinate tangent vectors.

The displayed three-dimensional parabolic coordinates use $u,v\ge 0$ and an azimuthal angle $\varphi$. Their coordinate tangent vectors have pairwise zero dot products, with scale factors of orthogonal coordinates $h_u=h_v=\sqrt{u^2+v^2}$ and $h_{\varphi }=uv$. The Jacobian determinant is $uv(u^2+v^2)$. Constant $u$ or $v$ gives a family of confocal circular elliptic paraboloids; the coordinate system degenerates on the axis where $uv=0$. For example, $0\le u\le v\le 1$ fills the cap $(x^2+y^2-1)/2\le z\le 0$, whose rim is the unit circle and whose vertex is $(0,0,-1/2)$.

This collects the gradient, divergence, curl and Laplacian formulas in cylindrical and spherical coordinates. Coordinate unit vectors vary with position, so taking their derivatives is essential when translating vector calculus operators from Cartesian coordinates. Curl in spherical coordinates is one of these formulas.

Hyperspherical coordinates extend polar coordinates to $n$ dimensions with one radius and $n-1$ angles.

Spherical coordinates $(r,\theta ,\varphi )$ in ${\mathbb{R}}^3$ use

With polar angle $\theta$ measured from positive $z$ and azimuth $\varphi$, the Cartesian formulas $x=r\sin \theta \cos \varphi$, $y=r\sin \theta \sin \varphi$, $z=r\cos \theta$ have Jacobian determinant $r^2\sin \theta$. For $r\ge 0$ and $0\le \theta \le \pi$ this gives the unsigned volume element shown. The inverse determinant is $1/(r^2\sin \theta )$ off the origin and polar axis.

In the usual physics convention for spherical coordinates, the polar angle is the angle between the positive $z$ axis and the position vector. For $r=\sqrt{x^2+y^2+z^2}>0$, it satisfies $\cos \theta =z/r$ with $0\le \theta \le \pi$. The north and south axes have $\theta =0$ and $\theta =\pi$, respectively; the plane $z=0$ has $\theta =\pi /2$. The angle is undefined at the origin. In plane polar coordinates, the same phrase usually denotes the oriented in-plane angle instead, so the reference axis and convention should always be stated. The azimuth supplies the other angular coordinate in three dimensions.

For radial mass density $\rho (r)=A{r}^{-2}(1-e^{-r/r_0})$ within radius $R>0$, integration against the spherical volume element gives $M=4\pi ARf(r_0/R)$. The dimensionless quantity $q=r_0/R>0$ hasThus $f$ is positive, strictly decreasing and convex, with $f(0^{+})=1$ and $f(\infty )=0$. Its asymptotic forms are $f(q)=1-q+q{e}^{-1/q}$ near zero and $f(q)=1/(2q)-1/(6q^2)+O(q^{-3})$ at infinity. The density's central $1/r$ divergence is integrable against the factor $r^{2}dr$.

At fixed angular coordinates, differentiating the Cartesian coordinates with respect to radial distance gives this unit vector. It is the local direction of increasing radius in the spherical coordinate system.

If orthogonal coordinates $q_i$ satisfythen $h_i=|\partial \mathbf{x}/\partial q_i|$ are their scale factors.

These factors determine the diagonal metric and volume element of orthogonal coordinates.

For $\mathbf{A}=A_r{\mathbf{e}}_r+A_{\theta }{\mathbf{e}}_{\theta }+A_{\varphi }{\mathbf{e}}_{\varphi }$, the spherical-coordinate curl follows from the orthogonal-coordinate scale factors $1,r,r\sin \theta$.

Its coordinate expression is part of del in cylindrical and spherical coordinates.

The radius-$R$ ball in ${\mathbb{R}}^n$ has volume ${\pi }^{n/2}{R}^n/\Gamma (n/2+1)$.

Let $V_d$ be the volume of an n-ball of unit radius in ${\mathbb{R}}^d$. A radial shell has volume $d{V}_d{r}^{d-1}dr$, by differentiating $V_d{r}^d$. Integrating shells, or using hyperspherical coordinates, therefore reduces a radial integrand over a Euclidean ball to the displayed one-dimensional integral. The constant $d{V}_d$ is the surface area of an n-sphere of unit radius, with sphere dimension $d-1$.

The boundary of the radius-$R$ ball has area $2{\pi }^{n/2}{R}^{n-1}/\Gamma (n/2)$.