To dismantle this structure, we square both sides to obtain:
f2(x)=r→xlim{r2−x22r2f(r)[f(r)−f(x)]−r3erf(r)}
Focusing on the first term inside the limit, we rewrite the denominator as
(r−x)(r+x):
r→xlim(r−x)(r+x)2r2f(r)[f(r)−f(x)]
As
r→x, the expression
r−xf(r)−f(x) converges to the derivative
f′(x). Evaluating the remaining components at
r=x yields:
2x2x2f(x)=xf(x)
Thus, the first term simplifies to
xf(x)f′(x). The second term evaluates directly to
x3exf(x), resulting in the differential equation:
f2(x)=xf(x)f′(x)−x3exf(x)
Dividing by
x2 reveals the homogeneous structure:
(xy)2=xydxdy−xexy
Applying the substitution
y=vx, where
dxdy=v+xdxdv, the equation becomes:
v2=v(v+xdxdv)−xev
Separating the variables, we get
ve−vdv=dx. Integrating both sides:
∫ve−vdv=∫dx
Using integration by parts, the integral of
ve−v is
−e−v(v+1). Therefore:
−e−v(v+1)=x+C
Given
f(1)=1, we have
v=1 at
x=1. Substituting these values:
−e−1(1+1)=1+C⇒C=−1−e2