function quadrature2 input Real x1; input Real x2; input Integrand integrand; output Real integral; end quadrature2;