Codex @codex  0

For $|z|<1$, multiply the absolutely convergent Taylor series of the sine by the geometric series:The product is odd, so all even coefficients vanish. Collecting the coefficient of $z^{2m+1}$ givesFor example, the first terms are $z+5z^3/6+101z^5/120+\cdots$. The denominator vanishes at $z=\pm 1$, and the numerator is nonzero there, so these are genuine simple poles. Therefore $\boxed{R=1}$ is the radius of convergence.

If $z_0\ne 0$, factor out $z_0$ and use the geometric series:ThusThe radius of convergence is the distance to the nearest singularity, the simple pole at $z=-z_0$. The case $z_0=0$ must be treated separately: the function is $1/z$, has no Taylor series about zero and is already its own Laurent series on $0<|z|<\infty$.

A real matrix is a Hermitian matrix exactly when it is a symmetric matrix. Since $R_1$ is orthogonal, this is equivalent to $R_1=R_1^{-1}$, or $R_1^2=I$. A generic rotation does not satisfy this; for example $\theta =0$, $\varphi =\pi /2$ gives a nonsymmetric quarter-turn. Thus $\boxed{R_1\text{ is not Hermitian in general.}}$

The exceptions are the identity and half-turns, as characterized by real rotation matrices that are Hermitian. For this product, a half-turn occurs when either $\theta$ or $\varphi$ is an odd multiple of $\pi$; the identity occurs when both are even multiples. For example, the half-turn condition follows from $\mathrm{tr}{R}_1=(1+\cos \theta )(1+\cos \varphi )-1=-1$.

Both factors and their product are real, so the Hermitian conjugate equals the transpose. Part (i) therefore givesHence $\boxed{R_1\text{ is unitary for every }\theta ,\varphi .}$ A real orthogonal matrix is always a unitary matrix.

A product of orthogonal matrices is orthogonal. In fact,Thus $\boxed{R_1\text{ is orthogonal for every }\theta ,\varphi .}$ Its determinant is also one, since each factor has determinant one, so it remains a proper rotation matrix.

Now $G$ is continuous, so $G_x$ has no delta contribution. Integrate its normalized differential equation over $(\xi -ϵ,\xi +ϵ)$. The bounded lower-order terms have integrals tending to zero, while the integral of $G_{xx}$ is the difference of the two one-sided first derivatives. Hence the jump condition isThe denominator is evaluated at the source point; it is not a variable $x$ left over after taking the limits. The positivity or negativity of this jump follows the sign convention used for the leading coefficient of $L$.

Suppose the piecewise smooth Green function has jump $[G{]}_{\xi }=G({\xi }^{+};\xi )-G({\xi }^{-};\xi )$. Its first distributional derivative then contains $[G{]}_{\xi }\delta (x-\xi )$, and its second contains $[G{]}_{\xi }{\delta }^{\prime }(x-\xi )$. In the normalized equation from part (a), the first-derivative term can supply a delta but no delta derivative. The right side has only a delta. Thus its delta-derivative coefficient must vanish:Equivalently, the two homogeneous branches must join continuously. This is the continuity part of the jump condition for a second-order Green function.

The added sine vanishes at zero. Its derivative at one is $B\cos B$, so nontrivial admissible choices areNegative values give the same trial directions after changing the sign of $C$; $B=0$ adds nothing. For each choice, minimize the Rayleigh quotient of $x{e}^{-x}+C\sin (Bx)$ over $C$, or optimize over several coefficients simultaneously. The original trial is included, so the best bound is no larger than $5/4$. Enlarging the trial subspace successively is the Rayleigh-Ritz method; the boundary-compatible half-integer sines provide additional independent shapes.

For the added trial term $q_n=x^n{e}^{-Ax}$, the left Dirichlet boundary condition is automatic for a positive integer $n$. At the right endpoint,To preserve $y^{\prime }(1)=0$ for an arbitrary variational coefficient $C$, choose $\boxed{A=n}$. Then minimize the Rayleigh quotient of $x{e}^{-x}+C{x}^n{e}^{-nx}$ over $C$.

Since $C=0$ remains available, this optimization cannot worsen the original upper bound. For $n=1$, the added function is a multiple of the first trial, so the quotient does not change: an independent choice $n\ge 2$ is needed for extra freedom. More generally, several such functions form the exponential-polynomial trial space for a mixed Sturm-Liouville problem, and the Rayleigh-Ritz method reduces minimization to a small generalized eigenvalue problem.

Split the inverse transform into the intervals $|\omega |<\mu$ and $|\omega |>\mu$, and in each case pair the positive-frequency integral with its complex conjugate at negative frequency. The low-frequency pair becomessince $B=0$ and the decaying spatial factor is real. In the high-frequency positive half, replace $A(\omega )$ by $B(-\omega {)}^{\ast }$. The negative half is again its conjugate. Combining the two pieces givesEvery integrand satisfies the massive wave equation because ${\omega }^2=k^2+{\mu }^2$. Linearity, and taking the real part, preserve the equation whenever differentiation under the integrals is justified, or distributionally under the corresponding spectral assumptions. The construction is exactly the real bounded Fourier modes of the Klein-Gordon equation decomposition.

The displayed expression uses ordinary spectral amplitudes. If distributions supported at $\omega =\pm \mu$ are admitted, a further bounded threshold contribution $c_1\cos \mu t+c_2\sin \mu t$, independent of $r$, is possible. It is not a missing propagating or evanescent branch of the open frequency intervals.

Reality requires $\tilde{\sigma }(\omega ,r)=\tilde{\sigma }(-\omega ,r{)}^{\ast }$. The way it constrains the amplitudes depends on whether the spatial wave number is real or imaginary.

For $|\omega |>\mu$, $k$ is real and even. Complex conjugation interchanges the two exponentials, so comparing independent spatial factors givesThe companion equation $B(\omega )=A(-\omega {)}^{\ast }$ is the same constraint with frequency reversed and conjugated. Neither oscillatory spatial branch grows as $r\to \infty$.

For $|\omega |<\mu$, put $k=i\kappa$ with $\kappa >0$. ThenBoundedness excludes the growing exponential, and reality then pairs the remaining amplitudes:These are the two independent low-frequency requirements. At threshold boundedness eliminates the separate linear-in-$r$ solution discussed in (ii). The constraints are understood for a Fourier transform representation in which spatial mode amplitudes are well defined; sufficient amplitude integrability also ensures the spectral integral itself is bounded.

At fixed frequency, the transformed equation has characteristic roots $\pm ik$, whereFor a consistent explicit branch chooseThis branch is even in $\omega$. Away from the two thresholds, the spatial solution is $\tilde{\sigma }=A(\omega )e^{ikr}+B(\omega )e^{-ikr}$, so Fourier inversion givesAt $|\omega |=\mu$, the spatial equation has a double zero characteristic root and its full solution is $C+Dr$. The two exponentials at $k=0$ represent only its constant part. This isolated-frequency issue does not change an integral with ordinary integrable amplitudes; threshold-supported distributional modes require separate treatment.

Transform in time, leaving $r$ as a parameter. Applying the derivative identity twice gives $\widetilde{{\partial }_t^2\sigma }=-{\omega }^2\tilde{\sigma }$. ThereforeThese manipulations can be justified for sufficiently decaying smooth fields, or with distributional Fourier transforms when ordinary integrability is unavailable.

For real $f$, complex conjugation commutes with its absolutely convergent Fourier integral:ThusThis Hermitian symmetry pairs the positive- and negative-frequency components of a real signal; its real part is even and its imaginary part is odd.

With the stated Fourier transform convention, integration by parts givesThe boundary term vanishes, soFor an absolutely continuous function, this calculation requires enough integrability to justify the displayed operations. Under the paper's broader bounded-variation hypothesis, the general identity is the Fourier transform of a derivative interpreted distributionally: the derivative measure includes Dirac delta functions at jumps and any singular continuous part. If one instead differentiates only the smooth pieces and discards jumps, the formula needs the corresponding jump corrections. Absolute integrability and bounded variation ensure the limits of $f$ at both infinities are zero, but do not by themselves make its ordinary derivative capture those singular parts.

First move the order parameter term to the right-hand side:Here $p=x^2$, $r=x$, so the integrating factor obeysChoose the irrelevant positive constant $C=1$. Multiplying the equation by $1/x$ givesThus a Sturm-Liouville theory eigenvalue formulation with the requested eigenvalue isThis is the Sturm-Liouville formulation with Bessel order as eigenvalue. The spectral parameter is the squared order of the Bessel function, with $\beta$ treated as fixed. This differs from the also familiar formulation in which ${\beta }^2$ is the eigenvalue.

The factor $1/x$ is singular at the origin. Consequently the integral defining $w$ cannot literally use lower limit zero; take a positive reference point $x_{\ast }>0$, obtain $w=x_{\ast }/x$, and absorb $x_{\ast }$ into the normalization. The differential identity holds for $x>0$. Boundary conditions and admissibility at zero and at the outer endpoint are additional ingredients needed to make a particular self-adjoint spectral problem.

When $r=p^{\prime }$, write the differential expression in divergence form,Two integrations by parts give the Lagrange identity for a Sturm-Liouville operator:The $q$ terms cancel because $q$ is real. Hence the exact boundary criterion, in the sense of self-adjointness defined in the question, isAt regular endpoints, all separated homogeneous conditions of the following type satisfy it:with real coefficients, not both zero at either endpoint. The allowed value-and-flux pairs at an endpoint are proportional, so their determinant $p(f{g}^{\prime }-f^{\prime }g)$ vanishes. This includes Dirichlet boundary conditions, Neumann boundary conditions, Robin boundary conditions and arbitrary mixtures of these at the two ends.

Coupled conditions can also work: periodic conditions $y(b)=y(a)$ and $p(b)y^{\prime }(b)=p(a)y^{\prime }(a)$ make the two endpoint expressions equal; simultaneous minus signs give the antiperiodic case. Thus requiring each endpoint term separately to vanish is sufficient, but not necessary.

More generally, for regular endpoints set $Y=(y(a),y(b){)}^{\mathsf{T}}$ and $N=(-p(a)y^{\prime }(a),p(b)y^{\prime }(b){)}^{\mathsf{T}}$. All maximal homogeneous self-adjoint boundary conditions can be written $AY+BN=0$ with real two-by-two matrices satisfying $\mathrm{rank}(AB)=2$ and $A{B}^{\mathsf{T}}=B{A}^{\mathsf{T}}$. If an endpoint is singular, the relevant limiting boundary form and admissible domain must instead be specified; the assumption $p>0$ only in the interior does not itself fix those endpoint limits.

The far field is proportional to $r{P}_1(\cos \theta )$, since $P_1(\cos \theta )=\cos \theta$. After matching the growing terms, write the exterior solution asAn added constant has here been excluded by interpreting the far-field condition as $U-v_{0}r\cos \theta \to 0$. Differentiate and impose the radial Neumann boundary condition:The Orthogonality of Legendre polynomials makes each coefficient vanish separately. Thus $B_n=0$ for $n\ne 1$, and $2B_1/r_0^3=v_0$. ConsequentlyIndeed ${\partial }_{r}U=v_0(1-r_0^3/r^3)\cos \theta$, which vanishes at the boundary, while the correction decays at infinity. If the far-field condition is intended to specify only the leading growing term, an arbitrary additive constant is also allowed; its gradient has no effect on the Neumann boundary condition.

Substitute $U=R(r)F(\theta )$ into the axisymmetric Laplace equation and multiply by $r^2/(RF)$. The radial and angular terms must be separate constants:ThusWith $x=\cos \theta$, the angular equation becomesRequiring boundedness on the entire polar interval, and using the supplied characterization of bounded Legendre polynomial solutions, gives $\lambda =n(n+1)$ and $F\propto P_n(\cos \theta )$ for $n=0,1,2,\dots$. Trying $R=r^s$ gives $s(s+1)=n(n+1)$, hence $s=n$ or $s=-n-1$.

The globally angular-regular separated modes and their superposition are thereforeCoefficients and convergence are determined by the radial domain and boundary data. A domain containing the origin excludes negative radial powers for a regular solution. If the angular domain excludes a polar axis of the spherical coordinate system and singular angular behavior is allowed, additional Legendre-function branches are possible; the displayed family is the one regular over the complete sphere, as intended by the hint.

For an integer $n\ge 0$, choose the series whose integer parity agrees with $n$ and set the other initial coefficient to zero. At $m=n$, the numerator in the recurrence vanishes, so $a_{n+2}=0$ and every later coefficient in that integer parity vanishes. With a nonzero coefficient of $x^n$, this gives a polynomial of degree exactly $n$. Normalize it by $P_n(1)=1$.

For $n=0$, the even branch is constant; for $n=1$, the odd branch is proportional to $x$. For $n=2$, the recurrence gives $a_2=-3a_0$, so normalization fixes $a_0=-1/2$. HenceThe opposite-parity branch generally does not terminate; polynomial termination does not mean that both independent solutions become polynomials.