# Set up
# Before this setup, please go and check "README.md" to install what are necessary for running the following code.
from dolo import *
import dolark
from dolark import HModel # The model is written in yaml file, and HModel reads the yaml file.
from dolark.equilibrium import find_steady_state
from dolark.perturbation import perturb #find variations
from dolo import time_iteration, improved_time_iteration #time iteration is backward induction
from matplotlib import pyplot as plt
import numpy as np
import pandas as pdOptions for the following code. Setting Fast to true will set initial capital levels to the steady state level in the solver, dramatically increasing run time.
Fast = TrueDeclare parameters that do not change. These correspond to the values chosen in Aiyagari (1994).
# Declare these consistent parameters.
β = 0.96
δ = 0.08
α = 0.36Define functions that will be useful during calculation.
σAdjust = lambda ρ, σ: σ*(1-ρ**2)**0.5
# We need to set steady state values of labor.
# Since we know ρ and σ we can calculate this using properties of the AR1 process and lognormal distribution:
ΣUnconditional = lambda ρ, σ: (σ**2)/(1-ρ**2)
LMean = lambda ρ, σ: np.exp(0 + ΣUnconditional(ρ, σAdjust(ρ, σ))/2)
rFunc = lambda k: α*(k)**(α - 1) - δ
wFunc = lambda k: (1-α)*(k)**(α)aggmodel = HModel('Aiyagari.yaml')
load_capital = [5.718341702347233, 5.998527915377517, 6.236659640412816, 5.727020642279832, 6.023762697282315, 6.277078594476787, 5.740838016083044, 6.066502939823299, 6.348284509611548, 5.746926103372787, 6.111043099717087, 6.452912003167148, 6.2559200120869, 6.811669970608172, 7.242157096783963, 6.316337432044385, 6.969322649437888, 7.476453776919827, 6.412728832934068, 7.238541044351281, 7.89109764156157, 6.500822965861592, 7.642509470947722, 8.67333741118939]
σ_values = [0.2, 0.4]
ρ_values = [0.0, 0.3, 0.6, 0.9]
θ_values = [1, 3, 5]
results = []
index = 0
for z in σ_values:
for y in ρ_values:
for x in θ_values:
aggmodel.model.set_calibration( θ = x )
aggmodel.model.set_calibration( ρ_param = y )
Σ_adjusted = σAdjust(y, z)**2
aggmodel.model.set_calibration( Σ_param = Σ_adjusted )
aggmodel.set_calibration( L = LMean(y, z) )
if Fast == True:
aggmodel.set_calibration( K = load_capital[index] )
eq = find_steady_state(aggmodel)
df = eq.as_df()
a = df['a']
e = df['e']
i = df['i']
μ = df['μ']
# calculate ratio of K/L
# note, we could use L = LMean(y, z) as before instead of np.inner(np.exp(e), μ)
k = np.inner(i,μ) / np.inner(np.exp(e), μ)
# calculate aggregate consumption using r and w
r = rFunc(k)
w = wFunc(k)
c = [(1+r)*a[j] + w*np.exp(e[j]) - i[j] for j in range(len(df))]
# calculate aggregate income using r and w
income = [(r+δ)*a[j] + w*np.exp(e[j]) for j in range(len(df))]
# aggregate consumption and innome
agg_c = np.inner(c,μ)
agg_inc = np.inner(income,μ)
saving = (1 - agg_c/agg_inc)*100 #convert to %
saving_rate = float("%.2f" % saving) #with 2 decimals
interest = r*100
interest_rate = float("%.4f" % interest) #with 4 decimals
steady_state_capital = np.inner(a,μ)
results.append([saving_rate, interest_rate, steady_state_capital]) # append the results
index = index + 1
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.099e+02 | 3.148e+00 | 5.146e-01 | 6.132e-02 | 1.874e-03 | 2.539e-01 | 12 | 0 |
| 1 | 1.142e+02 | 7.152e+00 | 7.501e-03 | 7.350e-04 | 1.648e-03 | 4.714e-01 | 18 | 0 |
| 2 | 5.484e+01 | 1.154e+01 | 7.129e-03 | 1.251e-03 | 1.614e-03 | 6.011e-01 | 34 | 0 |
| 3 | 2.613e+01 | 1.749e+01 | 7.118e-03 | 1.979e-03 | 1.608e-03 | 7.403e-01 | 55 | 0 |
| 4 | 1.383e+01 | 2.885e+01 | 7.337e-03 | 2.449e-03 | 1.602e-03 | 8.105e-01 | 70 | 0 |
| 5 | 7.448e+00 | 4.484e+01 | 7.425e-03 | 2.219e-03 | 1.712e-03 | 8.127e-01 | 62 | 0 |
| 6 | 3.830e+00 | 5.941e+01 | 7.206e-03 | 2.182e-03 | 1.745e-03 | 8.035e-01 | 62 | 0 |
| 7 | 1.793e+00 | 6.083e+01 | 6.897e-03 | 2.305e-03 | 1.532e-03 | 8.220e-01 | 67 | 0 |
| 8 | 6.982e-01 | 4.060e+01 | 6.729e-03 | 1.515e-03 | 1.511e-03 | 4.875e-01 | 44 | 0 |
| 9 | 1.835e-01 | 1.393e+01 | 7.006e-03 | 1.846e-03 | 1.637e-03 | 2.955e-01 | 47 | 0 |
| 10 | 2.613e-02 | 2.149e+00 | 7.315e-03 | 1.439e-03 | 1.510e-03 | 3.798e-01 | 42 | 0 |
| 11 | 1.429e-03 | 9.094e-02 | 7.050e-03 | 1.097e-03 | 1.687e-03 | 3.911e-01 | 30 | 0 |
| 12 | 6.271e-06 | 2.025e-04 | 7.159e-03 | 6.049e-04 | 1.600e-03 | 3.410e-01 | 14 | 0 |
| 13 | 3.312e-10 | 4.489e-09 | 7.201e-03 | 1.469e-04 | 1.657e-03 | 4.352e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 1.28 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.112e+02 | 3.211e+00 | 2.118e-02 | 1.575e-03 | 1.434e-02 | 2.539e-01 | 12 | 0 |
| 1 | 1.152e+02 | 7.277e+00 | 4.315e-02 | 8.631e-04 | 1.906e-03 | 4.784e-01 | 20 | 0 |
| 2 | 5.533e+01 | 1.173e+01 | 8.845e-03 | 1.379e-03 | 1.976e-03 | 6.172e-01 | 36 | 0 |
| 3 | 2.639e+01 | 1.780e+01 | 8.590e-03 | 2.208e-03 | 1.950e-03 | 7.616e-01 | 61 | 0 |
| 4 | 1.399e+01 | 2.934e+01 | 8.605e-03 | 3.314e-03 | 1.907e-03 | 8.507e-01 | 94 | 0 |
| 5 | 7.543e+00 | 4.551e+01 | 8.802e-03 | 3.841e-03 | 1.957e-03 | 8.771e-01 | 104 | 0 |
| 6 | 3.885e+00 | 6.023e+01 | 8.888e-03 | 3.044e-03 | 1.997e-03 | 7.183e-01 | 82 | 0 |
| 7 | 1.824e+00 | 6.173e+01 | 8.741e-03 | 4.022e-03 | 1.882e-03 | 8.804e-01 | 109 | 0 |
| 8 | 7.122e-01 | 4.130e+01 | 9.427e-03 | 3.897e-03 | 2.047e-03 | 8.398e-01 | 105 | 0 |
| 9 | 1.868e-01 | 1.452e+01 | 9.049e-03 | 3.970e-03 | 3.026e-03 | 6.729e-01 | 112 | 0 |
| 10 | 2.816e-02 | 3.010e+00 | 1.121e-02 | 4.210e-03 | 2.012e-03 | 6.893e-01 | 118 | 0 |
| 11 | 3.583e-03 | 4.086e-01 | 8.869e-03 | 3.663e-03 | 2.076e-03 | 7.100e-01 | 102 | 0 |
| 12 | 1.541e-04 | 1.275e-02 | 8.509e-03 | 2.415e-03 | 2.184e-03 | 7.377e-01 | 67 | 0 |
| 13 | 2.806e-07 | 1.438e-05 | 8.489e-03 | 9.589e-04 | 2.006e-03 | 7.418e-01 | 24 | 0 |
| 14 | 8.021e-10 | 1.243e-08 | 8.603e-03 | 1.574e-04 | 1.927e-03 | 7.153e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.38 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.115e+02 | 3.224e+00 | 2.703e-02 | 1.558e-03 | 2.555e-02 | 2.539e-01 | 12 | 0 |
| 1 | 1.154e+02 | 7.302e+00 | 8.227e-03 | 8.249e-04 | 1.904e-03 | 4.798e-01 | 21 | 0 |
| 2 | 5.542e+01 | 1.177e+01 | 8.029e-03 | 1.416e-03 | 3.023e-03 | 6.204e-01 | 37 | 0 |
| 3 | 2.644e+01 | 1.786e+01 | 1.120e-02 | 2.262e-03 | 1.986e-03 | 7.656e-01 | 62 | 0 |
| 4 | 1.402e+01 | 2.944e+01 | 8.463e-03 | 3.531e-03 | 2.039e-03 | 8.582e-01 | 100 | 0 |
| 5 | 7.563e+00 | 4.564e+01 | 8.280e-03 | 4.412e-03 | 2.146e-03 | 8.932e-01 | 124 | 0 |
| 6 | 3.897e+00 | 6.040e+01 | 9.661e-03 | 3.936e-03 | 2.013e-03 | 6.676e-01 | 105 | 0 |
| 7 | 1.831e+00 | 6.189e+01 | 8.576e-03 | 4.732e-03 | 3.142e-03 | 8.619e-01 | 134 | 0 |
| 8 | 7.156e-01 | 4.142e+01 | 1.029e-02 | 5.455e-03 | 2.010e-03 | 8.711e-01 | 148 | 0 |
| 9 | 1.882e-01 | 1.461e+01 | 8.487e-03 | 5.866e-03 | 2.462e-03 | 7.505e-01 | 162 | 0 |
| 10 | 2.876e-02 | 3.218e+00 | 9.179e-03 | 6.641e-03 | 2.298e-03 | 8.010e-01 | 181 | 0 |
| 11 | 4.268e-03 | 5.944e-01 | 8.852e-03 | 6.173e-03 | 1.889e-03 | 8.210e-01 | 173 | 0 |
| 12 | 3.670e-04 | 4.162e-02 | 8.305e-03 | 4.425e-03 | 1.862e-03 | 8.400e-01 | 128 | 0 |
| 13 | 3.866e-06 | 2.607e-04 | 8.285e-03 | 2.361e-03 | 1.927e-03 | 8.431e-01 | 64 | 0 |
| 14 | 1.448e-09 | 1.942e-08 | 8.265e-03 | 1.547e-04 | 1.884e-03 | 7.237e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.37 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.156e+02 | 3.202e+00 | 2.679e-02 | 1.533e-03 | 4.911e-03 | 2.499e-01 | 12 | 0 |
| 1 | 1.167e+02 | 7.290e+00 | 2.247e-02 | 2.297e-03 | 5.236e-03 | 4.767e-01 | 20 | 0 |
| 2 | 5.618e+01 | 1.168e+01 | 2.225e-02 | 1.358e-03 | 1.944e-03 | 6.112e-01 | 35 | 0 |
| 3 | 2.679e+01 | 1.778e+01 | 8.401e-03 | 2.080e-03 | 1.851e-03 | 7.437e-01 | 56 | 0 |
| 4 | 1.418e+01 | 2.939e+01 | 8.139e-03 | 2.598e-03 | 1.916e-03 | 8.129e-01 | 71 | 0 |
| 5 | 7.669e+00 | 4.551e+01 | 8.257e-03 | 1.925e-03 | 1.940e-03 | 7.007e-01 | 53 | 0 |
| 6 | 3.969e+00 | 5.978e+01 | 8.199e-03 | 2.234e-03 | 1.918e-03 | 8.035e-01 | 61 | 0 |
| 7 | 1.878e+00 | 6.107e+01 | 8.269e-03 | 2.328e-03 | 1.886e-03 | 8.284e-01 | 67 | 0 |
| 8 | 7.443e-01 | 4.140e+01 | 8.394e-03 | 1.683e-03 | 1.868e-03 | 9.382e-01 | 46 | 0 |
| 9 | 1.962e-01 | 1.414e+01 | 8.106e-03 | 1.753e-03 | 1.904e-03 | 3.210e-01 | 47 | 0 |
| 10 | 2.770e-02 | 2.178e+00 | 8.279e-03 | 1.497e-03 | 1.797e-03 | 3.831e-01 | 42 | 0 |
| 11 | 1.602e-03 | 9.041e-02 | 8.150e-03 | 1.093e-03 | 1.862e-03 | 4.037e-01 | 30 | 0 |
| 12 | 4.770e-06 | 1.854e-04 | 8.074e-03 | 5.572e-04 | 1.921e-03 | 4.255e-01 | 14 | 0 |
| 13 | 2.303e-10 | 4.423e-09 | 8.124e-03 | 1.349e-04 | 1.908e-03 | 3.338e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.33 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.169e+02 | 3.265e+00 | 2.254e-02 | 1.685e-03 | 4.803e-03 | 2.499e-01 | 12 | 0 |
| 1 | 1.176e+02 | 7.412e+00 | 2.108e-02 | 2.270e-03 | 4.556e-03 | 4.839e-01 | 20 | 0 |
| 2 | 5.665e+01 | 1.187e+01 | 8.860e-03 | 1.361e-03 | 1.912e-03 | 6.295e-01 | 37 | 0 |
| 3 | 2.705e+01 | 1.809e+01 | 7.988e-03 | 2.251e-03 | 1.868e-03 | 7.655e-01 | 62 | 0 |
| 4 | 1.434e+01 | 2.987e+01 | 8.136e-03 | 3.329e-03 | 1.929e-03 | 8.530e-01 | 96 | 0 |
| 5 | 7.764e+00 | 4.617e+01 | 8.178e-03 | 3.784e-03 | 1.900e-03 | 8.820e-01 | 110 | 0 |
| 6 | 4.025e+00 | 6.060e+01 | 8.154e-03 | 3.629e-03 | 1.887e-03 | 8.820e-01 | 106 | 0 |
| 7 | 1.911e+00 | 6.197e+01 | 8.146e-03 | 3.331e-03 | 1.882e-03 | 7.016e-01 | 97 | 0 |
| 8 | 7.602e-01 | 4.213e+01 | 7.992e-03 | 3.674e-03 | 1.906e-03 | 8.394e-01 | 106 | 0 |
| 9 | 1.996e-01 | 1.472e+01 | 8.055e-03 | 3.949e-03 | 1.845e-03 | 6.684e-01 | 114 | 0 |
| 10 | 3.128e-02 | 3.063e+00 | 8.087e-03 | 4.116e-03 | 2.047e-03 | 7.044e-01 | 118 | 0 |
| 11 | 3.957e-03 | 4.070e-01 | 8.443e-03 | 3.486e-03 | 1.928e-03 | 7.184e-01 | 102 | 0 |
| 12 | 1.676e-04 | 1.288e-02 | 8.238e-03 | 2.401e-03 | 1.870e-03 | 7.483e-01 | 67 | 0 |
| 13 | 2.879e-07 | 1.434e-05 | 8.002e-03 | 8.721e-04 | 1.898e-03 | 7.420e-01 | 24 | 0 |
| 14 | 7.409e-10 | 1.194e-08 | 7.883e-03 | 1.390e-04 | 1.868e-03 | 7.838e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.35 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.172e+02 | 3.278e+00 | 2.409e-02 | 1.728e-03 | 4.344e-03 | 2.499e-01 | 12 | 0 |
| 1 | 1.178e+02 | 7.437e+00 | 2.376e-02 | 2.351e-03 | 4.361e-03 | 4.854e-01 | 21 | 0 |
| 2 | 5.674e+01 | 1.191e+01 | 1.361e-02 | 1.381e-03 | 1.833e-03 | 6.331e-01 | 37 | 0 |
| 3 | 2.710e+01 | 1.815e+01 | 8.062e-03 | 2.281e-03 | 1.797e-03 | 7.695e-01 | 63 | 0 |
| 4 | 1.437e+01 | 2.996e+01 | 8.168e-03 | 3.579e-03 | 1.886e-03 | 8.603e-01 | 101 | 0 |
| 5 | 7.783e+00 | 4.630e+01 | 8.137e-03 | 4.534e-03 | 1.818e-03 | 8.968e-01 | 130 | 0 |
| 6 | 4.037e+00 | 6.076e+01 | 7.925e-03 | 4.791e-03 | 1.879e-03 | 9.078e-01 | 136 | 0 |
| 7 | 1.917e+00 | 6.213e+01 | 8.262e-03 | 4.887e-03 | 1.870e-03 | 8.792e-01 | 139 | 0 |
| 8 | 7.634e-01 | 4.224e+01 | 8.118e-03 | 5.136e-03 | 1.801e-03 | 8.804e-01 | 148 | 0 |
| 9 | 2.001e-01 | 1.479e+01 | 8.132e-03 | 5.646e-03 | 1.815e-03 | 7.885e-01 | 165 | 0 |
| 10 | 3.191e-02 | 3.279e+00 | 9.309e-03 | 6.326e-03 | 3.015e-03 | 8.107e-01 | 183 | 0 |
| 11 | 4.884e-03 | 5.956e-01 | 8.619e-03 | 6.188e-03 | 3.003e-03 | 8.265e-01 | 173 | 0 |
| 12 | 3.835e-04 | 4.228e-02 | 1.354e-02 | 4.687e-03 | 2.202e-03 | 8.417e-01 | 128 | 0 |
| 13 | 3.462e-06 | 2.638e-04 | 1.108e-02 | 2.330e-03 | 3.099e-03 | 8.535e-01 | 64 | 0 |
| 14 | 1.496e-09 | 2.079e-08 | 1.223e-02 | 1.762e-04 | 2.136e-03 | 7.859e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.38 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.214e+02 | 3.270e+00 | 2.286e-02 | 1.599e-03 | 5.267e-03 | 2.509e-01 | 13 | 0 |
| 1 | 1.194e+02 | 7.475e+00 | 2.407e-02 | 2.636e-03 | 6.127e-03 | 4.855e-01 | 20 | 0 |
| 2 | 5.752e+01 | 1.191e+01 | 2.226e-02 | 3.860e-03 | 3.319e-03 | 6.326e-01 | 36 | 0 |
| 3 | 2.745e+01 | 1.822e+01 | 1.246e-02 | 2.061e-03 | 3.063e-03 | 7.485e-01 | 56 | 0 |
| 4 | 1.451e+01 | 3.016e+01 | 1.322e-02 | 2.758e-03 | 2.523e-03 | 8.178e-01 | 74 | 0 |
| 5 | 7.860e+00 | 4.642e+01 | 1.287e-02 | 2.430e-03 | 2.762e-03 | 8.253e-01 | 67 | 0 |
| 6 | 4.082e+00 | 6.029e+01 | 1.312e-02 | 2.343e-03 | 3.060e-03 | 8.209e-01 | 63 | 0 |
| 7 | 1.944e+00 | 6.135e+01 | 1.174e-02 | 2.443e-03 | 1.886e-03 | 8.376e-01 | 67 | 0 |
| 8 | 7.809e-01 | 4.236e+01 | 8.147e-03 | 1.700e-03 | 1.887e-03 | 7.077e-01 | 47 | 0 |
| 9 | 2.120e-01 | 1.469e+01 | 1.215e-02 | 1.695e-03 | 1.875e-03 | 3.594e-01 | 47 | 0 |
| 10 | 2.554e-02 | 2.142e+00 | 8.062e-03 | 1.559e-03 | 1.922e-03 | 3.902e-01 | 42 | 0 |
| 11 | 1.397e-03 | 8.885e-02 | 8.607e-03 | 1.167e-03 | 2.014e-03 | 4.033e-01 | 30 | 0 |
| 12 | 6.464e-06 | 2.116e-04 | 8.350e-03 | 5.851e-04 | 1.859e-03 | 4.263e-01 | 14 | 0 |
| 13 | 2.978e-10 | 5.600e-09 | 8.315e-03 | 1.280e-04 | 1.911e-03 | 3.248e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.35 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.226e+02 | 3.334e+00 | 2.586e-02 | 1.664e-03 | 7.329e-03 | 2.512e-01 | 13 | 0 |
| 1 | 1.204e+02 | 7.597e+00 | 2.577e-02 | 2.428e-03 | 4.732e-03 | 4.930e-01 | 21 | 0 |
| 2 | 5.798e+01 | 1.209e+01 | 2.146e-02 | 1.572e-03 | 3.042e-03 | 6.545e-01 | 39 | 0 |
| 3 | 2.771e+01 | 1.853e+01 | 1.261e-02 | 2.335e-03 | 2.118e-03 | 7.757e-01 | 64 | 0 |
| 4 | 1.467e+01 | 3.063e+01 | 1.277e-02 | 3.498e-03 | 2.998e-03 | 8.581e-01 | 99 | 0 |
| 5 | 7.956e+00 | 4.707e+01 | 1.296e-02 | 4.214e-03 | 3.028e-03 | 8.891e-01 | 119 | 0 |
| 6 | 4.140e+00 | 6.109e+01 | 1.226e-02 | 4.191e-03 | 1.907e-03 | 8.956e-01 | 120 | 0 |
| 7 | 1.977e+00 | 6.226e+01 | 7.974e-03 | 3.763e-03 | 1.941e-03 | 7.640e-01 | 104 | 0 |
| 8 | 7.979e-01 | 4.311e+01 | 7.989e-03 | 3.729e-03 | 1.851e-03 | 8.079e-01 | 108 | 0 |
| 9 | 2.159e-01 | 1.524e+01 | 8.155e-03 | 4.048e-03 | 1.919e-03 | 7.166e-01 | 116 | 0 |
| 10 | 2.958e-02 | 3.061e+00 | 8.229e-03 | 4.136e-03 | 1.749e-03 | 7.181e-01 | 120 | 0 |
| 11 | 3.889e-03 | 4.031e-01 | 8.155e-03 | 3.611e-03 | 1.824e-03 | 7.268e-01 | 103 | 0 |
| 12 | 1.606e-04 | 1.288e-02 | 8.225e-03 | 2.414e-03 | 1.898e-03 | 7.458e-01 | 68 | 0 |
| 13 | 3.046e-07 | 1.459e-05 | 7.917e-03 | 9.449e-04 | 1.953e-03 | 7.405e-01 | 25 | 0 |
| 14 | 6.224e-10 | 9.958e-09 | 8.269e-03 | 1.352e-04 | 1.887e-03 | 7.695e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.38 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.229e+02 | 3.347e+00 | 2.613e-02 | 1.721e-03 | 6.519e-03 | 2.512e-01 | 13 | 0 |
| 1 | 1.206e+02 | 7.622e+00 | 3.213e-02 | 9.012e-04 | 2.865e-03 | 4.945e-01 | 21 | 0 |
| 2 | 5.808e+01 | 1.213e+01 | 1.327e-02 | 1.542e-03 | 2.977e-03 | 6.587e-01 | 40 | 0 |
| 3 | 2.776e+01 | 1.860e+01 | 1.323e-02 | 2.416e-03 | 3.030e-03 | 7.798e-01 | 66 | 0 |
| 4 | 1.471e+01 | 3.072e+01 | 1.345e-02 | 4.002e-03 | 3.088e-03 | 8.651e-01 | 105 | 0 |
| 5 | 7.976e+00 | 4.720e+01 | 1.347e-02 | 5.171e-03 | 2.747e-03 | 9.024e-01 | 138 | 0 |
| 6 | 4.151e+00 | 6.124e+01 | 1.293e-02 | 5.784e-03 | 2.162e-03 | 9.156e-01 | 155 | 0 |
| 7 | 1.984e+00 | 6.243e+01 | 9.187e-03 | 5.643e-03 | 2.200e-03 | 9.093e-01 | 155 | 0 |
| 8 | 8.012e-01 | 4.324e+01 | 8.666e-03 | 5.610e-03 | 1.980e-03 | 8.635e-01 | 155 | 0 |
| 9 | 2.166e-01 | 1.531e+01 | 9.150e-03 | 6.100e-03 | 1.961e-03 | 8.311e-01 | 171 | 0 |
| 10 | 3.057e-02 | 3.283e+00 | 8.113e-03 | 6.624e-03 | 1.996e-03 | 8.287e-01 | 186 | 0 |
| 11 | 4.524e-03 | 5.887e-01 | 8.269e-03 | 6.218e-03 | 1.959e-03 | 8.335e-01 | 175 | 0 |
| 12 | 3.247e-04 | 4.171e-02 | 8.135e-03 | 4.636e-03 | 1.820e-03 | 8.433e-01 | 129 | 0 |
| 13 | 3.119e-06 | 2.612e-04 | 8.225e-03 | 2.310e-03 | 1.869e-03 | 8.412e-01 | 64 | 0 |
| 14 | 1.303e-09 | 1.756e-08 | 8.086e-03 | 1.512e-04 | 1.909e-03 | 7.410e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.41 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.267e+02 | 3.335e+00 | 2.391e-02 | 1.706e-03 | 5.200e-03 | 2.846e-01 | 14 | 0 |
| 1 | 1.226e+02 | 7.669e+00 | 2.712e-02 | 2.480e-03 | 4.676e-03 | 4.910e-01 | 22 | 0 |
| 2 | 5.887e+01 | 1.221e+01 | 2.270e-02 | 3.652e-03 | 3.191e-03 | 6.549e-01 | 33 | 0 |
| 3 | 2.808e+01 | 1.887e+01 | 8.137e-03 | 1.910e-03 | 1.858e-03 | 7.445e-01 | 55 | 0 |
| 4 | 1.479e+01 | 3.130e+01 | 8.171e-03 | 2.517e-03 | 1.831e-03 | 8.111e-01 | 68 | 0 |
| 5 | 7.994e+00 | 4.776e+01 | 7.946e-03 | 2.423e-03 | 1.837e-03 | 8.133e-01 | 67 | 0 |
| 6 | 4.147e+00 | 6.106e+01 | 7.981e-03 | 2.426e-03 | 1.910e-03 | 8.116e-01 | 67 | 0 |
| 7 | 1.974e+00 | 6.140e+01 | 8.194e-03 | 2.372e-03 | 1.873e-03 | 8.410e-01 | 67 | 0 |
| 8 | 7.953e-01 | 4.322e+01 | 8.322e-03 | 1.785e-03 | 1.917e-03 | 5.643e-01 | 49 | 0 |
| 9 | 2.228e-01 | 1.566e+01 | 8.619e-03 | 1.832e-03 | 3.017e-03 | 3.901e-01 | 49 | 0 |
| 10 | 3.077e-02 | 2.316e+00 | 1.007e-02 | 1.586e-03 | 1.936e-03 | 3.977e-01 | 43 | 0 |
| 11 | 1.732e-03 | 8.986e-02 | 8.245e-03 | 1.169e-03 | 1.875e-03 | 3.945e-01 | 30 | 0 |
| 12 | 7.573e-06 | 2.051e-04 | 8.178e-03 | 6.013e-04 | 1.830e-03 | 4.516e-01 | 14 | 0 |
| 13 | 1.846e-10 | 4.656e-09 | 8.162e-03 | 1.349e-04 | 2.010e-03 | 7.057e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.33 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.279e+02 | 3.400e+00 | 2.489e-02 | 1.730e-03 | 2.775e-02 | 2.858e-01 | 14 | 0 |
| 1 | 1.235e+02 | 7.795e+00 | 2.192e-02 | 2.456e-03 | 4.974e-03 | 5.289e-01 | 22 | 0 |
| 2 | 5.935e+01 | 1.241e+01 | 1.202e-02 | 1.570e-03 | 1.866e-03 | 6.911e-01 | 42 | 0 |
| 3 | 2.834e+01 | 1.918e+01 | 8.370e-03 | 2.514e-03 | 1.845e-03 | 8.079e-01 | 68 | 0 |
| 4 | 1.495e+01 | 3.177e+01 | 8.434e-03 | 3.814e-03 | 2.943e-03 | 8.711e-01 | 105 | 0 |
| 5 | 8.091e+00 | 4.840e+01 | 1.124e-02 | 4.644e-03 | 1.933e-03 | 8.981e-01 | 130 | 0 |
| 6 | 4.206e+00 | 6.183e+01 | 8.204e-03 | 5.100e-03 | 2.052e-03 | 9.143e-01 | 145 | 0 |
| 7 | 2.007e+00 | 6.233e+01 | 8.512e-03 | 5.396e-03 | 1.913e-03 | 9.211e-01 | 147 | 0 |
| 8 | 8.112e-01 | 4.399e+01 | 8.319e-03 | 4.399e-03 | 1.911e-03 | 8.746e-01 | 125 | 0 |
| 9 | 2.267e-01 | 1.620e+01 | 8.008e-03 | 4.489e-03 | 1.895e-03 | 7.954e-01 | 127 | 0 |
| 10 | 3.113e-02 | 3.059e+00 | 7.997e-03 | 4.448e-03 | 1.918e-03 | 7.710e-01 | 126 | 0 |
| 11 | 3.717e-03 | 4.021e-01 | 8.113e-03 | 3.708e-03 | 1.882e-03 | 7.594e-01 | 106 | 0 |
| 12 | 1.623e-04 | 1.281e-02 | 8.112e-03 | 2.408e-03 | 1.908e-03 | 7.620e-01 | 69 | 0 |
| 13 | 2.918e-07 | 1.459e-05 | 8.201e-03 | 1.018e-03 | 2.147e-03 | 7.785e-01 | 25 | 0 |
| 14 | 7.724e-10 | 1.228e-08 | 8.155e-03 | 1.380e-04 | 1.755e-03 | 7.482e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.37 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.281e+02 | 3.413e+00 | 2.483e-02 | 1.713e-03 | 5.380e-03 | 2.861e-01 | 14 | 0 |
| 1 | 1.237e+02 | 7.820e+00 | 2.248e-02 | 2.517e-03 | 5.236e-03 | 5.224e-01 | 22 | 0 |
| 2 | 5.944e+01 | 1.245e+01 | 2.391e-02 | 4.342e-03 | 2.252e-03 | 6.952e-01 | 43 | 0 |
| 3 | 2.839e+01 | 1.925e+01 | 1.359e-02 | 2.724e-03 | 2.856e-03 | 8.154e-01 | 72 | 0 |
| 4 | 1.498e+01 | 3.186e+01 | 1.262e-02 | 4.031e-03 | 3.016e-03 | 8.843e-01 | 117 | 0 |
| 5 | 8.111e+00 | 4.853e+01 | 1.299e-02 | 5.602e-03 | 2.627e-03 | 9.178e-01 | 159 | 0 |
| 6 | 4.217e+00 | 6.197e+01 | 1.291e-02 | 6.637e-03 | 1.829e-03 | 9.326e-01 | 187 | 0 |
| 7 | 2.013e+00 | 6.250e+01 | 8.172e-03 | 6.980e-03 | 1.967e-03 | 9.376e-01 | 197 | 0 |
| 8 | 8.147e-01 | 4.410e+01 | 7.964e-03 | 6.740e-03 | 1.861e-03 | 9.231e-01 | 190 | 0 |
| 9 | 2.283e-01 | 1.628e+01 | 8.415e-03 | 6.677e-03 | 1.904e-03 | 8.923e-01 | 193 | 0 |
| 10 | 3.189e-02 | 3.257e+00 | 8.032e-03 | 6.889e-03 | 1.855e-03 | 8.743e-01 | 199 | 0 |
| 11 | 3.875e-03 | 5.905e-01 | 8.208e-03 | 6.349e-03 | 1.855e-03 | 8.635e-01 | 183 | 0 |
| 12 | 3.092e-04 | 4.218e-02 | 8.325e-03 | 4.596e-03 | 1.840e-03 | 8.582e-01 | 134 | 0 |
| 13 | 3.074e-06 | 2.746e-04 | 8.098e-03 | 2.376e-03 | 1.901e-03 | 8.617e-01 | 66 | 0 |
| 14 | 1.209e-09 | 1.825e-08 | 8.286e-03 | 1.531e-04 | 1.768e-03 | 9.268e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.38 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.147e+02 | 2.749e+00 | 3.264e-02 | 1.567e-03 | 6.278e-03 | 2.405e-01 | 12 | 0 |
| 1 | 1.199e+02 | 6.801e+00 | 1.283e-02 | 8.180e-04 | 3.019e-03 | 4.407e-01 | 19 | 0 |
| 2 | 5.790e+01 | 1.115e+01 | 1.257e-02 | 1.350e-03 | 3.089e-03 | 5.827e-01 | 33 | 0 |
| 3 | 2.746e+01 | 1.682e+01 | 1.343e-02 | 2.002e-03 | 3.078e-03 | 7.251e-01 | 52 | 0 |
| 4 | 1.443e+01 | 2.781e+01 | 1.335e-02 | 2.512e-03 | 2.893e-03 | 8.013e-01 | 69 | 0 |
| 5 | 7.790e+00 | 4.377e+01 | 1.180e-02 | 2.417e-03 | 1.836e-03 | 8.184e-01 | 68 | 0 |
| 6 | 4.025e+00 | 5.884e+01 | 8.346e-03 | 2.331e-03 | 1.884e-03 | 8.197e-01 | 64 | 0 |
| 7 | 1.902e+00 | 6.163e+01 | 7.928e-03 | 1.890e-03 | 1.901e-03 | 5.515e-01 | 52 | 0 |
| 8 | 7.540e-01 | 4.269e+01 | 7.728e-03 | 1.771e-03 | 2.958e-03 | 6.783e-01 | 49 | 0 |
| 9 | 2.062e-01 | 1.550e+01 | 1.093e-02 | 1.745e-03 | 1.859e-03 | 3.712e-01 | 48 | 0 |
| 10 | 3.057e-02 | 2.556e+00 | 8.080e-03 | 1.582e-03 | 1.822e-03 | 3.885e-01 | 43 | 0 |
| 11 | 1.971e-03 | 1.257e-01 | 8.028e-03 | 1.150e-03 | 1.957e-03 | 4.073e-01 | 31 | 0 |
| 12 | 1.198e-05 | 3.576e-04 | 8.173e-03 | 6.580e-04 | 1.914e-03 | 4.524e-01 | 15 | 0 |
| 13 | 4.555e-10 | 8.995e-09 | 8.041e-03 | 1.330e-04 | 1.826e-03 | 3.336e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.37 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.160e+02 | 2.806e+00 | 2.991e-02 | 1.594e-03 | 5.345e-03 | 2.399e-01 | 12 | 0 |
| 1 | 1.210e+02 | 6.924e+00 | 2.080e-02 | 8.850e-04 | 3.038e-03 | 4.464e-01 | 20 | 0 |
| 2 | 5.837e+01 | 1.133e+01 | 1.278e-02 | 1.304e-03 | 2.253e-03 | 5.934e-01 | 34 | 0 |
| 3 | 2.770e+01 | 1.712e+01 | 1.352e-02 | 2.081e-03 | 2.989e-03 | 7.417e-01 | 56 | 0 |
| 4 | 1.459e+01 | 2.830e+01 | 1.313e-02 | 3.111e-03 | 1.951e-03 | 8.305e-01 | 84 | 0 |
| 5 | 7.886e+00 | 4.443e+01 | 8.239e-03 | 3.813e-03 | 1.935e-03 | 8.697e-01 | 107 | 0 |
| 6 | 4.081e+00 | 5.966e+01 | 8.173e-03 | 4.170e-03 | 1.930e-03 | 8.836e-01 | 117 | 0 |
| 7 | 1.932e+00 | 6.250e+01 | 8.130e-03 | 4.013e-03 | 1.938e-03 | 8.748e-01 | 113 | 0 |
| 8 | 7.678e-01 | 4.345e+01 | 8.297e-03 | 3.882e-03 | 2.017e-03 | 7.805e-01 | 108 | 0 |
| 9 | 2.091e-01 | 1.608e+01 | 8.040e-03 | 4.123e-03 | 1.946e-03 | 7.266e-01 | 115 | 0 |
| 10 | 3.251e-02 | 3.424e+00 | 8.051e-03 | 4.146e-03 | 1.881e-03 | 7.222e-01 | 119 | 0 |
| 11 | 3.945e-03 | 4.889e-01 | 8.491e-03 | 3.748e-03 | 1.943e-03 | 7.270e-01 | 104 | 0 |
| 12 | 1.915e-04 | 1.786e-02 | 8.086e-03 | 2.519e-03 | 1.890e-03 | 7.352e-01 | 70 | 0 |
| 13 | 4.898e-07 | 2.835e-05 | 8.268e-03 | 1.091e-03 | 1.868e-03 | 7.932e-01 | 28 | 0 |
| 14 | 7.550e-10 | 1.209e-08 | 8.325e-03 | 1.223e-04 | 1.893e-03 | 7.469e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.39 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.163e+02 | 2.817e+00 | 2.299e-02 | 1.551e-03 | 4.860e-03 | 2.398e-01 | 12 | 0 |
| 1 | 1.212e+02 | 6.949e+00 | 2.176e-02 | 2.203e-03 | 4.742e-03 | 4.475e-01 | 20 | 0 |
| 2 | 5.846e+01 | 1.137e+01 | 3.625e-02 | 3.556e-03 | 5.047e-03 | 5.955e-01 | 34 | 0 |
| 3 | 2.775e+01 | 1.718e+01 | 1.370e-02 | 2.062e-03 | 1.939e-03 | 7.449e-01 | 57 | 0 |
| 4 | 1.462e+01 | 2.840e+01 | 7.935e-03 | 3.195e-03 | 1.899e-03 | 8.361e-01 | 88 | 0 |
| 5 | 7.904e+00 | 4.457e+01 | 8.143e-03 | 4.225e-03 | 1.940e-03 | 8.807e-01 | 118 | 0 |
| 6 | 4.092e+00 | 5.982e+01 | 8.055e-03 | 5.154e-03 | 1.971e-03 | 9.036e-01 | 145 | 0 |
| 7 | 1.937e+00 | 6.268e+01 | 8.230e-03 | 5.503e-03 | 1.877e-03 | 9.086e-01 | 156 | 0 |
| 8 | 7.698e-01 | 4.357e+01 | 8.173e-03 | 5.370e-03 | 1.858e-03 | 8.758e-01 | 155 | 0 |
| 9 | 2.092e-01 | 1.615e+01 | 8.448e-03 | 5.891e-03 | 1.898e-03 | 8.332e-01 | 169 | 0 |
| 10 | 3.257e-02 | 3.620e+00 | 8.246e-03 | 6.470e-03 | 1.996e-03 | 8.288e-01 | 184 | 0 |
| 11 | 4.543e-03 | 6.868e-01 | 8.339e-03 | 6.143e-03 | 1.976e-03 | 8.327e-01 | 175 | 0 |
| 12 | 3.989e-04 | 5.356e-02 | 8.026e-03 | 4.623e-03 | 3.061e-03 | 8.399e-01 | 132 | 0 |
| 13 | 4.810e-06 | 4.253e-04 | 1.090e-02 | 2.468e-03 | 2.054e-03 | 8.575e-01 | 69 | 0 |
| 14 | 1.682e-09 | 2.308e-08 | 8.242e-03 | 1.340e-04 | 1.931e-03 | 7.376e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.37 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.250e+02 | 2.837e+00 | 3.659e-02 | 1.821e-03 | 5.616e-03 | 2.289e-01 | 12 | 0 |
| 1 | 1.239e+02 | 7.026e+00 | 3.052e-02 | 2.306e-03 | 2.744e-03 | 4.490e-01 | 20 | 0 |
| 2 | 5.980e+01 | 1.142e+01 | 8.059e-03 | 1.294e-03 | 1.464e-03 | 5.910e-01 | 33 | 0 |
| 3 | 2.843e+01 | 1.731e+01 | 6.736e-03 | 1.977e-03 | 1.441e-03 | 7.301e-01 | 53 | 0 |
| 4 | 1.496e+01 | 2.875e+01 | 5.994e-03 | 2.644e-03 | 1.413e-03 | 8.045e-01 | 71 | 0 |
| 5 | 8.115e+00 | 4.497e+01 | 9.771e-03 | 2.831e-03 | 3.043e-03 | 8.257e-01 | 75 | 0 |
| 6 | 4.233e+00 | 5.957e+01 | 1.117e-02 | 2.602e-03 | 2.598e-03 | 8.280e-01 | 69 | 0 |
| 7 | 2.030e+00 | 6.172e+01 | 5.953e-03 | 1.829e-03 | 1.377e-03 | 3.100e-01 | 48 | 0 |
| 8 | 8.256e-01 | 4.371e+01 | 7.805e-03 | 1.900e-03 | 3.003e-03 | 6.183e-01 | 50 | 0 |
| 9 | 2.275e-01 | 1.578e+01 | 1.368e-02 | 1.883e-03 | 2.850e-03 | 4.180e-01 | 49 | 0 |
| 10 | 3.318e-02 | 2.592e+00 | 8.442e-03 | 1.670e-03 | 2.042e-03 | 4.148e-01 | 43 | 0 |
| 11 | 1.940e-03 | 1.190e-01 | 8.461e-03 | 1.276e-03 | 1.879e-03 | 4.080e-01 | 31 | 0 |
| 12 | 8.330e-06 | 3.686e-04 | 6.893e-03 | 7.076e-04 | 1.874e-03 | 4.433e-01 | 15 | 0 |
| 13 | 4.924e-10 | 9.610e-09 | 7.379e-03 | 1.132e-04 | 1.305e-03 | 3.746e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.35 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.262e+02 | 2.893e+00 | 2.894e-02 | 1.587e-03 | 8.205e-03 | 2.284e-01 | 12 | 0 |
| 1 | 1.249e+02 | 7.144e+00 | 1.514e-02 | 9.379e-04 | 3.035e-03 | 4.551e-01 | 21 | 0 |
| 2 | 6.024e+01 | 1.159e+01 | 1.329e-02 | 1.364e-03 | 2.984e-03 | 6.022e-01 | 35 | 0 |
| 3 | 2.867e+01 | 1.761e+01 | 1.356e-02 | 2.213e-03 | 2.966e-03 | 7.467e-01 | 57 | 0 |
| 4 | 1.511e+01 | 2.923e+01 | 1.375e-02 | 3.138e-03 | 3.007e-03 | 8.331e-01 | 86 | 0 |
| 5 | 8.211e+00 | 4.561e+01 | 1.357e-02 | 4.035e-03 | 3.014e-03 | 8.748e-01 | 112 | 0 |
| 6 | 4.289e+00 | 6.038e+01 | 1.344e-02 | 4.745e-03 | 3.083e-03 | 8.948e-01 | 131 | 0 |
| 7 | 2.062e+00 | 6.262e+01 | 1.032e-02 | 4.769e-03 | 1.956e-03 | 8.979e-01 | 132 | 0 |
| 8 | 8.420e-01 | 4.448e+01 | 1.366e-02 | 4.208e-03 | 2.181e-03 | 8.474e-01 | 120 | 0 |
| 9 | 2.308e-01 | 1.633e+01 | 8.530e-03 | 4.428e-03 | 2.200e-03 | 7.777e-01 | 121 | 0 |
| 10 | 3.789e-02 | 3.506e+00 | 8.957e-03 | 4.460e-03 | 2.998e-03 | 7.498e-01 | 122 | 0 |
| 11 | 4.986e-03 | 4.838e-01 | 1.110e-02 | 3.754e-03 | 2.103e-03 | 7.416e-01 | 105 | 0 |
| 12 | 2.384e-04 | 1.817e-02 | 8.243e-03 | 2.529e-03 | 1.939e-03 | 7.512e-01 | 70 | 0 |
| 13 | 5.010e-07 | 2.743e-05 | 8.515e-03 | 1.090e-03 | 1.896e-03 | 7.502e-01 | 28 | 0 |
| 14 | 6.332e-10 | 1.033e-08 | 9.210e-03 | 1.721e-04 | 2.126e-03 | 7.796e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.40 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.265e+02 | 2.904e+00 | 1.250e-02 | 6.409e-04 | 3.057e-03 | 2.282e-01 | 12 | 0 |
| 1 | 1.251e+02 | 7.168e+00 | 1.296e-02 | 8.969e-04 | 3.001e-03 | 4.563e-01 | 21 | 0 |
| 2 | 6.033e+01 | 1.163e+01 | 1.235e-02 | 1.401e-03 | 2.277e-03 | 6.044e-01 | 35 | 0 |
| 3 | 2.872e+01 | 1.767e+01 | 1.306e-02 | 2.159e-03 | 3.066e-03 | 7.499e-01 | 58 | 0 |
| 4 | 1.514e+01 | 2.933e+01 | 1.207e-02 | 3.232e-03 | 1.866e-03 | 8.385e-01 | 89 | 0 |
| 5 | 8.230e+00 | 4.574e+01 | 8.447e-03 | 4.362e-03 | 1.861e-03 | 8.847e-01 | 123 | 0 |
| 6 | 4.301e+00 | 6.054e+01 | 9.832e-03 | 5.579e-03 | 1.886e-03 | 9.123e-01 | 159 | 0 |
| 7 | 2.069e+00 | 6.279e+01 | 8.196e-03 | 6.396e-03 | 1.977e-03 | 9.239e-01 | 180 | 0 |
| 8 | 8.451e-01 | 4.461e+01 | 8.282e-03 | 6.279e-03 | 1.884e-03 | 9.144e-01 | 177 | 0 |
| 9 | 2.315e-01 | 1.640e+01 | 7.926e-03 | 6.233e-03 | 1.838e-03 | 8.724e-01 | 180 | 0 |
| 10 | 3.850e-02 | 3.711e+00 | 8.294e-03 | 6.553e-03 | 1.820e-03 | 8.530e-01 | 189 | 0 |
| 11 | 5.895e-03 | 6.825e-01 | 8.191e-03 | 6.311e-03 | 1.878e-03 | 8.465e-01 | 178 | 0 |
| 12 | 4.830e-04 | 5.431e-02 | 7.939e-03 | 4.665e-03 | 1.843e-03 | 8.472e-01 | 133 | 0 |
| 13 | 4.875e-06 | 4.205e-04 | 8.222e-03 | 2.430e-03 | 1.813e-03 | 8.563e-01 | 69 | 0 |
| 14 | 1.943e-09 | 3.605e-08 | 8.230e-03 | 1.700e-04 | 1.906e-03 | 7.843e-01 | 1 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.34 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.347e+02 | 3.058e+00 | 2.208e-02 | 1.491e-03 | 5.233e-03 | 2.297e-01 | 12 | 0 |
| 1 | 1.283e+02 | 7.431e+00 | 1.850e-02 | 8.130e-04 | 1.846e-03 | 4.664e-01 | 21 | 0 |
| 2 | 6.158e+01 | 1.194e+01 | 8.097e-03 | 1.352e-03 | 1.748e-03 | 6.136e-01 | 35 | 0 |
| 3 | 2.936e+01 | 1.830e+01 | 8.115e-03 | 2.132e-03 | 1.878e-03 | 7.444e-01 | 56 | 0 |
| 4 | 1.543e+01 | 3.052e+01 | 8.201e-03 | 2.681e-03 | 1.840e-03 | 8.126e-01 | 74 | 0 |
| 5 | 8.393e+00 | 4.702e+01 | 8.160e-03 | 3.065e-03 | 1.951e-03 | 8.365e-01 | 82 | 0 |
| 6 | 4.400e+00 | 6.070e+01 | 8.151e-03 | 2.887e-03 | 1.934e-03 | 8.421e-01 | 78 | 0 |
| 7 | 2.130e+00 | 6.195e+01 | 8.244e-03 | 2.347e-03 | 1.938e-03 | 8.340e-01 | 63 | 0 |
| 8 | 8.824e-01 | 4.490e+01 | 8.279e-03 | 1.962e-03 | 1.838e-03 | 5.949e-01 | 53 | 0 |
| 9 | 2.525e-01 | 1.649e+01 | 8.239e-03 | 1.968e-03 | 1.911e-03 | 4.856e-01 | 51 | 0 |
| 10 | 2.921e-02 | 2.445e+00 | 8.196e-03 | 1.635e-03 | 1.865e-03 | 4.411e-01 | 44 | 0 |
| 11 | 2.027e-03 | 1.168e-01 | 8.196e-03 | 1.198e-03 | 1.863e-03 | 4.248e-01 | 31 | 0 |
| 12 | 8.206e-06 | 3.359e-04 | 8.193e-03 | 6.371e-04 | 1.899e-03 | 4.273e-01 | 15 | 0 |
| 13 | 3.578e-10 | 6.377e-09 | 8.197e-03 | 1.531e-04 | 1.886e-03 | 4.100e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.32 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.357e+02 | 3.115e+00 | 1.663e-02 | 5.660e-04 | 2.025e-03 | 2.296e-01 | 12 | 0 |
| 1 | 1.293e+02 | 7.544e+00 | 8.841e-03 | 9.410e-04 | 2.015e-03 | 4.732e-01 | 22 | 0 |
| 2 | 6.204e+01 | 1.211e+01 | 8.452e-03 | 1.520e-03 | 2.047e-03 | 6.264e-01 | 37 | 0 |
| 3 | 2.961e+01 | 1.859e+01 | 8.462e-03 | 2.322e-03 | 1.901e-03 | 7.615e-01 | 61 | 0 |
| 4 | 1.559e+01 | 3.096e+01 | 8.316e-03 | 3.256e-03 | 2.082e-03 | 8.408e-01 | 90 | 0 |
| 5 | 8.491e+00 | 4.765e+01 | 8.673e-03 | 4.297e-03 | 1.996e-03 | 8.818e-01 | 119 | 0 |
| 6 | 4.459e+00 | 6.149e+01 | 8.135e-03 | 5.151e-03 | 2.036e-03 | 9.059e-01 | 147 | 0 |
| 7 | 2.165e+00 | 6.288e+01 | 8.144e-03 | 5.612e-03 | 2.051e-03 | 9.151e-01 | 157 | 0 |
| 8 | 9.013e-01 | 4.572e+01 | 8.323e-03 | 5.148e-03 | 1.927e-03 | 9.037e-01 | 144 | 0 |
| 9 | 2.582e-01 | 1.704e+01 | 8.227e-03 | 4.758e-03 | 1.861e-03 | 8.401e-01 | 133 | 0 |
| 10 | 3.542e-02 | 3.440e+00 | 8.244e-03 | 4.637e-03 | 2.010e-03 | 7.935e-01 | 128 | 0 |
| 11 | 4.409e-03 | 4.573e-01 | 8.334e-03 | 3.851e-03 | 1.924e-03 | 7.668e-01 | 108 | 0 |
| 12 | 2.089e-04 | 1.698e-02 | 8.222e-03 | 2.586e-03 | 1.871e-03 | 7.615e-01 | 71 | 0 |
| 13 | 4.891e-07 | 2.486e-05 | 8.232e-03 | 1.095e-03 | 1.886e-03 | 7.664e-01 | 27 | 0 |
| 14 | 7.582e-10 | 1.233e-08 | 8.324e-03 | 1.352e-04 | 1.966e-03 | 7.708e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.39 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.360e+02 | 3.127e+00 | 2.476e-02 | 5.240e-04 | 1.312e-03 | 2.296e-01 | 12 | 0 |
| 1 | 1.294e+02 | 7.567e+00 | 6.088e-03 | 9.089e-04 | 1.342e-03 | 4.745e-01 | 22 | 0 |
| 2 | 6.214e+01 | 1.215e+01 | 6.040e-03 | 1.440e-03 | 1.404e-03 | 6.289e-01 | 37 | 0 |
| 3 | 2.966e+01 | 1.865e+01 | 1.184e-02 | 2.410e-03 | 1.444e-03 | 7.647e-01 | 62 | 0 |
| 4 | 1.563e+01 | 3.105e+01 | 6.519e-03 | 3.385e-03 | 3.006e-03 | 8.460e-01 | 93 | 0 |
| 5 | 8.511e+00 | 4.777e+01 | 1.057e-02 | 4.705e-03 | 1.453e-03 | 8.904e-01 | 129 | 0 |
| 6 | 4.471e+00 | 6.165e+01 | 6.641e-03 | 6.091e-03 | 1.351e-03 | 9.197e-01 | 173 | 0 |
| 7 | 2.172e+00 | 6.306e+01 | 6.582e-03 | 7.430e-03 | 1.580e-03 | 9.349e-01 | 206 | 0 |
| 8 | 9.050e-01 | 4.586e+01 | 6.539e-03 | 7.381e-03 | 1.475e-03 | 9.372e-01 | 211 | 0 |
| 9 | 2.582e-01 | 1.709e+01 | 6.108e-03 | 7.234e-03 | 1.444e-03 | 9.138e-01 | 201 | 0 |
| 10 | 3.581e-02 | 3.659e+00 | 6.230e-03 | 7.045e-03 | 1.404e-03 | 8.854e-01 | 201 | 0 |
| 11 | 5.362e-03 | 6.563e-01 | 6.051e-03 | 6.516e-03 | 1.473e-03 | 8.681e-01 | 184 | 0 |
| 12 | 4.419e-04 | 5.262e-02 | 6.663e-03 | 4.822e-03 | 1.453e-03 | 8.591e-01 | 136 | 0 |
| 13 | 4.734e-06 | 4.055e-04 | 6.041e-03 | 2.561e-03 | 1.564e-03 | 8.576e-01 | 70 | 0 |
| 14 | 1.629e-09 | 3.240e-08 | 6.455e-03 | 1.521e-04 | 1.430e-03 | 8.168e-01 | 1 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.33 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.428e+02 | 3.443e+00 | 3.310e-02 | 2.239e-03 | 6.342e-03 | 3.074e-01 | 15 | 0 |
| 1 | 1.330e+02 | 8.076e+00 | 1.857e-02 | 1.014e-03 | 3.016e-03 | 5.122e-01 | 23 | 0 |
| 2 | 6.346e+01 | 1.277e+01 | 1.309e-02 | 1.545e-03 | 2.754e-03 | 6.659e-01 | 39 | 0 |
| 3 | 3.027e+01 | 2.012e+01 | 1.258e-02 | 2.232e-03 | 2.896e-03 | 7.654e-01 | 59 | 0 |
| 4 | 1.585e+01 | 3.353e+01 | 1.361e-02 | 2.738e-03 | 2.720e-03 | 8.185e-01 | 76 | 0 |
| 5 | 8.588e+00 | 5.037e+01 | 1.240e-02 | 3.134e-03 | 3.088e-03 | 8.455e-01 | 87 | 0 |
| 6 | 4.492e+00 | 6.222e+01 | 1.235e-02 | 3.067e-03 | 2.592e-03 | 8.548e-01 | 86 | 0 |
| 7 | 2.171e+00 | 6.181e+01 | 8.271e-03 | 2.671e-03 | 1.909e-03 | 8.403e-01 | 74 | 0 |
| 8 | 9.016e-01 | 4.589e+01 | 8.366e-03 | 2.248e-03 | 1.845e-03 | 6.682e-01 | 61 | 0 |
| 9 | 2.684e-01 | 1.782e+01 | 8.047e-03 | 2.018e-03 | 1.865e-03 | 5.607e-01 | 55 | 0 |
| 10 | 3.982e-02 | 2.726e+00 | 8.117e-03 | 1.738e-03 | 1.937e-03 | 4.888e-01 | 46 | 0 |
| 11 | 1.782e-03 | 1.058e-01 | 8.222e-03 | 1.221e-03 | 1.909e-03 | 4.548e-01 | 31 | 0 |
| 12 | 1.055e-05 | 2.808e-04 | 8.306e-03 | 6.840e-04 | 1.844e-03 | 4.999e-01 | 15 | 0 |
| 13 | 2.820e-10 | 3.679e-09 | 8.360e-03 | 1.230e-04 | 1.933e-03 | 4.894e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.37 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.438e+02 | 3.508e+00 | 2.527e-02 | 1.824e-03 | 9.751e-03 | 3.090e-01 | 15 | 0 |
| 1 | 1.340e+02 | 8.192e+00 | 2.522e-02 | 1.006e-03 | 3.031e-03 | 5.120e-01 | 24 | 0 |
| 2 | 6.397e+01 | 1.297e+01 | 1.264e-02 | 1.642e-03 | 2.511e-03 | 6.861e-01 | 43 | 0 |
| 3 | 3.055e+01 | 2.044e+01 | 8.072e-03 | 2.616e-03 | 1.874e-03 | 8.016e-01 | 71 | 0 |
| 4 | 1.602e+01 | 3.396e+01 | 8.130e-03 | 3.778e-03 | 1.820e-03 | 8.682e-01 | 105 | 0 |
| 5 | 8.693e+00 | 5.094e+01 | 8.164e-03 | 4.806e-03 | 1.928e-03 | 8.992e-01 | 136 | 0 |
| 6 | 4.555e+00 | 6.298e+01 | 7.998e-03 | 5.940e-03 | 1.858e-03 | 9.202e-01 | 170 | 0 |
| 7 | 2.208e+00 | 6.277e+01 | 8.284e-03 | 6.796e-03 | 1.871e-03 | 9.327e-01 | 195 | 0 |
| 8 | 9.209e-01 | 4.678e+01 | 8.071e-03 | 6.861e-03 | 1.792e-03 | 9.364e-01 | 195 | 0 |
| 9 | 2.736e-01 | 1.841e+01 | 7.991e-03 | 5.974e-03 | 1.835e-03 | 9.175e-01 | 171 | 0 |
| 10 | 3.848e-02 | 3.368e+00 | 8.146e-03 | 5.362e-03 | 1.795e-03 | 8.741e-01 | 151 | 0 |
| 11 | 3.474e-03 | 4.465e-01 | 8.195e-03 | 4.238e-03 | 1.852e-03 | 8.348e-01 | 121 | 0 |
| 12 | 1.684e-04 | 1.625e-02 | 8.188e-03 | 2.671e-03 | 1.793e-03 | 8.051e-01 | 77 | 0 |
| 13 | 4.364e-07 | 2.499e-05 | 8.181e-03 | 1.091e-03 | 1.807e-03 | 7.876e-01 | 29 | 0 |
| 14 | 7.739e-10 | 1.280e-08 | 8.110e-03 | 1.571e-04 | 1.803e-03 | 7.822e-01 | 0 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.39 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
Computing Initial Rule... ------------------------------------------------------------------------------------------------------
| Start Improved Time Iterations. |
------------------------------------------------------------------------------------------------------
| N | f_x | d_x | Time_residuals | Time_inversion | Time_search | Lambda_0 | N_invert | N_search |
------------------------------------------------------------------------------------------------------
| 0 | 2.440e+02 | 3.521e+00 | 2.630e-02 | 8.383e-04 | 4.908e-03 | 3.093e-01 | 15 | 0 |
| 1 | 1.342e+02 | 8.215e+00 | 3.108e-02 | 2.908e-03 | 5.473e-03 | 5.120e-01 | 25 | 0 |
| 2 | 6.408e+01 | 1.301e+01 | 2.147e-02 | 1.764e-03 | 3.060e-03 | 6.901e-01 | 44 | 0 |
| 3 | 3.060e+01 | 2.050e+01 | 1.291e-02 | 2.762e-03 | 3.048e-03 | 8.048e-01 | 72 | 0 |
| 4 | 1.606e+01 | 3.405e+01 | 1.254e-02 | 4.028e-03 | 3.003e-03 | 8.730e-01 | 110 | 0 |
| 5 | 8.714e+00 | 5.105e+01 | 1.286e-02 | 5.362e-03 | 3.071e-03 | 9.060e-01 | 147 | 0 |
| 6 | 4.567e+00 | 6.313e+01 | 1.354e-02 | 6.728e-03 | 2.125e-03 | 9.292e-01 | 193 | 0 |
| 7 | 2.215e+00 | 6.295e+01 | 8.193e-03 | 8.587e-03 | 1.962e-03 | 9.448e-01 | 240 | 0 |
| 8 | 9.249e-01 | 4.691e+01 | 8.116e-03 | 9.378e-03 | 1.851e-03 | 9.530e-01 | 269 | 0 |
| 9 | 2.761e-01 | 1.849e+01 | 8.013e-03 | 9.455e-03 | 1.937e-03 | 9.528e-01 | 265 | 0 |
| 10 | 3.814e-02 | 3.530e+00 | 8.278e-03 | 8.726e-03 | 1.903e-03 | 9.372e-01 | 246 | 0 |
| 11 | 4.509e-03 | 6.470e-01 | 8.181e-03 | 7.642e-03 | 2.149e-03 | 9.169e-01 | 214 | 0 |
| 12 | 4.581e-04 | 5.020e-02 | 8.207e-03 | 5.490e-03 | 1.951e-03 | 8.975e-01 | 153 | 0 |
| 13 | 7.583e-06 | 3.956e-04 | 7.834e-03 | 2.732e-03 | 1.929e-03 | 8.787e-01 | 76 | 0 |
| 14 | 2.222e-09 | 3.454e-08 | 8.079e-03 | 1.690e-04 | 1.902e-03 | 7.733e-01 | 1 | 0 |
------------------------------------------------------------------------------------------------------
| Elapsed: 0.41 seconds. |
------------------------------------------------------------------------------------------------------
done
Computing Steady State...done
aiyagari = [[23.67, 4.1666], [23.71, 4.1456], [23.83, 4.0858], [23.73, 4.1365], [23.91, 4.0432], [24.19, 3.9054], [23.82, 4.0912], [24.25, 3.8767], [24.86, 3.5857], [24.14, 3.9305], [25.51, 3.2903], [27.36, 2.5260], [23.87, 4.0649], [24.44, 3.7816], [25.22, 3.4177], [24.09, 3.9554], [25.22, 3.4188], [26.66, 2.8032], [24.5, 3.7567], [26.71, 2.7835], [29.37, 1.807], [25.47, 3.3054], [31.0, 1.2894], [37.63, -0.3456]]
conv = lambda x, y, xd, yd: ("%."+str(xd)+"f") % x + "/" + ("%."+str(yd)+"f") % y
row_names = {
0:'ρ = 0',
1:'ρ = 0.3',
2:'ρ = 0.6',
3:'ρ = 0.9'
}
def OutputTables(title, subtitle, operation_func, xd = 4, yd = 2):
index = 0
print(title)
for sigma in range(1, 3):
output = []
for rho in range(1, 5):
output_slice = []
for theta in range(1, 4):
r1, r2 = operation_func(results[index][0], results[index][1], aiyagari[index][0], aiyagari[index][1])
if not r2 is None:
output_slice.append(conv(r1, r2, xd, yd))
else:
output_slice.append(str(r1))
index = index + 1
output.append(output_slice)
df = pd.DataFrame(output, columns=["µ = 1", "µ = 3", "µ = 5"])
df = df.rename(index = row_names)
print(subtitle + " (σ = " + str(0.2*sigma) + ")")
print(df)
print("\n")
def aiyagari_only(agg_s, r, A_agg_s, A_r):
return A_r, A_agg_s
def results_only(agg_s, r, A_agg_s, A_r):
return r, agg_s
OutputTables("Table II, Aiyagari (1994)", "Net return to capital in % / aggregate saving rate in %", aiyagari_only)
OutputTables("Table II, Replication (DolARK)", "Net return to capital in % / aggregate saving rate in %", results_only)
def differences(agg_s, r, A_agg_s, A_r):
return r - A_r, agg_s - A_agg_s
OutputTables("DolARK minus Aiyagari (1994)", "Net return to capital in % / aggregate saving rate in %", differences)
def differences_percent(agg_s, r, A_agg_s, A_r):
return 100*(r - A_r)/A_r, 100*(agg_s - A_agg_s)/A_agg_s
OutputTables("DolARK minus Aiyagari (1994), Percent Error", "Net return to capital in % / aggregate saving rate in %", differences_percent, 2, 2)
Table II, Aiyagari (1994)
Net return to capital in % / aggregate saving rate in % (σ = 0.2)
µ = 1 µ = 3 µ = 5
ρ = 0 4.1666/23.67 4.1456/23.71 4.0858/23.83
ρ = 0.3 4.1365/23.73 4.0432/23.91 3.9054/24.19
ρ = 0.6 4.0912/23.82 3.8767/24.25 3.5857/24.86
ρ = 0.9 3.9305/24.14 3.2903/25.51 2.5260/27.36
Net return to capital in % / aggregate saving rate in % (σ = 0.4)
µ = 1 µ = 3 µ = 5
ρ = 0 4.0649/23.87 3.7816/24.44 3.4177/25.22
ρ = 0.3 3.9554/24.09 3.4188/25.22 2.8032/26.66
ρ = 0.6 3.7567/24.50 2.7835/26.71 1.8070/29.37
ρ = 0.9 3.3054/25.47 1.2894/31.00 -0.3456/37.63
Table II, Replication (DolARK)
Net return to capital in % / aggregate saving rate in % (σ = 0.2)
µ = 1 µ = 3 µ = 5
ρ = 0 3.9451/24.11 3.5850/24.86 3.2999/25.49
ρ = 0.3 3.9335/24.13 3.5539/24.93 3.2533/25.59
ρ = 0.6 3.9151/24.17 3.5017/25.04 3.1723/25.78
ρ = 0.9 3.9071/24.19 3.4480/25.16 3.0560/26.05
Net return to capital in % / aggregate saving rate in % (σ = 0.4)
µ = 1 µ = 3 µ = 5
ρ = 0 3.7117/24.59 3.0909/25.97 2.6643/27.01
ρ = 0.3 3.6399/24.74 2.9296/26.35 2.4492/27.56
ρ = 0.6 3.5276/24.98 2.6677/27.00 2.0944/28.53
ρ = 0.9 3.4274/25.20 2.3033/27.95 1.5019/30.31
DolARK minus Aiyagari (1994)
Net return to capital in % / aggregate saving rate in % (σ = 0.2)
µ = 1 µ = 3 µ = 5
ρ = 0 -0.2215/0.44 -0.5606/1.15 -0.7859/1.66
ρ = 0.3 -0.2030/0.40 -0.4893/1.02 -0.6521/1.40
ρ = 0.6 -0.1761/0.35 -0.3750/0.79 -0.4134/0.92
ρ = 0.9 -0.0234/0.05 0.1577/-0.35 0.5300/-1.31
Net return to capital in % / aggregate saving rate in % (σ = 0.4)
µ = 1 µ = 3 µ = 5
ρ = 0 -0.3532/0.72 -0.6907/1.53 -0.7534/1.79
ρ = 0.3 -0.3155/0.65 -0.4892/1.13 -0.3540/0.90
ρ = 0.6 -0.2291/0.48 -0.1158/0.29 0.2874/-0.84
ρ = 0.9 0.1220/-0.27 1.0139/-3.05 1.8475/-7.32
DolARK minus Aiyagari (1994), Percent Error
Net return to capital in % / aggregate saving rate in % (σ = 0.2)
µ = 1 µ = 3 µ = 5
ρ = 0 -5.32/1.86 -13.52/4.85 -19.23/6.97
ρ = 0.3 -4.91/1.69 -12.10/4.27 -16.70/5.79
ρ = 0.6 -4.30/1.47 -9.67/3.26 -11.53/3.70
ρ = 0.9 -0.60/0.21 4.79/-1.37 20.98/-4.79
Net return to capital in % / aggregate saving rate in % (σ = 0.4)
µ = 1 µ = 3 µ = 5
ρ = 0 -8.69/3.02 -18.26/6.26 -22.04/7.10
ρ = 0.3 -7.98/2.70 -14.31/4.48 -12.63/3.38
ρ = 0.6 -6.10/1.96 -4.16/1.09 15.90/-2.86
ρ = 0.9 3.69/-1.06 78.63/-9.84 -534.58/-19.45
df = eq.as_df()
#print(1 + df['r'][0])
#print(1 + 0.36*(1.0833/40)**(1-0.36) - 0.08)
#print(df['w'][0])
#print((1-0.36)*(40/1.0833)**(0.36))
# extract variables from the steady state solution
a = df['a']
#r = df['r']
#w = df['w']
e = df['e']
i = df['i']
μ = df['μ']
agg_i = np.inner(i,μ)
r = 0.36*(1.0833/agg_i)**(1-0.36) - 0.08
w = (1-0.36)*(agg_i/1.0833)**(0.36)
# calculate consumption
c = [(1+r)*a[j] + w*np.exp(e[j]) - i[j] for j in range(len(df))]
# calculate income
income = [(r+0.08)*a[j] + w*np.exp(e[j]) for j in range(len(df))] #depreciation is 8%
print(np.inner(a,μ))
print("aggregate L", np.inner(np.exp(e), μ))
print("aggregate K", np.inner(i, μ))
agg_c = np.inner(c,μ)
agg_inc = np.inner(income,μ)
print("aggregate cons", agg_c)
print("aggregate inc", agg_inc)
# saving rate
saving = 1 - agg_c/agg_inc
print("R", 1 + r)
print("agg saving rate", saving*100)8.67333574917682
aggregate L 1.0821562193840921
aggregate K 8.673335749176845
aggregate cons 1.595380105351616
aggregate inc 2.2892469652857876
R 1.0150829297421882
agg saving rate 30.309829846058133