anchor.diff:      REAL(our_dble), INTENT(OUT)     :: rslt
   
       REAL(our_dble), INTENT(IN)      :: periods_payoffs_systematic(num_periods, max_states_period, 4)
  +    REAL(our_dble), INTENT(IN)      :: periods_draws_prob(num_periods, num_draws_prob, 4)
       REAL(our_dble), INTENT(IN)      :: data_array(num_agents_est * num_periods, 8)
  -    REAL(our_dble), INTENT(IN)      :: periods_draws_prob(:, :, :)
       REAL(our_dble), INTENT(IN)      :: shocks_cholesky(4, 4)
       REAL(our_dble), INTENT(IN)      :: periods_emax(:, :)
   
  
positive.code:     !/* internal objects        */
  
      REAL(our_dble)                  :: crit_val(num_agents_est * num_periods)
      REAL(our_dble)                  :: shocks_cov(4, 4)
  
      INTEGER(our_int)                :: edu_lagged
      INTEGER(our_int)                :: counts(4)
  
negative.code:     INTEGER(our_int)                :: k
      INTEGER(our_int)                :: j
  
      LOGICAL                         :: is_deterministic
      LOGICAL                         :: is_working
  
  