Solving numerically the equation of motion of D7 brane perturbation

I want to solve this equation

$$\partial_{\rho}^{2}\phi+\frac{3}{\rho}\partial_{\rho}\phi+\left(\frac{M^{2}}{(1+\rho^{2})^{2}}-\frac{l(l+2)}{\rho^{2}}\right)\phi=0$$

numerically.

I know that this equation can be transformed into the hypergeometric equation through the transformation $$\phi(\rho) = \rho^l (1+\rho^2)^{-\alpha} P(\rho)$$ (in which $P$ is some function) whose exact solution is the well known function see here $$_2 F_1(a,b;c;\rho) = \sum_{n=0}^\infty \frac{(a)_n (b)_n}{(c)_n} \frac{\rho^n}{n!}.$$

The crucial characteristic of this function is that if $a$ or $b$ are negative integers, then the series is finite.

However, I'm interested in exploring a numerical solution for this equation and I would like to know how to obtain numerically the finite series solutions.

Any idea?

Thanks.

migrated from physics.stackexchange.comNov 26 '13 at 8:16

This question came from our site for active researchers, academics and students of physics.

• You might want to try Math SE. Take a look at the non-linear example here: en.wikipedia.org/wiki/… – user6972 Nov 25 '13 at 18:00
• You might be able to identify the form/solution with this book books.google.com/… – user6972 Nov 25 '13 at 18:06