anchor.diff:          END
   
   
  -        
  +
   C       **********************************
   
           SUBROUTINE CH12N(N,Z,NM,CHF1,CHD1,CHF2,CHD2)
  
positive.code:         END
  
  
          
  C       **********************************
  
          SUBROUTINE CBK(M,N,C,CV,QT,CK,BK)
  
negative.code:                  BJ0=SR*(PU0*DCOS(T0)-QU0*DSIN(T0))
                   BJ1=SR*(PU1*DCOS(T1)-QU1*DSIN(T1))
  C                Forward recurrence for J_|v| (Abm & Stg 9.1.27)
  C                It's OK for the limited range -8.0 ≤ v ≤ 12.5,
  C                since x >= 20 here; but would be unstable for v <~ -20
                   BF0=BJ0
                   BF1=BJ1
  