Introduction - If you have any usage issues, please Google them yourself
QUAD Numerically evaluate integral, adaptive Lobatto quadrature.function FUN from A to B to within an error of 1.e-6 using high order
recursive adaptive quadrature. FUN is a function handle. The function
Y=FUN(X) should accept a vector argument X and return a vector result
Y, the integrand evaluated at each element of X.