anchor.diff:    ! IMSL library in the original work.
     !*****************************************************************************
   
  +  IMPLICIT NONE
  +
  +  PUBLIC
  +
     CONTAINS
   
   !*******************************************************************************
   !*******************************************************************************
  -SUBROUTINE LFCDS(a, matrix, b , factor, c, d)
  +SUBROUTINE LFCDS(a, matrix, b, factor, c, d)
   ! I WANT THE INTERFACE CLOSER TO THE ORGINAL IMSL IF POSSIBLE
   !This subroutine generates the cholesky factor of the argument 'matrix' and stores
   !the result in the argument 'factor'.
  
positive.code:   ! The input and output connections are opened and the commenting
    ! style is changed to comply with the FORTRAN95 standard. Line
    ! breaks are removed.
    !*****************************************************************************
  
  	!Interface to IMSL replacements
  	USE IMSL_REPLACEMENTS
    
  !********************************************************
  !*  PROGRAM TO CONSTRUCT MONTE CARLO DATA FOR DP MODEL  *
  
negative.code:       INTEGER X1,X2,E,T
  
  
    ! Open files
    open(9,file='in1.txt')
      
    open(10,file='otest.txt')
      
    open(11,file='ftest.txt')
  
  
        READ(9,1500) NPER,NPOP,DRAW,DRAW1,TAU
  