Table of Contents
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