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

Print/export