MINICERN/mathlib/gen/c/leg1.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:01:58  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE LEG1
  11 #include "legbl.inc"
  12      1,SIGMA,FAC
  13      2,CLOGOK,CDIGAM
  14       DATA GAMMA/0.5772156649/
  15       A=-VV
  16       B=VV+1.0
  17       C=(1.0,0.0)
  18       CALL HYPGM(A,B,C,Z1,PP,ACC,NCVG)
  19       SIGMA=(0.0,0.0)
  20       FAC=(1.0,0.0)
  21       A=0.5*CLOGOK((ZZ+1.0)/(ZZ-1.0),NFRIG,1)-GAMMA-CDIGAM(B)
  22       QQ=A
  23       DO 17 L=1,50
  24       EL=L
  25       SIGMA=SIGMA+1.0/EL
  26       FAC=-FAC*(VV+EL)*(VV-EL+1.0)*Z1/(EL*EL)
  27       B=(A+SIGMA)*FAC
  28       QQ=QQ+B
  29       IF(CRIT(QQ,B,ACC))80,17,17
  30    17 CONTINUE
  31       NCVG=NCVG+2
  32    80 RETURN
  33       END