====== GetMultipoleBesselIntegral ======

###
GetMultipoleBesselIntegral(q,k,f1,f2,n) calculates $\int f_1(r)j_k(qr)f_2(r)\mathrm{d}r$, where $j_k(qr)$ is a spherical Bessel function. The functions f1 and f2 are interpolating functions and the domain of the integral is the domain of interpolating function f1. It uses Gaussian quadrature of order n for each interval between the knots of the interpolating function f1.
###

===== Input =====

  * q : real
  * k : integer
  * f1: InterpolatingFunction
  * f2: InterpolatingFunction
  * n : unsigned int - optional, default is 3, maximum is 7

===== Output =====

  * real

===== Example =====

###
description text
###

==== Input ====
<code Quanty Example.Quanty>
-- some example code
</code>

==== Result ====
<file Quanty_Output>
text produced as output
</file>

===== Table of contents =====
{{indexmenu>.#1}}