PROGRAM D6R20 C Driver for routine BESSJ CHARACTER TEXT*18 OPEN(5,FILE='FNCVAL.DAT',STATUS='OLD') 10 READ(5,'(A)') TEXT IF (TEXT.NE.'Bessel Function Jn') GOTO 10 READ(5,*) NVAL WRITE(*,*) TEXT WRITE(*,'(1X,T5,A,T12,A,T20,A,T33,A)') * 'N','X','Actual','BESSJ(N,X)' DO 11 I=1,NVAL READ(5,*) N,X,VALUE WRITE(*,'(1X,I4,F8.2,2E15.6)') N,X,VALUE,BESSJ(N,X) 11 CONTINUE CLOSE(5) END