This directory contains programs for the numerical results in
``The optimum quantity of debt,'' which is joint work with
Rao Aiyagari.  A brief description of the files is given below.
For more details, see the paper or the comments in optimal.f.



compss                Compilation instructions
ilut.f                Sparse matrix routines called by optimal.f
matvec.f              Sparse matrix routines called by optimal.f
optimal.f             Main FORTRAN77 code 
optimal.f.bench       Backup of main FORTRAN77 code used for reported results
optimal.inp.bench     Input file for benchmark economy 
optimal.inp.beta      Input file for perturbation in beta (so r=6%)
optimal.inp.del5      Input file for perturbation in delta (delta=5%)
optimal.inp.del6      Input file for perturbation in delta (delta=6%)
optimal.inp.eta       Input file for perturbation in eta (so elasticity=1)
optimal.inp.lump      Input file for example with lump-sum taxes
optimal.inp.mu2       Input file for perturbation in mu (mu=2)
optimal.inp.mu3       Input file for perturbation in mu (mu=3)
optimal.inp.rho       Input file for perturbation in rho (rho=.5)
optimal.inp.sig       Input file for perturbation in sigma (sigma=2.5)
optimal.inp.tra       Input file for example with transfers=0
optimal.out.bench     Output file for benchmark economy 
optimal.out.beta      Output file for perturbation in beta (so r=6%)
optimal.out.del5      Output file for perturbation in delta (delta=5%)
optimal.out.del6      Output file for perturbation in delta (delta=6%)
optimal.out.eta       Output file for perturbation in eta (so elasticity=1)
optimal.out.lump      Output file for example with lump-sum taxes
optimal.out.mu2       Output file for perturbation in mu (mu=2)
optimal.out.mu3       Output file for perturbation in mu (mu=3)
optimal.out.rho       Output file for perturbation in rho (rho=.5)
optimal.out.sig       Output file for perturbation in sigma (sigma=2.5)
optimal.out.tra       Output file for example with transfers=0
qcksrt3.f             Subroutine for sorting arrays
qgausl.f              Subroutine for Legendre quadrature