Chapter 9 Solow O. Afonso, P. B. Vasconcelos Computational Economics: a concise introduction O. Afonso, P. B. Vasconcelos Computational Economics 1 / 27 Overview 1 Introduction 2 Economic model 3 Computational implementation 4 Numerical results and simulation 5 Highlights 6 Main references O. Afonso, P. B. Vasconcelos Computational Economics 2 / 27 Introduction In 1956, Robert Solow published a seminal paper on economic growth and development. The proposed model, also known has the Solow–Swan, ignores some important aspects of macroeconomics, such as short-run fluctuations in employment, and makes several assumptions to describe the long-run path of the economy. It remains highly influential even today and despite its relative simplicity conveys a number of very useful insights about the dynamics of the growth process. It is also worth teaching from a methodological perspective. Numerical methods for solving differential equations will be introduced and the Euler method implemented to solve initial value problems (Dahlquist and Björck (2008) and Süli and Mayers (2003)). O. Afonso, P. B. Vasconcelos Computational Economics 3 / 27 Economic model This is a model of capital accumulation in a pure production economy. Countries produce and consume only a single, homogeneous good (output, GDP or real income). There are no prices because there is no need for money. Everyone works all the time, saves a fixed portion of income, invests and owns the firm (consumer side is not modelled). There is no government, thus no taxation nor subsidies. It is a closed economy model. This model captures the pure impact that savings (investment) has on the long-run standard of living (per capita income). The model is built around two equations: a production function and a capital accumulation equation. O. Afonso, P. B. Vasconcelos Computational Economics 4 / 27 Economic model Production function The neoclassic production function is assumed to have the Cobb–Douglas form: α 1 α Y (t) = F(K (t); L(t)) = K (t) L(t) − (1) where K (t) is capital input, L(t) is labour input and α, between 0 and 1, is the capital share in production. This production function is neoclassic since it exhibits constant returns to scale, presents positive and diminishing marginal returns to factor accumulation, satisfies the Inada conditions. O. Afonso, P. B. Vasconcelos Computational Economics 5 / 27 Economic model Production function In per capita terms, the production function (1) is: Y (t) K (t)αL(t)1 α = − , y(t) = k(t)α (2) L(t) L(t) Y (t) K (t) where y = L(t) and k(t) = L(t) From (1), taking logs and differentiating with respect to time on both sides, denoting by Y_ (t) and L_ (t) the derivatives of Y and L with respect to t, the results is Y_ (t) K_ (t) L_ (t) = α + (1 − α) : Y (t) K (t) L(t) The growth rate of per capita output is then simply y_ (t) k_ (t) = α : y(t) k(t) Thus, the source of increase in per capita output is the capital deepening. O. Afonso, P. B. Vasconcelos Computational Economics 6 / 27 Economic model Capital accumulation The second key equation describes the path of capital accumulation K_ (t) = Y (t) − C(t) − δK (t) (3) i.e. the addition to the capital stock each period, K_ (t), depends positively on savings and negatively on depreciation, which takes place at rate δ. Since a fraction s of output is saved, Y (t) − C(t) = sY (t), the labour input L_ (t) grows at rate L(t) = n, the path of per capita capital accumulation is k_ (t) = sy(t) − (δ + n)k(t): (4) Thus, the change in per capita capital each period is determined by three terms: per capita investment, sy(t); per capita depreciation, δk(t); and population growth, nk(t). This is a differential equation and, for an initial stock of capital k0, it defines an initial value problem. O. Afonso, P. B. Vasconcelos Computational Economics 7 / 27 Economic model Transitional dynamics: the Solow diagram The Solow diagram has two curves, plotted as functions of k(t), and can be used to understand per capita output evolves over time f(k) (δ+n)k consumption sf(k) investment k∗ Solow diagram O. Afonso, P. B. Vasconcelos Computational Economics 8 / 27 Economic model Transitional dynamics: the Solow diagram The curve y(t) = f (k(t)) = k(t)α is the production function, and the curve sy(t) = sk(t)α depicts the amount of per capita investment, which shifts y(t) down by the factor s. In turn, the line (δ + n)k(t) represents the amount of new per capita investment required to keep the amount of per capita capital constant. When the difference between the curve sy(t) and the line (δ + n)k(t) is positive (negative) the economy is increasing (decreasing) its per capita capital: capital deepening (widening) is occurring. when sy(t) = (n + δ)k(t), so that k_ (t) = 0, the amount of per capita capital remains constant, and such a point is a steady state. To sum up, we have the following. From the diagram it is found the steady-state value of per capita capital. Then, from the production function, the also constant steady-state value of output per worker is obtained. In turn, the steady state per capita consumption is given by the difference between steady state per capita values of output and investment. O. Afonso, P. B. Vasconcelos Computational Economics 9 / 27 Economic model Steady state Since, in steady state, k_ (t) = 0, (4) and (2) can be used to solve the steady-state values of per capita capital and per capita output; thus, α α k_ (t) = sk(t) − (δ + n)k(t) ) 0 = sk ∗ − (δ + n)k ∗ and so the steady-state values are: 1 s 1−α k ∗ = (5) δ + n α s 1−α y ∗ = : (6) δ + n Ceteris paribus, countries that have high savings/investment rates will tend to be richer since they accumulate more per capita capital and, as a result, have more per capita output. In turn, countries that have high population growth rates will tend to be poorer. O. Afonso, P. B. Vasconcelos Computational Economics 10 / 27 Economic model Golden rule The optimal capital accumulation leads to the golden rule savings rate that maximises the steady-state level of consumption. The per capita consumption c(t) is given by c(t) = (1 − s)y(t) and s can be y(t) c(t) written as s = y−(t) . _ At the steady state, k(t) = 0, that is, sy ∗ = (δ + n)k ∗. The latter can be written in terms of c(t), as c∗ = y ∗ − (δ + n)k ∗. dc To find k that maximises c, the first order condition, dk = 0, produces 1 α 1−α k gr = ; (7) δ + n gr considering A = 1. Note that k ∗ > k whenever s > α. O. Afonso, P. B. Vasconcelos Computational Economics 11 / 27 Economic model Golden rule f(k) (δ+n)k sf(k) sgrf(k) gr k k∗ Golden rule savings rate O. Afonso, P. B. Vasconcelos Computational Economics 12 / 27 Economic model Linear approximation The Solow model has a closed solution and can be demonstrated that the convergence process is (globally) stable. Taking (4), considering k_ (t) = φ(k(t)), a first order Taylor approximation for k_ (t), being φ(k(t)) at least twice differentiable, is given by _ dφ(k(t)) k(t) ' φ(k ∗) + dk(t) jk ∗ (k(t) − k ∗) ' sf 0(k ∗) − (n + δ) (k(t) − k ∗) since φ(k ∗) = 0. _ _ Thus, whenever k(t) < k ∗ (k(t) > k ∗) then k(t) > 0 (k(t) < 0) and physical capital will accumulate (diminish) towards the steady state. O. Afonso, P. B. Vasconcelos Computational Economics 13 / 27 Economic model Differential equation Since most differential equations are not analytically soluble, numerical solution of ordinary differential equations is a fundamental technique. A differential equation involves an unknown function, y(t), and its derivatives. A first–order ordinary differential equation, ODE, has the form dy y_ (t) = (t) = f (t; y(t)) (8) dt where f : R × R ! R and y(t):[t0; tT ] ⊂ R ! R. Equation (8) is non-autonomous but often economic problems are time-autonomous dy y_ (t) = (t) = f (y(t)): (9) dt The solution of (8) (or (9)) is a family of functions determined by a constant. O. Afonso, P. B. Vasconcelos Computational Economics 14 / 27 Economic model Initial value problems A particular solution is computed by requiring that it goes through a specific point, the initial condition, (t0; y0 = y(t0)). The problem specified both by (8) (or (9)) and the initial condition is called an initial value problem, IVP. Solving the IVP is to predict the path that a quantity will take during a certain time interval, given the initial quantity. Problems involving ODEs of a higher order can be reduced to a system of first ODE equations by introducing new variables. Seldom do these equations have solutions that can be expressed in a closed form or the analytical form is often too cumbersome; solution techniques are generally unable to deal with large and nonlinear systems of equations that arise in real problems. O. Afonso, P. B. Vasconcelos Computational Economics 15 / 27 Economic model Initial value problems To solve a continuous problem in a computer, a discretisation process is required. Initial value problems can be numerically solved using finite difference methods and recursive procedures. The numerical procedures are based on approximations y0; y1; ··· ; yT to the exact solution y(t0); y(t1); ··· ; y(tT ) at the grid points: t0 < t1 < ··· < tT .