anchor.diff:    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
  
positive.code:   print *, 'My check:', mPolicyFunction(1000,3)
    print *, ' '
  
    total = etime(elapsed)
  
    print *, 'Elapsed time is ', elapsed(1)
  
  end program RBC_F90
  
negative.code: !============================================================================
  
  program RBC_F90
    !----------------------------------------------------------------
    ! 0. variables to be defined
    !----------------------------------------------------------------
   
    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(8) :: aalpha, bbeta, capitalSteadyState, outputSteadyState, consumptionSteadyState
    real(8) :: valueHighSoFar, valueProvisional, consumption, capitalChoice	
    real(8) :: maxDifference,diff,diffHighSoFar
  