Computation program files for "Pareto Distribution of Income in Neoclassical Growth Models" by Makoto Nirei and Shuhei Aoki, to appear in Review of Economic Dynamics

November 4, 2015

The computation program for the baseline case consists of the following source files.

main.f90
eqfuns.f90
global_vars.f90
myprocs.f90
params_V5.txt

The program execution file is compiled using Intel Fortran.  To compile, type "V5" for a command argument so that parameter values in "params_V5.txt" are incorporated. The compiled exe file stores output values in "results_V5.dat" and computed wealth distributions in "invdistk_V5.dat". The execution takes less than a day with an 8-core Windows Server PC. MATLAB is used to produce Figure 1 and distribution plots in Figures 3-7, using "plot_inv.m".

Various distributions shown in Figures 3-7 are produced similarly with altered command arguments or modified parameters in "global_vars.f90", as explained below.

Figure 3 is produced with altered arguments "V5rH", "V5rL", and "V5r0".

Figure 4 is produced with "V5whi" and "V5wlo".

Figure 5 is produced with "V5rip97" and "V5ripv7".

Figure 6 is produced by modifying the parameter values described in "global_vars.f90". Parameters mu, phib, or gamma need to be modified before the program is recompiled with altered arguments "V5mu0", "V5phiH", or "V5g4".

Figure 7 is produced with altered argument "V5t28".

There are no random numbers used for the computation results.