Update evaluateFarfieldSecondOrder.m

src/evaluateFarfieldSecondOrder.m

-function [F] = evaluateFarfieldSecondOrder(k, sampling, objects, par, R,q)
-%   Calculates a discretized version of the Born farfield operator by using
-%   the composite trapezoidal rule for the outer integral and calculations as suggested
-%   by Vainikko in order to approximate the individual Born farfields.
-%   INPUT:  kappa - wave number
-%           sampling - cell array that contains all information about the
-%                      discretization of the two spheres
-%                      (detectors, illumination directions)
-%           object - String, that determines the shape of the refractive
-%                    index n=1-q
-%           par - Vector, that contains all parameters for calculating q
-%           R,z - Ass.: supp(q) is contained in B_R(z)
-%   OUTPUT: F - Born farfield operator, nxhat*nd array
+function [F] = evaluateFarfieldSecondOrder(k, sampling, objects, par, R, q)
+
+% Calculates a discretized version of the Born farfield operator of second order 
+% by using the composite trapezoidal rule for the outer integral and calculations
+% as suggested by Vainikko in order to approximate the individual Born farfields.
+%
+% INPUT: kappa      Wave number, >0.
+%        sampling   Structure containing information about the discretization.
+%        objects    Structure including two strings that set the shape-types of
+%                   the two scatterers.
+%        par        Vector containing all further information about the scatterers
+%                   shape, 2*6-array.
+%        R          Sizes of the two scatterers, i.e. radii of balls containing them,
+%                   vector of 2.
+%        q          Values of the contrast function within the individual scatterers,
+%                   >-1.
+%   
+% OUTPUT: F     Born farfield operator of second order, nxhat*nd-array.
+%
+% ************************************************************************************
 
 kR = k * R;
 % kz = k * z;