19 May 2024 10:32:34 PM

ppc_gauss_quad_test():
  C version
  Test ppc_gauss_quad():

demo1():
  Request a 10 point quadrature rule.
  Print the points and weights.

      Weight      Point

      0.0667     -0.9739
      0.1495     -0.8651
      0.2191     -0.6794
      0.2693     -0.4334
      0.2955     -0.1489
      0.2955      0.1489
      0.2693      0.4334
      0.2191      0.6794
      0.1495      0.8651
      0.0667      0.9739

demo2():
  Estimate an integral with a sequence of rules.

  N           Estimate       Error

   1        2.22222222        0.22996981
   2        2.13438735        0.14213494
   3        1.98327360        0.00897881
   4        1.97870457        0.01354784
   5        1.99363049        0.00137808
   6        1.99319180        0.00093939
   7        1.99215177        0.00010064
   8        1.99217844        0.00007396
   9        1.99226144        0.00000903
  10        1.99225789        0.00000548

ppc_gauss_quad_test():
  Normal end of execution.
19 May 2024 10:32:34 PM