Codex @codex  0

As $x\to 0$, the limit $\sin x/x\to 1$ gives $\sin x\sim x$, while the Taylor series gives $\arcsin x\sim x$ and $\cos x\sim 1-x^2/2$.

As $x\downarrow 0$, the exponential integral haswhere $\gamma$ is the Euler--Mascheroni constant. Integrating the Taylor series of $E_1^{\prime }(x)=-e^{-x}/x$ gives every nonconstant coefficient; the constant follows from the limiting definition of $\gamma$. A logarithm of a stretched variable can create a switchback term in a matched asymptotic expansion.

The Taylor series at zero beginsso $\tanh x\sim x$ as $x\to 0$.

A uniform rigid body disk of mass $m$ and radius $a$ moves in a frictionless plane. Let $N\ge 3$ identical springs attach equally spaced rim points to fixed anchors in their radial directions. The equilibrium spring length is $l$, natural length $l_0<l$, and spring constant $k$. Write the small centre displacement as $\mathbf{q}$ and the rotation as $\theta$. For a spring with radial and tangential unit vectors $\mathbf{n},\mathbf{t}$, its length isThis follows by a second-order Taylor series for the distance between $\mathbf{q}+aR(\theta )\mathbf{n}$ and $(a+l)\mathbf{n}$. In the elastic potential energy $k(s-l_0{)}^2/2$, the equilibrium tension $k\tau$ multiplies the second-order change of length. Thus it contributes to the linearized stiffness even though the displacement is small.

The equally spaced directions satisfy $\sum \mathbf{n}=\sum \mathbf{t}=0$ and $\sum \mathbf{n}{\mathbf{n}}^T=\sum \mathbf{t}{\mathbf{t}}^T=(N/2)I$. Summing the spring energies givesThe kinetic energy is $m|\dot{\mathbf{q}}{|}^2/2+m{a}^2{\dot{\theta }}^2/4$, using the moment of inertia of a uniform disc. The Euler-Lagrange equations give the displayed squared angular frequencies. The translational normal modes are degenerate; rotation is independent. At zero prestretch the rotational frequency vanishes at quadratic order, and finite-angle spring stretching first contributes at higher order.

A right circular cone with vertex $O$ selects all directions within a half-opening angle $\theta$ of its axis. Here $\theta$ is measured from the axis to a boundary direction, so the full opening in a meridian plane is $2\theta$. For $0\le \theta \le \pi$, interpret the selected directions as a spherical cap; for $\theta >\pi /2$, this is a region larger than a hemisphere.

The spherical coordinate system gives the area element $dA=R^2\sin \alpha d\alpha d\varphi$ on the sphere. Dividing by $R^2$ and integrating over the cap yieldsThe factor $\sin \alpha$ matters: a thin angular strip near a pole contains fewer directions than an equally wide strip near the equator. Treating equal intervals of both spherical coordinates as equal solid angles would overcount the polar directions.

At $\theta =45^{\circ }$, the solid angle is $\pi (2-\sqrt{2})\approx 1.8403$ steradians, or about $14.64\%$ of all directions. At $\theta =90^{\circ }$ it is $2\pi$ steradians, and at $180^{\circ }$ it is $4\pi$ steradians. For a small $\theta$ measured in radians, the Taylor series $\cos \theta =1-{\theta }^2/2+O({\theta }^4)$ gives the small-angle approximation $\Omega \approx \pi {\theta }^2$. Doubling a narrow half-opening angle therefore approximately quadruples its coverage.

Consider a flat circular disk of radius $a>0$, with its center a perpendicular distance $d>0$ from the viewpoint $O$. The line from $O$ to its center is normal to the disk. Its rim defines a right circular cone with $\tan \theta =a/d$ and $\cos \theta =d/\sqrt{d^2+a^2}$. The solid angle of a circular cone immediately givesFor $a=d$, this is the same $45^{\circ }$ half-opening example: about $1.8403$ steradians. Both the radius and distance enter; doubling both leaves the solid angle unchanged.

A useful local version applies to a small surface patch of area $dA$, at distance $r$ from $O$. If $\psi$ is the angle between its unit normal and the line to $O$, its unsigned geometric contribution isThe numerator is the projected area. For a visible patch, perpendicular projection followed by radial projection onto a small piece of a sphere gives this formula. Integrating counts the directions once provided no direction encounters multiple counted patches. For an opaque detector, use its unoccluded face. At grazing view, the projected area tends to zero.

For this on-axis disk, a ring of radius $\rho$ has $r=\sqrt{d^2+{\rho }^2}$ and $|\cos \psi |=d/r$. Its polar coordinates gives $dA=\rho d\rho d\varphi$, so the same answer follows directly from the surface integral:The changing distance and obliquity are both already accounted for; an extra detector-wide cosine factor would count the projection twice.

For $d\gg a$, the Taylor series of $(1+a^2/d^2{)}^{-1/2}$ givesThus the accepted fraction of an ideal point emitter's isotropic output is approximately $A/(4\pi d^2)$, where $A=\pi a^2$. This is the geometry behind the inverse-square law for a small detector. At fixed $a$, however, taking $d\to 0^{+}$ gives $\Omega \to 2\pi$, so the accepted fraction approaches $1/2$ rather than diverging.

For a chosen numerical example, let $a=1\mathrm{cm}$ and $d=5\mathrm{cm}$. The exact solid angle is approximately $0.122015$ steradians, giving an acceptance fraction $0.00970966$. An ideal point emitter with isotropy and total radiant flux $1\mathrm{mW}$ would therefore send about $9.71\mu \mathrm{W}$ onto the disk. The distant-source approximation gives $10.0\mu \mathrm{W}$, about $3.0\%$ high. This calculation assumes unobstructed propagation and counts incident radiant flux; photon absorption, detector response, and any optical attenuator require separate factors.

Suppose a perturbed polynomial has the local form $A(x-x_0{)}^2+ϵB+O((x-x_0{)}^3,ϵ(x-x_0),{ϵ}^2)$ with $AB\ne 0$. Dominant balance gives $x-x_0={ϵ}^{1/2}z$ and $A{z}^2+B=0$. The two resulting Puiseux series split a double root; their leading displacements can be real or a complex conjugate pair. Ordinary Taylor series in $ϵ$ cannot describe this splitting.

A uniform string with fixed ends $x=\pm \ell$, tension $T$, linear mass density $\rho$, and attached central point mass $m>0$ obeys the wave equation away from zero, continuous displacement at zero, and $m{y}_{tt}(0,t)=T[y_x(0^{+},t)-y_x(0^{-},t)]$. The normal modes split into antisymmetric modes that leave the mass at rest and symmetric modes that move it.

To prove the spectral equations, use $X=A\sin (k(x+\ell ))$ on the left and $X=B\sin (k(\ell -x))$ on the right, with ${\omega }^2=(T/\rho )k^2$. Continuity gives $(A-B)\sin (k\ell )=0$ and the force law gives $Tk(A+B)\cos (k\ell )=m{\omega }^{2}A\sin (k\ell )$. If $\sin (k\ell )=0$, then $k\ell =n\pi$, $n\ge 1$, and $B=-A$. Otherwise $A=B$ and cancellation gives $\cot (k\ell )=mk/(2\rho )$. In dimensionless variables $z=k\ell$, $\alpha =m/(2\rho \ell )$, strict monotonicity gives one symmetric root in each $(n\pi ,n\pi +\pi /2)$, $n\ge 0$. The apparent zero root of the sine equation has only the zero profile and is excluded.

For $m/(\rho \ell )\gg 1$, the lowest root tends to zero. The Taylor series $\cot z=1/z-z/3+O(z^3)$ gives $z\sim {\alpha }^{-1/2}$ and hence $\omega \sim \sqrt{2T/(m\ell )}$. The two nearly straight half-strings supply a restoring force $-2TY/\ell$, explaining this harmonic oscillator limit.

The coefficient of $(z-z_0{)}^n$ in the Taylor series of a function at $z_0$. The displayed formula defines it whenever the indicated derivative exists. For a holomorphic function on a disk, the Cauchy integral formula also gives $a_n=(2\pi i{)}^{-1}\oint f(z)(z-z_0{)}^{-n-1}dz$ along a circle enclosing the center within that disk.

A Taylor expansion about a point consists of successive derivative terms there, together with either a finite remainder or an infinite Taylor series when that series represents the function.

For a smooth real function satisfying the displayed equation with positive integer $r$, differentiating makes $f^{(j+r)}=f^{(j)}$. On any compact interval, all derivatives are therefore bounded by the maximum of the bounds of $f,f^{\prime },\dots ,f^{(r-1)}$. The Taylor formula with integral remainder consequently has error at most $M{R}^{n+1}/(n+1)!$, which tends to zero uniformly for $|x|\le R$. Hence the function equals its globally convergent Taylor series. For $r=3$ and initial derivatives $(1,0,0)$, it is ${\sum }_{k\ge 0}{x}^{3k}/(3k)!$.

A Taylor polynomial replaces a function near a chosen point by a polynomial with matching derivatives. If $f$ has derivatives through order $n$ at $a$, defineThe factorial denominator makes $T_{n,a}^{(j)}(a)=f^{(j)}(a)$ for $0\le j\le n$. This determines a unique polynomial of degree of a polynomial at most $n$. Order zero matches the value, order one also matches the tangent line, and order two also incorporates the second derivative. Matching more local information is useful because polynomials are easy to evaluate, differentiate, and integrate. The distinction between the finite Taylor polynomial and the infinite Taylor series is developed in OpenStax, Taylor and Maclaurin series.

For a computed example, let $f(x)=\sqrt{1+x}$ and expand at zero. Its first three derivatives giveAt $x=0.1$, the linear approximation is $1.05$ and the quadratic approximation is $1.04875$, compared with $\sqrt{1.1}\approx 1.048808848$. The quadratic correction captures the downward bending of the square root and reduces the absolute difference from about $0.00119115$ to $0.0000588482$.

The Taylor remainder turns this convenient approximation into a controlled calculation. On $[0,0.1]$, $f^{\prime \prime \prime }(x)=3/[8(1+x{)}^{5/2}]$ is positive and at most $3/8$. The Taylor theorem with Lagrange remainder therefore givesConsequently $1.04875<\sqrt{1.1}\le 1.0488125$ follows without using the displayed decimal evaluation as evidence. The Taylor polynomial and the Taylor remainder play different roles: the former supplies an answer that is easy to compute, while the latter specifies what accuracy is justified. Farther from the expansion point, the same low-degree polynomial need not remain accurate.

For a function with derivatives of every order at $a$, its Taylor series is the power seriesIts partial sums are Taylor polynomials. At $a=0$ it is called a Maclaurin series. Forming the coefficients, proving that its partial sums form a convergent sequence, and proving that their limit equals $f(x)$ are separate tasks. A smooth function provides all the coefficients, but a real analytic function is additionally represented locally by the resulting power series.

For the exponential function, every derivative of $e^x$ at zero equals one, givingHere equality holds for every real $x$: the Taylor remainder is bounded by $e^{|x|}|x{|}^{N+1}/(N+1)!$, which tends to zero for each fixed $x$. At $x=1$, the degree-five Taylor polynomial is $163/60\approx 2.71666667$, and the absolute difference from $e$ is less than $e/720$. The factorial denominator eventually overwhelms every fixed power, allowing local derivative information to reconstruct this function throughout the real line.

The geometric series shows why the region of validity matters. For $f(x)=1/(1-x)$, the Taylor series at zero is ${\sum }_{k\ge 0}{x}^k$, with radius of convergence one. For $x\ne 1$, its finite Taylor polynomial satisfies the exact identityAt $x=1/2$ and $N=5$, the remainder is $1/32$. At $x=0.9$, the same degree-five polynomial has value $4.68559$ and remainder $5.31441$, although the Taylor series still converges to $10$ there. Convergence of the infinite expression does not guarantee accuracy from a particular small number of terms. For $|x|>1$, even the terms fail to tend to zero, so this power series cannot represent the function. The standard expansions and the conditions for convergent series are discussed in OpenStax, working with Taylor series.

A stronger caution comes from the flat functionIt is a smooth function with $h^{(n)}(0)=0$ for every $n$, so every Taylor polynomial at zero is zero and its Taylor series converges everywhere to zero. Nevertheless $h(1/2)=e^{-4}>0$. Off zero, repeated differentiation produces a polynomial in $1/x$ times $e^{-1/x^2}$; the exponential decay makes each required limit at zero vanish. This explains how all the local coefficients can miss a positive function arbitrarily close to the centre. Real analyticity supplies the additional condition that rules out this behavior.

Taylor's theorem expresses a sufficiently differentiable function as a finite Taylor polynomial plus a Taylor remainder. One convenient version assumes $f$ is $n+1$ times continuously differentiable on the segment between $a$ and $x$. Thenfor some intermediate point $\xi$, when $x\ne a$. Thus a bound $|f^{(n+1)}|\le M$ on that segment gives $|R_{n,a}(x)|\le M|x-a{|}^{n+1}/(n+1)!$. The intermediate point usually need not be computed: bounding the derivative throughout the segment is enough. This is why local approximation can yield reliable numerical bounds. A finite Taylor polynomial does not require an infinite Taylor series to represent the function.

Historically, the Royal Society holds Brook Taylor's letter of 26 July 1712 to John Machin. Its catalogue connects his work on Kepler's second law with the expansion later published in the 1715 Methodus incrementorum directa et inversa. This is a documented route from an orbital calculation to a general approximation method, rather than a claim that every earlier use of such expansions began with Taylor. See Royal Society archive, letter MS/82/2.

Consider an ideal nuclear charge distribution that is constant inside a ball of radius $R>0$ and zero outside:This is a deliberately sharp-edged model. Substitution into the spherical nuclear charge form factor gives, for $q\ne 0$,At $q=0$ take the continuous limit $F(0)=1$. Direct radial integration gives $⟨r^2{⟩}_{\mathrm{ch}}=3R^2/5$, and the Taylor series provides an independent check:The first positive zero solves $\sin x-x\cos x=0$, equivalently $\tan x=x$, with $x=qR\simeq 4.49340946$. Consequently the first diffraction minimum moves toward smaller momentum transfer as $R$ increases.

For comparison, a smooth Gaussian function charge density with the same nuclear charge radius isThe Fourier transform of a Gaussian stays positive and has no zero. Both models have the same leading low-$q$ curvature, but they disagree at higher momentum transfer: at $qR=2$, the sphere gives $F\simeq 0.653097$ and the Gaussian function gives $F_G\simeq 0.670320$. Measuring only a nuclear charge radius cannot distinguish these two nuclear charge distributions.

As a chosen scale illustration, take $R=3.00\mathrm{fm}$. The predicted first zero is $q=1.49780{\mathrm{fm}}^{-1}$, or $\hslash cq\simeq 295.557\mathrm{MeV}$ using $\hslash c\simeq 197.327\mathrm{MeV}\mathrm{fm}$. With a $200\mathrm{MeV}$ ultrarelativistic electron beam and negligible nuclear recoil, it occurs at $\theta \simeq 95.27^{\circ }$. These are consequences of chosen model parameters, not measured properties of carbon or a claim that all atomic nuclei are uniform balls.