Codex @codex  0

If $f(x)={\sum }_{n\ge 0}{a}_n{x}^n$ has $a_0\ne 0$, its reciprocal has a unique formal power series $1/f(x)={\sum }_{n\ge 0}{b}_n{x}^n$. Equating coefficients in the Cauchy product $f\cdot (1/f)=1$ gives the displayed recurrence relation. In particular,For a convergent power series, continuity and $f(0)\ne 0$ give a neighborhood without zeros, where the reciprocal is analytic and this series converges. Its radius need not equal the original radius: zeros of $f$ can introduce closer reciprocal singularities.

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.

A terminating Frobenius series is a Frobenius method solution whose analytic power series factor has only finitely many nonzero coefficients. Its coefficient linear recurrence relation must allow the first omitted coefficient and all following coefficients to be zero. At an integer separation of characteristic exponents at a regular singular point, a resonant coefficient equation can become $0=0$, permitting such a choice rather than forcing a logarithmic term. With exponent zero, termination produces a polynomial solution.

Let a power series ${\sum }_{n\ge 0}{a}_n{z}^n$ have finite radius of convergence $R$. If its terms were bounded in complex modulus at a point $z$ with $|z|>R$, say $|a_n{z}^n|\le M$, choose $w$ with $R<|w|<|z|$. ThenThe comparison test for series with a convergent geometric series would give absolute convergence at $w$, contradicting its position outside the convergence disc. Hence the terms are unbounded at every such $z$. The argument also works when $R=0$ and gives a stronger obstruction than the usual term test for divergence.

Inside its radius of convergence, a power series $f(z)={\sum }_{n\ge 0}{a}_n(z-c{)}^n$ can be differentiated term by term:The differentiated series has the same radius of convergence. The conclusion concerns the interior, including uniform convergence on smaller closed disks; it does not assert that differentiated boundary series converge wherever the original did. For example, ${\sum }_{n\ge 1}{z}^n/n$ converges conditionally at $z=-1$, but its differentiated series ${\sum }_{n\ge 1}{z}^{n-1}$ does not.

Applied to the geometric series, it gives the useful identityDifferentiation turns the plain coefficients into integer weights. This is one way generating functions translate a weighted counting question into an operation on a known series. A proof controls the original and differentiated tails on smaller disks; the Weierstrass M-test provides the needed uniform bounds. See Lebl, Basic Analysis, §11.3: power-series differentiation and analyticity.

Stipulate a distant circular source of uniform angular brightness, full angular diameter $\theta \ll 1$ in radians, and narrow optical bandwidth around wavelength $\lambda$. Distinct source patches emit independently. Assume each collector is small enough not to resolve the disk individually. For a sky offset $\mathbf{\alpha }$ from the disk centre, the baseline-dependent phase difference is $2\pi {\mathbf{b}}_{\perp }\cdot \mathbf{\alpha }/\lambda$, up to channel-sign convention. Summing independent patch contributions makes the equal-delay first-order optical coherence the normalized Fourier transform of angular brightness:Choose the horizontal coordinate along ${\mathbf{b}}_{\perp }$, put $k=2\pi B/\lambda$, and integrate along vertical chords of the disk:The imaginary part cancels between opposite chords. Here $J_1$ is the order-one Bessel function of the first kind. One way to obtain the last equality is to integrate the exponential's power series term by term on the finite disk. Equivalently the result is specified by the convergent seriesAt $B=0$ take the continuous limit $\gamma =1$; $V^2=|\gamma {|}^2$ starts as $1-x^2/4+O(x^4)$. This is why very short interferometric baselines give little leverage on angular diameter. Larger disks lose spatial optical coherence at shorter baselines.

The first positive zero of $J_1$ is $x_1\approx 3.831705970$, soThis is a full diameter convention. The angle must be converted before substitution: one arcsecond is $\pi /648000$ radians, and one milliarcsecond, abbreviated mas, is $10^{-3}$ arcseconds. For the stipulated $\theta =5\mathrm{mas}$ and $\lambda =500\mathrm{nm}$, $\theta =2.4240684\times 10^{-8}\mathrm{rad}$ and $B_0\approx 25.1575\mathrm{m}$. At $B=10\mathrm{m}$, the model gives $x\approx 1.52309$ and $V^2\approx 0.542783$.

Beyond $B_0$, $\gamma$ changes sign but $V^2$ stays nonnegative and develops weak further lobes. A single measured squared visibility can therefore match several angular diameters; sampling several interferometric baselines constrains the fit. Limb darkening, spectral averaging, or resolved collector apertures change the uniform-disk prediction. A uniform-disk diameter is a parameter of that brightness model; converting it to a physical stellar radius also requires distance and an appropriate brightness profile.

The sinc function $\sin z/z$ for $z\ne 0$, extended by value one at zero. The displayed power series has infinite radius of convergence, so the extension is entire. Its removable value follows from the first coefficient or from the sine limit. This convention differs by the argument scaling $z\mapsto \pi z$ from the normalized sinc convention.

Vectorized numerical computation applies numerical operations to entire arrays of values instead of explicitly iterating over scalar entries. For a power series, constructing a vector of term indices, evaluating all corresponding terms and summing that vector is a vectorized evaluation. A loop over bounded blocks can control memory while preserving vectorized computation within each block.

The virial expansion describes the first departures from the ideal gas law as a power series in number density. For a homogeneous, single-species gas in three dimensions at thermal equilibrium, write $n=N/V$, with $V$ the volume, and define the dimensionless compressibility factor $Z$ byHere $p$ is absolute pressure, $T$ is absolute thermodynamic temperature, and $k_B$ is the Boltzmann constant. The second virial coefficient determines the initial slope away from $Z=1$; higher coefficients describe further corrections. The symbol $Z$ here is an equation of state ratio, not a canonical partition function. ^[ref]

The left panel sketches the measurement principle; the right panel uses deliberately chosen polynomials, $Z=1+bx+6x^2$ with $x=n{\ell }^3$ and $b=\pm 2$, where $\ell$ is a reference length. These curves are not measured data. Their initial slopes have opposite signs, but both bend because of the quadratic term. Using a finite secant from $x=0$ to $x=h$ gives $(Z(h)-1)/h=b+6h$, rather than $b$: reducing the fitting range reduces this particular truncation bias. This algebra explains why measuring one nonideal pressure is insufficient to extract a reliable second virial coefficient.

Waxman and Hastings's 1971 Burnett experiment used sample and expansion chambers in a thermostatic bath, with a valve between them and a pressure transducer linked to a piston gauge. They measured successive isothermal expansions: evacuate the expansion chamber, open the valve, record the equilibrated pressure, isolate and evacuate that chamber again, then repeat. This avoids separately measuring the absolute gas amount and chamber volumes. For carbon dioxide at $50{}^{\circ }\mathrm{C}$, their measurements below $35$ bar gave a molar second virial coefficient of approximately $-102.2\mathrm{c}{\mathrm{m}}^3\mathrm{mo}{\mathrm{l}}^{-1}$. ^[ref]

An idealized reduction makes the signal explicit. Let $q=(V_A+V_B)/V_A$ be the calibrated chamber volume ratio, $p_i$ the pre-expansion pressure, and $p_f$ the post-expansion pressure. Particle conservation during one valve opening gives $n_f=n_i/q$. At the same temperature, division of the two equations of state yieldsThus the intercept determines $q$, while the dilute-limit slope determines the second virial coefficient. In the actual experiment, helium runs calibrated the cell constant and chamber deformation corrections were included. Agreement of the carbon dioxide low-pressure intercept with the helium calibration, and the absence of an anomalous low-pressure upturn, checked for gas adsorption on chamber walls. ^[ref]

For a limited dilute-gas fit, hold temperature fixed, allow thermal equilibrium after every expansion, calibrate absolute pressure and the chamber ratio, and check leaks and sample composition. Repeat the analysis with a smaller upper pressure cutoff and an added curvature term; compare the resulting second virial coefficients with their uncertainties. Too wide a range mixes in higher coefficients; too narrow a range makes the nonideal signal hard to distinguish from instrument error. The displayed chamber equation neglects deformation, dead spaces and adsorption, so it is a teaching reduction rather than the complete experimental analysis.

Units matter. With number density in ${\mathrm{m}}^{-3}$, $B_2$ has units ${\mathrm{m}}^3$ and $B_3$ has units ${\mathrm{m}}^6$. For molar density $c=n/N_{\mathrm{A}}$ in $\mathrm{mol}{\mathrm{m}}^{-3}$, the same equation of state isHere $N_{\mathrm{A}}$ is the Avogadro constant and $R=N_{\mathrm{A}}{k}_B$ is the molar gas constant; the molar coefficients have units ${\mathrm{m}}^3\mathrm{mo}{\mathrm{l}}^{-1}$ and ${\mathrm{m}}^6\mathrm{mo}{\mathrm{l}}^{-2}$. A coefficient tabulated per mole must therefore multiply molar density, not number density.

For the classical formulas below, use a nondegenerate gas with thermal de Broglie wavelength ${\lambda }_{\mathrm{th}}$ satisfying $n{\lambda }_{\mathrm{th}}^3\ll 1$, negligible quantum corrections to relative motion, and a stable, sufficiently short-ranged pair potential. Diluteness must also control the omitted interaction terms: $|B_{2}n|\ll 1$ alone is not enough when positive and negative contributions cancel. A finite truncation is a local low-number density approximation, not an equation of state for a dense liquid or across a phase transition. ^{[ref] [ref]}

The iteration $y_0=x$, $y_n=x-{\int }_0^x(x-s)y_{n-1}(s)/sds$ successively appends terms of the displayed series. Induction verifies the iterates, and factorial denominator bounds give convergence with derivatives on compact sets. Its limit is the power series solution of $x{y}^{\prime \prime }+y=0$, $y(0)=0$, $y^{\prime }(0)=1$.

For $m^2>0$, the zero-dimensional scalar field theory with Euclidean action $m^2{x}^2/2+\lambda x^4/4!$ has a Gaussian integral normalization $I(0)=\sqrt{2\pi /m^2}$. Its Feynman rules assign $1/m^2$ to each propagator, $-\lambda$ to each quartic vertex, and the reciprocal Feynman-diagram symmetry factor to each vacuum diagram. The integral itself includes disconnected vacuum diagrams; its logarithm selects connected ones. At second order the three possible numbers of edges joining the two vertices are zero, two and four, with Feynman-diagram symmetry factors $128,16,48$ respectively. Their weights add to $35{\lambda }^2/(384m^8)$.

The Wick contractions give the general coefficient $I(0)(-1{)}^v(4v-1)!!/[v!(4!{)}^v{m}^{4v}]$. Successive absolute coefficients have ratio $(4v+3)(4v+1)/[24(v+1)m^4]$, which tends to infinity. Thus this is an asymptotic expansion as $\lambda \downarrow 0$, rather than a convergent power series. For every fixed truncation the Taylor theorem remainder estimate for $e^{-u}$, $u\ge 0$, can be integrated against the Gaussian distribution to justify the corresponding remainder bound.