

   NNuummeerriiccaall IInntteeggrraattiioonn ooff aa FFuunnccttiioonn

        int(f, a, b, eps=1.0e-6, k=5)

   AArrgguummeennttss::

          f: The function to integrate.

          a: Lower bound. (If non-numeric, taken to be -infty.)

          b: Upper bound. (If non-numeric, taken to be infty.)

        eps: Precision.

          k: Number of points for polynomial interpolation.

   DDeessccrriippttiioonn::

        `int' performs numerical integration of a given func-
        tion using the Romberg method.

   EExxaammpplleess::

        f <- function(x) sin(x)+cos(x)-x^2
        int(f,0,2)
        #
        f <- function(x) exp(-(x-2)^2/2)
        int(f,0,"infty")

