using ApproxFun, SingularIntegralEquations
Following the mastery material, we know that we can re-cast the Cauchy transform in terms of a two-sheeted version of the Cauchy transform. In particular, define C1u(z)=Cu(z)+i√−zC[u/√⋄](z)2C2u(z)=Cu(z)−i√−zC[u/√⋄](z)2
x = Fun(0..100)
u = exp(-x)
C1(u,z) = (cauchy(u,z) + im*sqrt(-z) * cauchy(u/sqrt(x), z))/2
C2(u,z) = (cauchy(u,z) - im*sqrt(-z) * cauchy(u/sqrt(x), z))/2
C1(u,0.1+eps()*im) - C2(u,0.1-eps()*im) ≈ u(0.1)
true
C2(u,0.1+eps()*im) - C1(u,0.1-eps()*im)
2.3703261575747092e-14 + 0.0im
By summing over all inverses of the map z↦z3/2 we can recover the Cauchy transform, with some care on which sheet we map to. In the upper-half plane, z2/3 and e4πi3z2/3 are still on the first sheet as the phase has not crossed 2π, while the third inverse e8πi3z2/3 lives on the second sheet. Therefore, for g(x)=f(x3/2) we have Cf(z)=C1g(z2/3)+C1g(e4πi3z2/3)+C2g(e8πi3z2/3)
g = exp(-x)
C_ex = 0.27264293583876564 - 0.20980563924551332im # "exact" Cf(0.1+0.1im) from Mathematica
z = 0.1+0.1im
C_ex ≈ C1(g, z^(2/3)) + C1(g, exp(4π*im/3)z^(2/3)) + C2(g, exp(8π*im/3)z^(2/3))
true
As we pass to the lower half plane over the negative real axis, the second root passes to the second sheet. However, simealtaneously the roots swap plasses: for x<0 we have for λ1(z):=z2/3, λ2(z):=e4iπ/3z2/3, λ3(z)=e8iπ/3z2/3: λ+1(x)=x2/3+=e4iπ/3x2/3−=λ−2(x)λ+2(x)=e4iπ/3x2/3+=e8iπ/3x2/3−=λ−3(x)λ+3(x)=e8iπ/3x2/3+= x2/3−=λ−1(x)
and we therefore get: Cf(z)=C1g(e4iπ/3z2/3)+C2g(e8πi3z2/3)+C2g(z2/3)
C_ex_m = -0.27264293583876564 - 0.20980563924551332im # "exact" Cf(0.1-0.1im) from Mathematica
z = 0.1-0.1im
C_ex_m ≈ C1(g, exp(4π*im/3)z^(2/3)) + C2(g, exp(8π*im/3)*z^(2/3)) + C2(g, z^(2/3))
true
In other words, in terms of g and h(x)=g(x)/√(x) we have: Cf(z)=12{Cg(λ1(z))+i√−λ1(z)Ch(λ1(z))+Cg(λ2(z))+i√−λ2(z)Ch(λ2(z))+Cg(λ3(z))−i√−λ3(z))Ch(λ3(z)))ℑz>0Cg(λ1(z))−i√−λ1(z)Ch(λ1(z))+Cg(λ2(z))+i√−λ2(z)Ch(λ2(z))+Cg(λ3(z))−i√−λ3(z))Ch(λ3(z)))ℑz<0
We specialize the above for f(x)=e−x2/3 so that g(x)=e−x and h(x)=e−x/√x. From Problem Sheet 3, Problem 4.2 we have for 0<α<1 C[⋄αe−⋄](z)=1Γ(−α)(−z)αe−zΓ(−α,−z)e−iπα−eiπα
For α=0, we use the alternative for x<0 Γ(0,x+)−Γ(0,x−)=−∮e−ζζdζ=−2πi