Herein, the prime denotes the derivative with respect to 7. The general solution of (9a) is h(r) = C e’ with a suitable constant C, which can be determined by imposing the ‘boundary condition’ (9b): C = - eP’, so the following iteration formula can be obtained