Fri Mar  6 13:20:11 2020 

gaussint_test
  R version 3.4.4 (2018-03-15) 
  gaussint evaluates a Gauss quadrature rule.


  w <- gauss weights

[1] 1 1

  x <- gauss abscissas

[1] -0.5773503  0.5773503

  f <- function ( x ) { x^3 + x + 1.0 }

  gaussint ( f, x, w ) =  2 

  f <- function ( x ) { cos ( x ) }

  gaussint ( f, x, w ) =  1.675824 

gaussint_test
  Normal end of execution.
Fri Mar  6 13:20:11 2020