CTITLEMBESSL -- BESSEL FUNCTION FOR KAISER WINDOW                       00000100
C***********************************************************************
C                 COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991            *
C***********************************************************************
CA    AUTHOR                      UNKNOWN                               00000200
CA    DESIGNER                    UNKNOWN                               00000300
CA    LANGUAGE                    FORTRAN                               00000400
CA    SYSTEM                      IBM AND CRAY                          00000500
CA    WRITTEN                                                           00000600
CA    REVISED                                                           00000700
CA                                                                      00000800
CA                                                                      00000900
CA    CALL MBESSL (X, Z)                                                00001000
CA                                                                      00001100
CA    IN   X     R4   INPUT ARGUMENT                                    00001200
CA    OUT  Z     R4   RETURNED VALUE                                    00001300
CA                                                                      00002200
C                                                                       00061200
C     FUNCTION INO -- BESSEL FUNCTION FOR KAISER WINDOW                 00061300
C                                                                       00061400
C     REAL FUNCTION INO(X)                                              00061500
      SUBROUTINE MBESSL (X, Z)                                          00061510
      Y  = X/2.0                                                        00061600
      T  = 1.0E-08                                                      00061700
      E  = 1.0                                                          00061800
      DE = 1.0                                                          00061900
      DO 10 I = 1, 25                                                   00062000
         XI  = I                                                        00062100
         DE  = DE*Y/XI                                                  00062200
         SDE = DE*DE                                                    00062300
         E   = E + SDE                                                  00062400
         IF( E*T .GT. SDE ) GOTO                                     20 00062500
   10    CONTINUE                                                       00062600
C  20 INO = E                                                           00062700
   20 Z   = E                                                           00062701
      RETURN                                                            00062800
      END                                                               00062900