MINICERN/mathlib/gen/c/cbsja64.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:09  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10 #if defined(CERNLIB_DOUBLE)
  11       SUBROUTINE WBSJA(Z,A,NMAX,ND,CB)
  12 #endif
  13 #if !defined(CERNLIB_DOUBLE)
  14       SUBROUTINE CBSJA(Z,A,NMAX,ND,CB)
  15 #endif
  16 #include "gen/imp64.inc"
  17 #include "gen/defc64.inc"
  18      +    Z,ZA,CB,I,CBA,RR,F,C,S,R,GCONJG,ZZ
  19       REAL R0,D,T,Q,U,V,TC(11)
  20       CHARACTER*80 ERRTXT
  21       CHARACTER NAME*(*)
  22 #if defined(CERNLIB_DOUBLE)
  23       PARAMETER (NAME = 'CBSJA/WBSJA')
  24 #endif
  25 #if !defined(CERNLIB_DOUBLE)
  26       PARAMETER (NAME = 'CBSJA')
  27 #endif
  28       LOGICAL LER
  29       DIMENSION CB(0:*),CBA(0:100),RR(0:100)
  30
  31 #if defined(CERNLIB_QF2C)
  32 #include "defdr.inc"
  33 #endif
  34
  35 #if defined(CERNLIB_DOUBLE)
  36       PARAMETER (I = (0D0,1D0))
  37 #endif
  38 #if !defined(CERNLIB_DOUBLE)
  39       PARAMETER (I = (0E0,1E0))
  40 #endif
  41       PARAMETER (Z0 = 0, Z1 = 1, HF = Z1/2, Z10 = 10)
  42
  43       DATA TC / 5.7941 E-5,-1.76148E-3, 2.08645E-2,-1.29013E-1,
  44      1          8.5777 E-1, 1.0125 E+0, 7.75   E-1, 2.3026 E+0,
  45      2          1.3863 E+0, 7.3576 E-1, 1.3591 E+0/
  46
  47 #if defined(CERNLIB_DOUBLE)
  48       GCONJG(ZZ)=DCONJG(ZZ)
  49 #endif
  50 #if !defined(CERNLIB_DOUBLE)
  51       GCONJG(ZZ)=CONJG(ZZ)
  52 #endif
  53 #include "c344cod1.inc"
  54 #if defined(CERNLIB_DOUBLE)
  55       F=EXP(YA+I*(A*ATAN2(YA,X)-X))*(HF*AZ)**A/DGAMMA(1+A)
  56 #endif
  57 #if !defined(CERNLIB_DOUBLE)
  58       F=EXP(YA+I*(A*ATAN2(YA,X)-X))*(HF*AZ)**A/GAMMA(1+A)
  59 #endif
  60 #include "c344cod2.inc"