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;