ADAS Subroutine r8yip

      FUNCTION R8YIP( XI , DELTA )
C
C
C-----------------------------------------------------------------------
C
C  **************** FORTRAN77 REAL*8 FUNCTION: R8YIP  ******************
C
C  PURPOSE: EVALUATES IMPACT PARAMETER CROSS-SECTION SECOND BESSEL   
C           INTEGRAL Y [BURGESS AND SUMMERS: MNRAS (1976) 
C           172,345 - EQN C14]   
C
C  CALLING PROGRAMS: EIQIP, CXEIQP
C
C  INPUT:  (R*8)  XI      = Z[1/KN-1/KN1]/A0 WITH Z TARGET CHARGE, 
C                              KN, KN1 INITIAL AND FINAL ELECTRON WAVE
C                              NUMBERS AND A0 THE BOHR RADIUS
C  INPUT:  (R*8)  DELTA   = RC[KN-KN1] WITH RC THE CLOSEST APPROACH
c
C  OUTPUT: (R*8)  R8YIP   = Y(XI,DELTA) 
C
C ROUTINES: NONE
C
C AUTHOR:  JONATHAN NASH (TESSELLA SUPPORT SERVICES PLC)
C          K1/0/81
C          JET EXT. 5183
C
C DATE:    07/10/93
C
C VERSION 1.1                             DATE: 07/10/93
C MODIFIED: JONATHAN NASH
C           - FIRST VERSION
C VERSION 1.2                             DATE: 17-04-07
C MODIFIED: HUGH SUMMERS
C           - COMPLETED COMMENT BLOCK DESCRIPTION
C
C-----------------------------------------------------------------------
C
C-----------------------------------------------------------------------
      REAL*8              DELTA,       XI