mod_csout Module

This subroutine does an outwards integration of $l_{lmax}$ functions f with an extended 3-point-Simpson:

$$\begin{equation} f_{int}\left(ll,i\right) = \int_{ir}^{ir_{max}}f\left(ll,i'\right)di' \end{equation}$$

The starting value for this integration at $ir_{min}+1$ is determined by a 4 point lagrangian integration, coefficients given by m. abramowitz and i.a. stegun, handbook of mathematical functions, nbs applied mathematics series 55 (1968)

B. Drittler Mar. 1989

Modified for functions with kinks - at each kink the integration is restarted

B. Drittler july 1989