Wed May 21 22:33:37 2025 

gaussint_test
  R version 4.1.2 (2021-11-01) 
  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.
Wed May 21 22:33:37 2025