−Table of Contents
GetMultipoleBesselIntegral
GetMultipoleBesselIntegral(q,k,f1,f2,n) calculates ∫f1(r)jk(qr)f2(r)dr, where jk(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
- Example.Quanty
-- some example code
Result
text produced as output