FUNCTION EBESI(X)
C
C       BÛ×ÈCËEHÈE ÔÓHKÖÈÈ BESI(X)/EXP(X)
C       AËÃOPÈTM BESI(X):
C         X < 1.   - PßÄ TEÉËOPA
C  1.   < X < 3.75 - ÏOËÈHOM 1
C  3.75 < X        - ÏOËÈHOM 2
C
      REAL C1(6),C2(9)
      DATA C1/3.5156229,3.0899424,1.2067492,0.2659732
     *,0.036076,0.0045813/
      DATA C2/0.39894228,0.01328592,0.00225219,-0.00157565
     *,0.00916281,-0.02057606,0.0263553,-0.01647633,0.00392377/
C
      T1=1.
      IF (X.GE.3.75) GOTO 30
      S=1.
      IF (X.GE.1) GOTO 20
C
C  PßÄ TEÉËOPA
      T=X*X*.25
      I=0
 11    I=I+1
       II=I*I
       T1=T1*T/II
       S=S+T1
       IF (T1.GE.1E-6) GOTO 11
 12   EBESI=S*EXP(-X)
      RETURN
C
C  ÏOËÈHOM 1
 20   T=X/3.75
      T=T*T
      DO 21 I=1,6
      T1=T1*T
 21   S=S+C1(I)*T1
      GOTO 12
C
C  ÏOËÈHOM 2
 30   T=3.75/X
      S=C2(1)
      DO 31 I=2,9
       T1=T1*T
 31    S=S+C2(I)*T1
      EBESI=S/SQRT(X)
      RETURN
      END