program nlces
         version 8.0
         if "`1'"=="?" {
                 global S_1 "beta eta rho delta"
                 global beta = 1
                 global  eta = 0.5
                 global  rho = -0.5
                 global  delta = 1
                 exit
         }
         replace `1'=$beta - $eta/$rho*ln($delta*l^(-$rho) + (1-$delta)*k^(-$rho))
end