function surfaceQuadrature input Real x1; input Real x2; input Real y1; input Real y2; input SurfaceIntegrand integrand; output Real integral; end surfaceQuadrature;