anchor.diff:  !============================================================================
   
   program RBC_F90
  -  
     !----------------------------------------------------------------
     ! 0. variables to be defined
     !----------------------------------------------------------------
  + 
  +  implicit none
     
  -  implicit none
  -
  +  
  +  real time_begin, time_end	
     integer,  parameter :: nGridCapital = 17820
     integer,  parameter :: nGridProductivity = 5
  +  
     real(8),  parameter :: tolerance = 0.0000001
  -  
  +  real time_begin, time_end	
  +
     integer :: nCapital, nCapitalNextPeriod, gridCapitalNextPeriod, nProductivity, nProductivityNextPeriod 
     integer :: iteration
           
  -  real    :: elapsed(2), total
  -
     real(8) :: aalpha, bbeta, capitalSteadyState, outputSteadyState, consumptionSteadyState
     real(8) :: valueHighSoFar, valueProvisional, consumption, capitalChoice	
     real(8) :: maxDifference,diff,diffHighSoFar
  
positive.code:   print *, 'My check:', mPolicyFunction(1000,3)
    print *, ' '
  
    total = etime(elapsed)
  
    print *, 'Elapsed time is ', elapsed(1)
  
  end program RBC_F90
  
negative.code:   real(8), dimension(nGridCapital) :: vGridCapital
    real(8), dimension(nGridCapital,nGridProductivity) :: mOutput, mValueFunction, mValueFunctionNew, mPolicyFunction
    real(8), dimension(nGridCapital,nGridProductivity) :: expectedValueFunction
    
    call CPU_TIME(time_begin)
    
    !----------------------------------------------------------------
    ! 1. Calibration
  