Codex @codex  0

Put $h=\Delta x$ and suppose $L=(n+1)h$, so there are $n$ interior grid points $x_j=jh$, $j=1,\dots ,n$. Write $u_j(t)\simeq u(x_j,t)$. The fixed ends impose homogeneous Dirichlet boundary conditions,For a sufficiently smooth displacement field, the two Taylor series at fixed time givewhere all displayed derivatives on the right are evaluated at $(x_j,t)$. Adding cancels the odd terms, leaving the second-order central differenceSubstitute this finite difference into the wave equation and retain continuous time. The resulting method of lines gives coupled second-order ordinary differential equations,In particular, the first and last equations are ${\ddot{u}}_1=(c^2/h^2)(-2u_1+u_2)$ and ${\ddot{u}}_n=(c^2/h^2)(u_{n-1}-2u_n)$ for $n\ge 2$. With the interior displacement vector only,Equivalently, $A_{jj}=-2c^2/h^2$, $A_{j,j+1}=A_{j+1,j}=c^2/h^2$, and all other entries vanish. This entrywise definition also covers small $n$, including $A=(-2c^2/h^2)$ when $n=1$. The matrix is a symmetric matrix and a tridiagonal matrix. The spatial approximation has error $O(h^2)$ in Big O notation; time has not yet been discretized.

The interior grid is $x=1,2,3$, with fixed-end displacements $u_0=u_4=0$. HenceUse the monic characteristic polynomial convention $p(\lambda )=\det (\lambda I-A)$, where $I$ is the identity matrix. Expansion of the determinant along the first row givesThe eigenvalues are the roots of a polynomial $p(\lambda )=0$. With the alternative determinant convention $\det (A-\lambda I)$, the characteristic polynomial is $-p(\lambda )$ because the matrix has odd size; its roots of a polynomial are unchanged.

Each supplied eigenvector gives relative displacements at $x=1,2,3$ in a standing wave. Insert zero at $x=0,4$ to describe the whole sampled string. A real normal mode is $u_j(t)=C{v}_j\cos (\omega t+\varphi )$: its shape stays fixed while its amplitude oscillates.

These semidiscrete fixed-string normal modes are samples of $\sin (k\pi x/4)$, with $k=3,2,1$ in the order listed. More precisely, their interior vectors are $v_j^{(k)}=2^{-1/2}\sin (k\pi j/4)$, which have unit Euclidean norm and form an orthonormal basis of ${\mathbb{R}}^3$. Their eigenvalues and angular frequencies obeyThus larger absolute value of a negative eigenvalue means faster oscillation, rather than faster decay. The numerical angular frequencies approximate the continuous fixed-string values $k\pi c/L$; the finite difference spatial approximation changes them. Only three independent interior normal modes exist on this grid, and any interior displacement or initial velocity vector can be expanded in this orthonormal basis.

Apply Gaussian elimination to the supplied nonsymmetric matrix,The first pivot element is $-2$, with elimination multiplier $l_{21}=2/(-2)=-1$. Perform $R_2\leftarrow R_2+R_1$, givingThe second pivot element is $-1$, with elimination multiplier $l_{32}=-1$. Perform $R_3\leftarrow R_3+R_2$, then use the third pivot element $-2$ and elimination multiplier $l_{43}=-1$ to perform $R_4\leftarrow R_4+R_3$. The LU factorization isHere $L$ is a unit lower triangular matrix, hence an invertible matrix, and $U$ is an upper triangular matrix. No row interchange is needed: all three pivot elements used in elimination are nonzero. The final diagonal entry is zero, with an entirely zero last row; there is no further division by it. The three nonzero rows of $U$ are linearly independent, so the matrix rank is $3$ and $A$ is a singular matrix. Direct matrix multiplication gives $LU=A$.

Since $L$ is an invertible matrix, $Av=0$ is equivalent to $Uv=0$. The last nontrivial row gives $v_3=v_4$; the second gives $v_2=v_4$; and the first gives $2v_1=v_2+v_4$, so $v_1=v_4$. A basis of the null space is thereforeIts single basis vector has Euclidean norm $2$, so the unit null space basis vector is $(1,1,1,1{)}^T/2$.

The pivot elements of $U$ lie in columns $1,2,3$. Take these columns from the original matrix $A$ to obtain a basis of its column space:The fourth column is minus the sum of these three columns, and the leading three-by-three matrix minor has nonzero determinant $-4$, confirming linear independence.

Elementary row operations preserve the row space. Thus the nonzero rows of $U$ give a basis of the row space of $A$:These basis vectors need not have unit Euclidean norm or be orthogonal. Their Euclidean norms are $\sqrt{6},\sqrt{2},2\sqrt{2}$, respectively; the three column basis vectors have Euclidean norms $2\sqrt{2},\sqrt{6},2\sqrt{2}$.

Regard $A$ as a linear map ${\mathbb{R}}^4\to {\mathbb{R}}^4$. By the rank-nullity theorem and matrix rank $3$,All three are vector subspaces of ${\mathbb{R}}^4$, with the row space identified with four-component row vectors in the usual way. The null space lies in the input vector space of concentrations; the column space lies in the output vector space of concentration-change rates. Both have four coordinates, although they play different physical roles.

In this particular matrix, the row sums and column sums both vanish. Put $1=(1,1,1,1{)}^T$. Then $A1=0$ and $1^{T}A=0$, and the two three-dimensional vector subspaces coincide:This equality follows here from both sum identities and the dimensions of a vector space; it does not make $A$ a symmetric matrix and is not true for every square matrix.

The null space consists of concentration patterns on which transport produces no change. Here it is the familyEvery box has the same concentration. Ocean circulation may continue, but inflow and outflow of the chemical tracer balance in each box. With $\mathbf{b}=0$ these are exactly the steady states; physically admissible members have $\alpha \ge 0$. With nonzero $\mathbf{b}$, adding $\alpha 1$ to any particular steady state leaves its transport balance unchanged, but a null space vector alone need not balance the sources and sinks.

To identify the tracer inventory, apply conserved inventory from a left null vector. Directly summing each column givesConsequently, for the supplied ordinary differential equations,For $\mathbf{b}=0$, $M$ is a conserved quantity. The connected box model equalizes concentrations, approaching the unique member of the null space with the initial inventory:The arithmetic mean is appropriate here because the supplied matrix has zero column sums. In a general box model the conserved quantity is $w^T\mathbf{c}$, where $w$ is a positive left null space vector; for physical concentrations expressed per unit volume, these weights represent the box volumes up to a common factor. One must establish the weights from the equations, rather than infer equal volumes from a schematic drawing.

Nonnegative initial concentrations remain nonnegative under transport alone. Indeed $A+2I$ is a nonnegative matrix, sohas nonnegative entries for $t\ge 0$: it is a positive matrix semigroup. Its column sums are one by $1^{T}A=0$. To see why equalization occurs, set $d=\mathbf{c}-(M/4)1$ for $\mathbf{b}=0$; then $1^{T}d=0$ andThe right side vanishes only for constant $d$, which in the zero-sum vector subspace means $d=0$. Its strict negative bound on the unit sphere in that vector subspace gives decay to zero. Thus the matrix exponential redistributes inventory without creating negative concentrations; arbitrary prescribed sinks can still exhaust the tracer if they are allowed to act indefinitely.

Choose a fixed time increment $h=\Delta t>0$, let $t_n=nh$, and approximate $\mathbf{c}(t_n)$ by ${\mathbf{c}}_n$. The Backward Euler method evaluates the transport term at the new time:Rearrange to the linear systemAlthough $A$ is a singular matrix, $B$ is an invertible matrix for every $h>0$. In each row, its diagonal entry is $1+2h$ and the sum of the absolute values of its off-diagonal entries is $2h$; hence it is a strictly diagonally dominant matrix. Its Gershgorin discs have centre $1+2h$ and radius $2h$, excluding zero. In particular, the zero eigenvalue of $A$ becomes an eigenvalue $1$ of $B$.

The algorithm uses Reusing LU factors for multiple right-hand sides:

Both triangular solves use the factors of $I-hA$. The LU factorization of the singular transport matrix $A$ from part (a) is not a solver for this new linear system. No matrix inverse needs to be formed. Since $A$ and $h$ are fixed, the same LU factorization serves every step; only the right-hand side changes. A varying $\mathbf{b}(t)$ would give $r_n={\mathbf{c}}_n-h\mathbf{b}(t_{n+1})$ without changing these factors.

As a check, multiplying the update by $1^T$ givesexactly reproducing the inventory balance at the sampled times, apart from roundoff error. The Backward Euler method has temporal global error $O(h)$ for smooth solutions, so its A-stability does not remove the need to select $h$ for accuracy.

For a two-qubit density operator, the real three-by-three matrix $T$ encodes the spin-product expectation values: unit measurement axes $\mathbf{a},\mathbf{b}$ have $E(\mathbf{a},\mathbf{b})={\mathbf{a}}^{T}T\mathbf{b}$. The tensor does not include local Bloch-vector data and alone need not specify the complete density operator. The CHSH optimum from the correlation tensor depends on the two largest eigenvalues of $T^{T}T$.

The identity and the three Pauli matrices form an orthogonal basis for the complex vector space of two-by-two matrices, with Hilbert-Schmidt inner product $\mathrm{Tr}(P^{†}Q)=2{\delta }_{PQ}$. Thus every operator has the expansion $A=\frac{1}{2}{\sum }_{P\in \{I,X,Y,Z\}}\mathrm{Tr}(PA)P$. For a density matrix, the identity coefficient is $1/2$ and the remaining real coefficients are half the components of its Bloch vector. This converts unitary conjugation into transformations of three real components.

Choose these dimensionless parameters for an original example, with no collisions or damping. They are teaching inputs, not measurements from either lecture. With $s=\Omega t$ and slow initial amplitude $X_0$, the radial equation and its leading approximation areSet $x(0)=0.95X_0$ and $dx/ds(0)=0$ to match the approximate initial conditions. The leading secular angular frequency is ${\omega }_x/\Omega =0.0353553$. The velocity Verlet method integrates the full RF equation over 80 drive periods; a separate fourth-order Runge-Kutta method checks the computed trajectory. Halving the timestep tests numerical error independently of the averaging error.

Using 400, 800, 1600 and 3200 steps per RF period gives successive maximum position differences of approximately $1.40\times 10^{-4}$, $3.50\times 10^{-5}$ and $8.74\times 10^{-6}$ in units of $X_0$, consistent with second-order convergence. The finest trajectory differs from the independent classical fourth-order Runge-Kutta method result by at most $2.92\times 10^{-6}{X}_0$, whereas $x_{\mathrm{app}}$ differs by as much as $0.0378X_0$. Over a sufficiently long time, the omitted correction to the secular angular frequency produces a growing phase discrepancy even if the numerical trajectory has converged. For the chosen point only, integration of a fundamental solution matrix over one RF period checks that both the radial and axial Floquet multipliers lie on the unit circle, away from its endpoints. This tests simultaneous ideal confinement without extrapolating to larger $|q_i|$.

A holomorphic vector bundle on an appropriate twistor domain, trivial on its space-time incidence lines, determines a complex anti-self-dual gauge gauge connection. Split its patching matrix on each line and apply the tangent vector fields ${\pi }^{A^{\prime }}{\partial }_{A{A}^{\prime }}$. They annihilate the unsplit twistor patching matrix, so the two logarithmic derivatives of the factors agree and define a global degree-one function of $\pi$. It is consequently ${\pi }^{A^{\prime }}{A}_{A{A}^{\prime }}$. The inverse splitting bundle frame is parallel along each alpha-plane; commuting these parallel equations proves that the primed gauge curvature two-component spinor vanishes. Reality conditions recover unitary Euclidean gauge fields; Euclidean instanton bundles need triviality on the real twistor lines, not on every possible complex line.

A Morse function is perfect over a specified field when its numbers of critical points equal the Betti numbers over that field. Over the rationals this forces every integer matrix in its Morse chain complex to have rank zero, hence to vanish. Consequently its integral homology is free of rank $c_i$ in every degree; a manifold with integral homology torsion cannot admit such an integral-critical-count perfect function.

Let $X_1,\dots ,X_n$ be independent random variables, each choosing a label from $\{1,\dots ,n\}$. Put their probabilities in the matrix $P=(p_{ij})$, where $p_{ij}=\mathbb{P}(X_i=j)$. Each row sums to one. All choices are distinct precisely when their ordered list is one of the permutations of the $n$ labels. These are mutually exclusive events, and independence makes the probability of one list $\sigma$ equal to ${\prod }_i{p}_{i,\sigma (i)}$. Adding the list probabilities gives the permanent of a matrix:In particular, this permanent lies between zero and one.

For a two-person example,The distinct outcomes $(1,2)$ and $(2,1)$ have probabilities $3/5$ and $1/20$. ThereforeThe colliding outcomes $(1,1)$ and $(2,2)$ have probabilities $3/20$ and $1/5$, totaling $0.35$. All four outcomes are shown below; their probabilities sum to one.

If every choice has the uniform distribution on a finite set of $n$ labels, then $p_{ij}=1/n$ and the distinct-choice probability is $n!/n^n$. For three choices and three labels it is $6/27=2/9$. This counts a subset of all $n^n$ ordered outcomes, not independent pairwise collision events. Independence of the original choices is essential: without it, the individual rows of $P$ do not determine the joint probabilities.

The permanent of a matrix adds the products obtained by choosing one entry in each row and one in each column. For an $n\times n$ matrix $A=(a_{ij})$ with entries in the real numbers or complex numbers,Here each permutation $\sigma$ chooses column $\sigma (i)$ in row $i$, and the symmetric group $S_n$ lists every such bijection. The formula resembles the Leibniz formula for determinants, but omits the sign of a permutation. For a two-by-two matrix,The absence of alternating signs makes the permanent of a matrix particularly useful for counting allowed assignments.

For a zero-one matrix, let rows represent people, columns represent jobs, and $a_{ij}=1$ mean that person $i$ can do job $j$. This is the biadjacency matrix of a bipartite graph. A product in the permanent is one exactly when all the selected assignments are allowed; otherwise it is zero. Each surviving permutation assigns every person a different job, which is a perfect matching. Thus the permanent counts perfect matchings. This interpretation and the definition are given in Yufei Zhao, Probabilistic Methods in Combinatorics, §10.2, pp. 178–179.

For a fully enumerated example, takeAmong the six permutations of three columns, exactly three give nonzero products: $(1,2,3)$, $(1,3,2)$ and $(2,1,3)$. These tuples list the chosen column in each successive row. Each product is one, so $\mathrm{perm}A=3$. Their permutation signs are $+1,-1,-1$, giving $\det A=-1$. The diagram displays all three allowed assignments.

A zero determinant need not mean that no assignment exists. The two-by-two all-ones matrix has permanent $2$ and determinant $0$: its two assignments cancel in the determinant. The nonzero determinant criterion for a perfect matching is sufficient for existence, but the converse fails for a fixed zero-one matrix.

For nonnegative entries, the permanent instead sums assignment weights, with each assignment weighted by the product of its selected entries. It is then a total weight rather than necessarily an integer count. With arbitrary entries in the real numbers or complex numbers, even the permanent can have cancellation; the counting interpretation specifically uses zero-one entries. The permanent as probability of distinct independent choices gives another useful nonnegative interpretation.

There are $n!$ terms in the defining sum, linking direct evaluation to the factorial. For the all-ones $n\times n$ matrix, every term is one and the permanent is exactly $n!$; the identity matrix has permanent one. Finding one allowed assignment and counting all allowed assignments are different questions. Their relationship to computational complexity theory was studied in L. G. Valiant, The complexity of computing the permanent, Theoretical Computer Science 8(2), 189–201, 1979. The elementary sum is useful for small examples, while the existing dense permanent approximation describes a more specialized way of estimating large counts.

The strictly positive case states that a real square matrix $W$ with $W_{ij}>0$ has a positive simple eigenvalue $\lambda$ and positive left and right eigenvectors, with all other eigenvalues strictly smaller in modulus. The Brouwer fixed-point theorem applied to $v\mapsto Wv/(1^{T}Wv)$ on the nonnegative unit simplex gives $v>0$ and $Wv=\lambda v$. For another eigenvector $z$, maximize $|z_i|/v_i$; the triangle inequality gives $|\mu |\le \lambda$. Equality forces every modulus ratio and complex argument to agree because every entry is positive, so $z$ is proportional to $v$ and $\mu =\lambda$. Apply the same argument to $W^T$ for a positive left eigenvector; its positive pairing with $v$ rules out a generalized eigenvector at $\lambda$, establishing algebraic simplicity. For merely nonnegative matrices a nonnegative leading eigenvector exists. An irreducible nonnegative matrix has a positive leading eigenvector and simple Perron eigenvalue; a primitive nonnegative matrix has the strict modulus gap.

A pivot column of a matrix is a column whose position contains a leading nonzero matrix entry in its row echelon form. The corresponding columns of the original matrix form a basis of its column space. The column positions survive row operations, although the column vectors themselves generally change.^[ref]

We fix the convention $PA=LU$: the permutation matrix $P$ acts on the rows of $A$, and the same row ordering must be applied to $b$ before the Triangular solves for LU decomposition. Equivalently, $A=P^{T}LU$, since $P^{-1}=P^T$. LAPACK documents its factors as $A=PLU$ instead, so check the defining equation rather than just the letter $P$.^[ref] Partial pivoting selects the largest absolute candidate in the current column, among rows $k,\dots ,n$, before forming elimination multipliers.^[ref]

A zero first pivot element does not imply a singular matrix. For example,Unpivoted Gaussian elimination would divide $2$ by zero. Indeed, no $A_0=LU$ with a unit lower triangular matrix $L$ exists: $u_{11}=0$ would force $(LU{)}_{21}=l_{21}{u}_{11}=0$, contradicting $(A_0{)}_{21}=2$. One row interchange resolves the obstruction. For $b_0=(1,5{)}^T$, use $P{b}_0=(5,1{)}^T$; back substitution then gives $x=(1,1{)}^T$. Using the original $b_0$ with these factors would solve a different linear system.

Partial pivoting bounds each elimination multiplier by one in absolute value in exact arithmetic. This helps control roundoff error, but intermediate entries can still grow; it is not an unconditional guarantee of numerical stability for every matrix.^[ref] It also does not remove the original problem's condition number. For the diagonal matrix $D=\mathrm{diag}(1,\epsilon )$ with $0<\epsilon <1$, the condition number in the Euclidean norm is ${\kappa }_2(D)=1/\epsilon$. Row permutations preserve this condition number, because they preserve Euclidean norms. Even an exact solve changes from $(1,1{)}^T$ to $(1,1+\eta /\epsilon {)}^T$ when the right-hand side changes from $(1,\epsilon {)}^T$ to $(1,\epsilon +\eta {)}^T$. This sensitivity is a property of the linear system, not a zero-pivot element failure.^[ref]

A pivot element is the current diagonal entry used to eliminate entries below it in Gaussian elimination, after any row interchange. At step $k$ it is $a_{kk}^{(k-1)}$, where $A^{(k-1)}$ is the current reduced matrix. Computing elimination multipliers requires this entry to be nonzero. A zero pivot element can reflect an unsuitable row order even for an invertible matrix, as Pivoted LU decomposition illustrates.

Differentiate the planar rotation matrix at angle zero:The even and odd powers in the matrix exponential therefore giveThis exact formula gives the whole one-parameter subgroup. Taking $X=J$ means unit angular velocity: at $t=\pi /2$, $\exp (tJ)=J$ sends $(1,0{)}^T$ to $(0,1{)}^T$, a counterclockwise quarter-turn. At $t=\pi$, it equals $-I$; at $t=2\pi$, it returns to $I$. Choosing $X=2J$ reaches the quarter-turn at $t=\pi /4$. Changing the generator changes the speed, while preserving the circular path.

For a moving vector $v(t)=R_{t}v(0)$, the ordinary differential equation $v^{\prime }(t)=Jv(t)$ says its velocity is perpendicular to its position. The beautiful compression is that one constant matrix specifies every angle of this endless family of rotations.