Application of symbolic mathematics to compute high precision Newton, Cotes, Gauss-legendre-Lobatto-Bond Quadrature formulas