quadratureLobatto(function f(), a, b); quadratureLobatto(function f(), a, b, tolerance=100*Modelica.Constants.eps);
Compute definite integral over function f(u,...) from u=a up to u=b using the adaptive Lobatto rule according to:
- Walter Gander:
- Adaptive Quadrature - Revisited. 1998. ftp://ftp.inf.ethz.ch/pub/publications/tech-reports/3xx/306.ps
See the examples in Modelica.Math.Nonlinear.Examples.
function quadratureLobatto extends Modelica.Icons.Function; input Modelica.Math.Nonlinear.Interfaces.partialScalarFunction f "Integrand function"; input Real a "Lower limit of integration interval"; input Real b "Upper limit of integration interval"; input Real tolerance = 100 * Modelica.Constants.eps "Relative tolerance for integral value"; output Real integral "integral value"; end quadratureLobatto;