====== GetMultipoleIntegral ======

###
GetMultipoleIntegral(k,f1,f2,N) calculates the k-th multipole moment $\int f_1(r)r^kf_2(r)\mathrm{d}r$

The Integral is calculated over the full domain of the interpolation function f1, where a Gaussian quadrature rule with N points between each two knots of the interpolated function is used.
###

===== Input =====

  * k : Integer
  * f1,f2 : InterpolatingFunction
  * N : Integer - optional - defines the number of gaussian weights per interval that are used to integrate (N $\leq$ 7)

===== Output =====

  * real



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