integrate(f,a,b) -- integrate f from a to b numerically, using Gaussian quadrature.
i1 : integrate(sin,0,pi) o1 = 2. o1 : RR