Nonlinear ODE Stemming from the P-Laplace Equation on Schwarzchild Space I am working on a project studying the p-Laplace equation, and I think it will be a good idea to look at the solution to the corresponding ODE on Schwarzchild space (given by the metric $\tilde{g}=(1+\frac{m}{2r})^{4}[dr^2+r^2(d\theta^2+ \sin^2(\theta)d\varphi^2)])$. I am looking for boundary conditions: initially a constant c and eventually 0 at infinity. After I set up the ODE using local coordinates for the gradient and divergent, I am left with the following incredibly messy ODE:
$$f''(r)+f''(r)^{p-1} f'(r)^{2-p} \frac{1}{p-2}+f'(r)\left(\frac{2}{(p-2)(r+\frac{m}{2})}\right)+f'(r)\left(\frac{m}{r^2}(1+\frac{m}{2r})^{-1}\right)=0$$
Should I bother with continuing to look for a closed form solution to this? It is nonlinear in the worst possible way and is not even separable. One thing I have tried is taking as an ansatz $f'(r)\neq 0$, so that we get
$$\frac{y''}{y'}+\left(\frac{y''}{y'}\right)^{p-1} \frac{1}{p-2}+\left(\frac{\frac{2}{p-2}+\frac{m}{r}}{r+\frac{m}{2}} \right)=0$$
Can this only be solved numerically? I am more interested in the properties of the solution than the closed form (what regularity does it have, etc).
 A: Let $m=1$, $f'(r)=e^{G(r)}$, $G'=g$, $p=2+\eta$, and $\eta>0$. Then the first equation is a polynomial in $g$ of degree $\eta+1$.
$$
\frac{2}{r}\frac{2r+\eta}{2r+1}+g(\eta+g^\eta)=0
$$
When $\eta$ is even, the polynomial is of odd degree thus has at least one real root. Consider the large $r$ behaviour
$$
g(\eta+g^\eta)\sim-\frac{2}{r} \qquad ;\qquad r\to\infty
$$
The left hand side is positive if $g \in\mathbb{R}$ and $\eta$ is odd, the right hand side is negative; thus there are only complex solutions when $\eta$ is odd. The consistent dominant balances are
$$
g \sim -\frac{2}{\eta r}  \qquad ; \qquad r \to \infty 
$$
For the real root when $\eta$ is even, and
$$
g \sim (-\eta)^{1/\eta} \qquad ; \qquad r \to \infty \tag{*}
$$
For the complex roots. Since you are interested in the solutions that remain finite for large $r$, we require $\Re(g)<0$ as $r \to \infty$. In this way you can determine from eq. (*) the imaginary part of the allowed complex roots in terms of $\eta$, which will determine the large $r$ oscillatory behavior of $f'$.
You may perform a similar analysis for $r\to 0$.
Finally, you could consider a perturbation series $g=\sum\limits_{n=0}^\infty g_n\eta^n$. Note $(g_0+\eta g_1+\eta^2g_2+\cdots)^\eta=1+\eta \ln(g_0)+\eta^2(g_1/g_0+\ln^2(g_0)/2)+\cdots$. Matching powers of $\eta$ we find for $\eta^0$
$$
g_0=-\frac{2}{r}\frac{2r}{2r+1}
$$
Matching powers of $\eta^1$
$$
g_1=-\frac{2}{r}\frac{1}{2r+1}-\ln(g_0)
$$
And so on. You'd need to carefully keep track of the different roots.
