MINICERN/mathlib/gen/d/rgs56p.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:14  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE RGS56P(F,A,B,RES,ERR)
  11       DOUBLE PRECISION E5,E6
  12
  13       PARAMETER (R1 = 1, HF = R1/2)
  14       DIMENSION X5(5),W5(5),X6(6),W6(6)
  15
  16       DATA (X5(I),W5(I),I=1,5)
  17      1/4.6910077030668004D-02, 1.1846344252809454D-01,
  18      2 2.3076534494715846D-01, 2.3931433524968324D-01,
  19      3 5.0000000000000000D-01, 2.8444444444444444D-01,
  20      4 7.6923465505284154D-01, 2.3931433524968324D-01,
  21      5 9.5308992296933200D-01, 1.1846344252809454D-01/
  22
  23       DATA (X6(I),W6(I),I=1,6)
  24      1/3.3765242898423989D-02, 8.5662246189585178D-02,
  25      2 1.6939530676686775D-01, 1.8038078652406930D-01,
  26      3 3.8069040695840155D-01, 2.3395696728634552D-01,
  27      4 6.1930959304159845D-01, 2.3395696728634552D-01,
  28      5 8.3060469323313225D-01, 1.8038078652406930D-01,
  29      6 9.6623475710157601D-01, 8.5662246189585178D-02/
  30
  31       RANG=B-A
  32       E5=0
  33       E6=0
  34       DO 1 I = 1,5
  35       E5=E5+W5(I)*F(A+RANG*X5(I))
  36       E6=E6+W6(I)*F(A+RANG*X6(I))
  37     1 CONTINUE
  38       E6=E6+W6(6)*F(A+RANG*X6(6))
  39       RES=HF*(E6+E5)*RANG
  40       ERR=ABS((E6-E5)*RANG)
  41       RETURN
  42       END