FUNCTION BESSJ(M,X)
      DATA E/0.0001/,P/3.14159/
C
C     ÏOÄÏPOÃPAMMA BÛ×ÈCËEHÈß ÔÓHKÖÈÈ ÁECCEËß 'JM'
C     ÄËß M=0 È M=1
C     TO×HOCTÜ = E
C     ÄËß X>10  JM  BÛ×ÈCËßETCß ÏO ACCÈMÏTOTÈKE
C
      IF(X-10.) 40,40,41
C
C     BÛ×ÈCËEHÈE ÏO ACCÈMÏTOTÈKE
   41 A0=SQRT(2./(P*X))
      IF(M) 42,42,43
   42 D=-.25
      A=X+D*P
      GOTO 44
   43 D=.75
      A=X-D*P
C
   44 BESSJ=A0*(COS(A)-SIN(A)*D/(2.*X))
      RETURN
C
C     BÛ×ÈCËEHÈE ÏO PßÄÓ
   40 A=1
      AM=0.
      IF(M) 45,45,46
   46 A=X*.5
      AM=1.
   45 AK=1.
      BESSJ=A
      D=-X*X*.25
   50 Q=D/(AK*(AK+AM))
      A=A*Q
      BESSJ=BESSJ+A
      AK=AK+1
      IF(Q-1.) 47,50,50
   47 IF(ABS(A)-E) 48,50,50
   48 RETURN
      END