MINICERN/mathlib/gen/d/simpscod.inc

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:13  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 *
  10 * simpscod.inc
  11 *
  12       DIMENSION F(0:*)
  13
  14       IF(N2 .LE. 0 .OR. 2*(N2/2) .NE. N2) THEN
  15        H=0
  16        WRITE(ERRTXT,101) N2+1
  17        CALL MTLPRT(NAME,'D101.1',ERRTXT)
  18       ELSE
  19 *      S1=0
  20        S1=F(N2-1)
  21        S2=0
  22 *      DO 1 N = 1,N2-1,2
  23 *   1  S1=S1+F(N)
  24 *      DO 2 N = 2,N2-2,2
  25 *   2  S2=S2+F(N)
  26
  27        DO 1 N = 1,N2-3,2
  28          S1=S1+F(N)
  29          S2=S2+F(N+1)
  30     1  CONTINUE
  31 *      S1=S1+F(N2-1)
  32
  33        S1=S1+S1+S2
  34        H=(F(0)+F(N2)+S1+S1)*(B-A)/(3*N2)
  35       ENDIF