Muy bien, esta respuesta es un tiro en la oscuridad, pero aquí va.
Primero, transforme el ODE de segundo orden en un sistema de dos ODE. Dejar
φ1φ2=ψ,=ψ˙,
donde los puntos sobre las funciones corresponden a la diferenciación con respecto a la variable independiente (en este caso, ).ξ
Entonces el ODE implícito de segundo orden
ψ¨(ξ)+2ξ−1ψ˙(ξ)ψ(0)ψ˙(0)=e−ψ(ξ)=0=0
puede expresarse como el ODE explícito de primer orden
φ˙1(ξ)φ˙2(ξ)φ1(0)φ2(0)=φ2(ξ)=−2ξ−1φ2(ξ)+e−φ1(ξ)=0=0.
Al principio, parecería que no podemos evaluar el lado derecho de este sistema ODE explícito en , como requiere un integrador numérico. Si existe una solución para este sistema, entonces debe ser diferenciable. Suponiendo que existe una solución, tome el límite del lado derecho como .ξ → 0ξ=0ξ→0
Primero, sabemos que
limξ→0φ2(ξ)=0,
porque asumimos que existe una solución, entonces es diferenciable, lo que significa que debe ser continua. El límite de una función continua en un punto es su valor en ese punto, y sabemos el valor de porque es una condición inicial.φ2φ2(0)
También sabemos que
limξ→0e−φ1(ξ)=1
por razones similares; hemos asumido que es diferenciable, por lo que es continuo, y porque es una condición inicial.φ1φ1(0)=0
Finalmente,
limξ→0−2φ2(ξ)ξ=limξ→0−2φ˙2(ξ),
usando la regla de L'Hôpital en la forma indeterminada .0/0
Para continuar, tenemos que hacer otra suposición: es continua en . Entonces se sigue queφ˙2ξ=0
limξ→0−2φ˙2(ξ)=−2φ2˙(0).
Revisando el ODE de primer orden y evaluando el lado derecho en , podemos ver que tenemos:ξ=0
φ˙1(0)=0φ˙2(0)=−2φ˙2(0)+1,
de donde se deduce que .φ˙2(0)=1/3
Usando este análisis, podría conectar una if
declaración que devuelva estos valores de la función del lado derecho en , lo que debería superar la singularidad. Dicho esto, este análisis requiere un par de suposiciones sobre la continuidad que pueden o no mantenerse, así que tome la solución resultante con un grano de sal.ξ=0