Core idea

For ν>0\nu>0 and smooth FF near Y=0Y=0, the regular solution of

YG+νG=F,G(0)=0YG''+\nu G'=F,\qquad G(0)=0

is

(JνF)(Y)=0Y01sν1F(st)dsdt.(J_\nu F)(Y)=\int_0^Y\int_0^1s^{\nu-1}F(st)\,ds\,dt.

The integral is smooth through zero. For Y>0Y>0, differentiating YνGY^\nu G' verifies the equation. The homogeneous derivative is a constant times YνY^{-\nu}, so C1C^1 regularity and G(0)=0G(0)=0 give uniqueness on the nonnegative side.

Coefficients

For an analytic input F(Y)=α0FαYαF(Y)=\sum_{\alpha\ge0}F_\alpha Y^\alpha, termwise integration gives

(JνF)0=0,(JνF)α+1=Fα(α+1)(α+ν).(J_\nu F)_0=0,\qquad (J_\nu F)_{\alpha+1}=\frac{F_\alpha}{(\alpha+1)(\alpha+\nu)}.

The operator raises the smallest radial degree by one and divides by two degree factors. Its boundedness in any chosen coefficient norm must be checked against that norm's weights.