Function for accessing the external GSL CQUAD doubly-adaptive
integration routine, which can handle most types of singularities
in the region of integration.
Keep in mind that the 'integrand' must be "created"
beforehand: see the 'CreateIntegrand' model and the
'integrand_setup' function.
function cquad
extends Modelica.Icons.Function;
input Integer n_par;
input Real par[n_par], a, b, rel_err;
input Integer n_int;
output Real y, err;
end cquad;