MINICERN/mathlib/gen/d/minfccod.inc

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:20  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 *
  10 * minfccod.inc
  11 *
  12       PARAMETER (W5 = 2.23606 79774 99790D0)
  13       PARAMETER (HV = (3-W5)/2, HW = (W5-1)/2, R1 = 1, HF = R1/2)
  14
  15       N=-1
  16       IF(A .NE. B) N=NINT(2.08*LOG(ABS((A-B)/EPS)))
  17       C=A
  18       D=B
  19       IF(A .GT. B) THEN
  20        C=B
  21        D=A
  22       ENDIF
  23       LLT=.TRUE.
  24       LGE=.TRUE.
  25
  26     1 H=D-C
  27       IF(N .LT. 0) THEN
  28        X=HF*(C+D)
  29        Y=F(X)
  30        LLM=ABS(X-A) .GT. DELTA .AND. ABS(X-B) .GT. DELTA
  31        RETURN
  32       ENDIF
  33       IF(LLT) THEN
  34        V=C+HV*H
  35        FV=F(V)
  36       ENDIF
  37       IF(LGE) THEN
  38        W=C+HW*H
  39        FW=F(W)
  40       ENDIF
  41       IF(FV .LT. FW) THEN
  42        LLT=.TRUE.
  43        LGE=.FALSE.
  44        D=W
  45        W=V
  46        FW=FV
  47       ELSE
  48        LLT=.FALSE.
  49        LGE=.TRUE.
  50        C=V
  51        V=W
  52        FV=FW
  53       ENDIF
  54       N=N-1
  55       GO TO 1
  56       END