Integrate function y = integrand(x) from x1 to x2
function quadrature input Real x1; input Real x2; input Integrand integrand; output Real integral; end quadrature;