Optimization with the value function
A method for solving constrained maximization programs is to use the Bellman dynamic programming approach. This method was developed for the optimal control case in order to evaluate public policy, or to solve problems in forward or backward induction in finite time. The model can be written as a Bellman equation of the form:
![Rendered by QuickLaTeX.com \[ V\left( k_{t},\lambda _{t}\right) =\max_{c_{t},h_{t}}\left[ \ln c_{t}+A\ln \left( 1-h_{t}\right) +\beta E_{t}\left\{ V\left( k_{t+1},\lambda _{t+1}\right) |\lambda _{t}\right\} \right] \]](http://vermandel.fr/wp-content/ql-cache/quicklatex.com-3e530706947697b7926873a4df5a4d10_l3.png)
under the constraints of the three equations mentioned above. Note that the Bellman function also contains
, this means that expectations are conditional on the achievement of
, this comes from the stochastic nature of
. The Bellman equation becomes:
![Rendered by QuickLaTeX.com \[ V\left( k_{t},\lambda _{t}\right) =\max_{k_{t+1},h_{t}}\left[ \ln \left( \lambda _{t} k_{t} ^{\theta } h_{t} ^{1-\theta }-k_{t+1}+\left( 1-\delta \right) k_{t}\right) +A\ln \left( 1-h_{t}\right) +\beta E_{t}\left\{ V\left( k_{t+1},\lambda _{t+1}\right) |\lambda _{t}\right\} \right] \]](http://vermandel.fr/wp-content/ql-cache/quicklatex.com-434b91a77605a50462688cbd48eccdd5_l3.png)
with the control variables
et
. The first order conditions are:

et

Applying the envelope theorem of BenVista-Scheinkman gives:
![Rendered by QuickLaTeX.com \[ \frac{\partial V\left( k_{t},\lambda _{t}\right) }{\partial k_{t}}=\frac{% \theta \lambda _{t}\left( k_{t}\right) ^{\theta -1}\left( h_{t}\right) ^{1-\theta }+\left( 1-\delta \right) }{\lambda _{t}\left( k_{t}\right) ^{\theta }\left( h_{t}\right) ^{1-\theta }-k_{t+1}+\left( 1-\delta \right) k_{t}} \]](http://vermandel.fr/wp-content/ql-cache/quicklatex.com-9eb431d30671b97040e6189285d32cc0_l3.png)
This allows to write the first order condition as follows,

The first order conditions are:
![Rendered by QuickLaTeX.com \[ \left\{ \begin{array}{c} \frac{c_{t+1}}{\beta c_{t}}r_{t}+\left( 1-\delta \right) \\ Ac_{t}=\left( 1-h_{t}\right) w_{t}% \end{array}% \right. \]](http://vermandel.fr/wp-content/ql-cache/quicklatex.com-90570187bbfa85f40761877fb777cc78_l3.png)