using the GSL QAG routine and a comparison with the known result: erf(b). Several values of b from 0.01 to 1 are explored and the integral is reevaluated for each one. The ratio of I(b)/erf(b) is also calculated and is always consistent with 1.
IMPORTANT!!!: The integrand-setup model 'CreateIntegrand' (or a duplicate) must be run before using GSL integration routines within a larger Modelica model (like this one). (Alternatively, one could edit integrand.c directly and then compile "by hand" in a command window:
gcc -fPIC -c integrand.c integration_mi.c
gcc -shared integration.o integration_mi.o -o libgsl_integration_MI.so [.dll for windows])
For the relevant meanings of the 'rel_err', 'limit', 'key', etc., see the 'qag' function and:
.GNU_ScientificLibrary.Examples.integration.ExampleIntegral.ExampleIntegral
