Test Shock Persistence and Variance Default Low Cmin (Save + Borr Distribution)
back to Fan's Dynamic Assets Repository Table of Content.
@seealso
- test speed: fsi_abz_ds_vecsv_speed
- test joint RANDOM: fsi_abz_ds_vecsv_joint_rand
- test price no default CROSS: fsi_abz_ds_vecsv_price_nbc_cross
- test price default CROSS: fsi_abz_ds_vecsv_price_default_cross
- test interest rate no default: fsi_abz_ds_vecsv_nbc
- test interest rate no default CROSS: fsi_abz_ds_vecsv_nbc_cross
- test interest rate no default GRID: fsi_abz_ds_vecsv_nbc_grid
- test interest rate default: fsi_abz_ds_vecsv_default
- test interest rate default V(a,z) Comparison: fsi_abz_ds_vecsv_default_compaz
- test interest rate default CROSS: fsi_abz_ds_vecsv_default_cross
- test interest rate default GRID: fsi_abz_ds_vecsv_default_grid
- test shock default (very low cmin): fsi_abz_ds_vecsv_shk_default_lowcmin
- test shock no default: fsi_abz_ds_vecsv_shk_nbc
- test shock no default CROSS: fsi_abz_ds_vecsv_shk_nbc_cross
- test shock default: fsi_abz_ds_vecsv_shk_default
- test shock default CROSS: fsi_abz_ds_vecsv_shk_default_cross
- test preference no default: fsi_abz_ds_vecsv_pref_nbc
- test preference no default CROSS: fsi_abz_ds_vecsv_pref_nbc_cross
- test preference default: fsi_abz_ds_vecsv_pref_default
- test preference default CROSS: fsi_abz_ds_vecsv_pref_default_cross
- test preference default (very low cmin): fsi_abz_ds_vecsv_pref_default_lowcmin
Contents
Set Shared Parameters
close all; clear all; ar_fl_z_wage_rho = [0.65, 0.80, 0.95]; ar_fl_z_wage_sig = [0.05, 0.20, 0.35]; % Accuracy % ar_it_a_n_hg = [750, 1250, 1750, 2250]; % ar_it_z_n_hg = [15, 19, 23, 27]; ar_it_a_n_hg = [750]; % Borrowing/Savings Parameters bl_default = true; fl_c_min = 10^-300; % cmin so low default exists but never chosen fl_b_bd = -20; fl_r_save = 0.02;
Simulate Model with Low Persistence
for fl_z_wage_sig = ar_fl_z_wage_sig disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp(['fl_z_wage_rho = ' num2str(ar_fl_z_wage_rho(1))]); disp(['fl_z_wage_sig = ' num2str(fl_z_wage_sig)]); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp(''); disp(''); disp(''); disp(''); % Call Default Parameters <https://fanwangecon.github.io/CodeDynaAsset/m_az/paramfunc/html/ffs_abz_set_default_param.html ffs_abz_set_default_param> bl_input_override = true; it_param_set = 9; [param_map, support_map] = ffs_abz_set_default_param(it_param_set); % Shock Parameters param_map('fl_z_wage_rho') = ar_fl_z_wage_rho(1); param_map('fl_z_wage_sig') = fl_z_wage_sig; % Borrowing Parameters param_map('bl_default') = bl_default; param_map('fl_c_min') = fl_c_min; param_map('fl_b_bd') = fl_b_bd; % Interest Rates param_map('fl_r_save') = fl_r_save; % Display Parameters support_map('bl_display') = false; support_map('bl_display_final') = false; support_map('bl_time') = true; support_map('bl_profile') = false; for it_accuracy = 1:length(ar_it_a_n_hg) % Accuracy Regular param_map('it_a_n') = ar_it_a_n_hg(it_accuracy); it_z_n = param_map('it_z_n'); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp(['it_a_n = ' num2str(ar_it_a_n_hg(it_accuracy)) ', it_z_n = ' num2str(it_z_n)]); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); % Call Grid Generator <https://fanwangecon.github.io/CodeDynaAsset/m_az/paramfunc/html/ffs_abz_get_funcgrid.html ffs_abz_get_funcgrid> [armt_map, func_map] = ffs_abz_get_funcgrid(param_map, support_map, bl_input_override); % Call Dynamic Programming Problem <https://fanwangecon.github.io/CodeDynaAsset/m_az/solve/html/ff_abz_vf_vecsv.html ff_abz_vf_vecsv> result_map = ff_abz_vf_vecsv(param_map, support_map, armt_map, func_map); % Call Distribution CProgram result_map = ff_az_ds_vecsv(param_map, support_map, armt_map, func_map, result_map, bl_input_override); end % Snap snapnow; end % close all close all;
xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx fl_z_wage_rho = 0.65 fl_z_wage_sig = 0.05 xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx it_a_n = 750, it_z_n = 55 xxxxxxxxxxxxxxxxxxxxxxxxxxx Elapsed time is 16.556212 seconds. Elapsed time is 1.549792 seconds. mean sd coefofvar min max pYis0 pYls0 pYgr0 pYisMINY pYisMAXY p0_1 p1 p5 p10 p15 p20 p25 p35 p50 p65 p75 p80 p85 p90 p95 p99 p99_9 fl_cov_cl_mt_val fl_cor_cl_mt_val fl_cov_cl_mt_pol_a fl_cor_cl_mt_pol_a fl_cov_cl_mt_coh fl_cor_cl_mt_coh fl_cov_cl_mt_pol_c fl_cor_cl_mt_pol_c fracByP0_1 fracByP1 fracByP5 fracByP10 fracByP15 fracByP20 fracByP25 fracByP35 fracByP50 fracByP65 fracByP75 fracByP80 fracByP85 fracByP90 fracByP95 fracByP99 fracByP99_9 ___________ ____________________ ____________________ _______ ______ _______ __________ ________ ___________ ___________ ________ ________ ________ _________ _________ ________ ______ ______ ______ ______ ______ ______ ______ ______ ______ ________ _______ ________________ ________________ __________________ ___________________ ________________ ___________________ __________________ __________________ __________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ ___________ cl_mt_val 2.6491e+115 0+5.863e+132i 0+2.2132e+17i -2e+150 14.088 0 1.4958e-30 1 -3.0169e-36 3.3242e-37 3.1966 3.3682 3.5148 3.6058 3.6434 3.6451 3.7345 3.7345 3.8246 3.9136 3.9136 3.9136 4.0015 4.0015 4.0881 4.1737 4.2782 -3.4374e+265 1+0i -2.9764e+116 0+1.1287e-15i -4.2009e+116 0+1.3998e-15i -3.2861e+115 0+6.835e-17i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 cl_mt_pol_a -0.013279 0.044977+0i -3.3871+0i -15.371 46.845 0.77461 0.18686 0.038527 -3.7937e-37 3.3242e-37 -0.36872 -0.20999 -0.1256 -0.045282 -0.044558 0 0 0 0 0 0 0 0 0 0 0.039164 0.12616 -2.9764e+116 0+1.1287e-15i 0.0020229 1+0i 0.00199 0.86436+0i 0.0018909 0.5127+0i 0.057408 0.21728 0.69715 0.87818 1.1433 1.1433 1.1433 1.1433 1.1433 1.1433 1.1433 1.1433 1.1433 1.1433 1.1433 1.041 1.0112 cl_mt_coh 0.096617 0.051188+0i 0.5298+0i -16.345 50.122 0 0.055434 0.94457 -4.3681e-37 -5.0766e-37 -0.28215 -0.12707 -0.01593 0.050409 0.072809 0.072809 0.0768 0.0992 0.1216 0.1216 0.1216 0.1216 0.1216 0.1216 0.1216 0.16155 0.25029 -4.2009e+116 0+1.3998e-15i 0.00199 0.86436+0i 0.0026202 1+0i 0.0017442 0.41554+0i -0.0063696 -0.021302 -0.040506 -0.019263 0.051693 0.051693 0.099474 0.30572 0.93941 0.93941 0.93941 0.93941 0.93941 0.93941 0.93941 0.98985 0.99764 cl_mt_pol_c 1.2787 0.082002+0i 0.064127+0i 1e-300 4.7379 0 0 1 -3.0169e-36 2.484e-37 1.1021 1.1246 1.1682 1.179 1.1802 1.2199 1.2277 1.2277 1.2771 1.3285 1.3285 1.3285 1.382 1.382 1.4377 1.4564 1.5166 -3.2861e+115 0+6.835e-17i 0.0018909 0.5127+0i 0.0017442 0.41554+0i 0.0067244 1+0i 0.00089531 0.0099513 0.046265 0.096063 0.16324 0.18352 0.36365 0.36365 0.59148 0.78796 0.78796 0.78796 0.9107 0.9107 0.96756 0.98922 0.99991
xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx fl_z_wage_rho = 0.65 fl_z_wage_sig = 0.2 xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx it_a_n = 750, it_z_n = 55 xxxxxxxxxxxxxxxxxxxxxxxxxxx Elapsed time is 18.272003 seconds. Elapsed time is 2.310810 seconds. mean sd coefofvar min max pYis0 pYls0 pYgr0 pYisMINY pYisMAXY p0_1 p1 p5 p10 p15 p20 p25 p35 p50 p65 p75 p80 p85 p90 p95 p99 p99_9 fl_cov_cl_mt_val fl_cor_cl_mt_val fl_cov_cl_mt_pol_a fl_cor_cl_mt_pol_a fl_cov_cl_mt_coh fl_cor_cl_mt_coh fl_cov_cl_mt_pol_c fl_cor_cl_mt_pol_c fracByP0_1 fracByP1 fracByP5 fracByP10 fracByP15 fracByP20 fracByP25 fracByP35 fracByP50 fracByP65 fracByP75 fracByP80 fracByP85 fracByP90 fracByP95 fracByP99 fracByP99_9 ___________ ____________________ ___________________ _______ ______ ________ ________ _______ ___________ ___________ _______ ________ _______ _______ _______ ________ ________ ________ ________ ________ _______ _______ _______ _______ ______ ______ ______ ________________ ________________ __________________ __________________ ________________ ___________________ __________________ ___________________ __________ _________ ________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ ___________ cl_mt_val 2.4059e+117 0+6.1511e+133i 0+2.5566e+16i -2e+150 14.292 0 0.023778 0.97622 -8.1842e-34 -5.9934e-36 -3.7853 -0.96069 0.80488 1.5045 1.9242 2.2095 2.464 2.8216 3.2369 3.6334 3.8731 4.0625 4.2164 4.4584 4.791 5.4046 6.0554 -3.7837e+267 1+0i -1.3242e+118 0+2.1529e-16i -4.6913e+118 0+7.1853e-16i -2.7858e+117 0+1.8825e-16i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 cl_mt_pol_a -0.16902 0.99993+0i -5.9161+0i -19.512 47.92 0.089091 0.51344 0.39747 -2.57e-35 -5.9934e-36 -4.6089 -3.1225 -1.963 -1.3878 -1.1095 -0.8535 -0.61718 -0.3414 -0.08535 0.091625 0.3665 0.45813 0.64138 0.91625 1.3744 2.3823 3.5734 -1.3242e+118 0+2.1529e-16i 0.99987 1+0i 1.0369 0.977+0i 0.21829 0.9074+0i 0.031755 0.22635 0.78862 1.2743 1.6365 1.9386 2.1331 2.4588 2.6003 2.5743 2.3764 2.3076 2.1506 1.8802 1.545 1.1549 1.0243 cl_mt_coh -0.090576 1.0614+0i -11.719+0i -19.968 50.122 0 0.46915 0.53085 -2.0908e-34 4.5609e-35 -4.9251 -3.2653 -2.0279 -1.3961 -1.0934 -0.81298 -0.57741 -0.25223 0.028142 0.21506 0.47303 0.58889 0.77581 1.0562 1.5235 2.5291 3.7665 -4.6913e+118 0+7.1853e-16i 1.0369 0.977+0i 1.1266 1+0i 0.21218 0.83091+0i 0.062273 0.44635 1.5358 2.4607 3.1341 3.6798 4.0187 4.5309 4.6431 4.4253 4.0176 3.7277 3.3737 2.7897 2.093 1.3344 1.0474 cl_mt_pol_c 1.2473 0.24058+0i 0.19288+0i 1e-300 4.7989 0 0 1 -8.1842e-34 -5.9934e-36 0.63559 0.7853 0.90122 0.96272 1.0065 1.041 1.067 1.1087 1.2293 1.3353 1.4427 1.4586 1.519 1.5649 1.6579 1.8087 1.957 -2.7858e+117 0+1.8825e-16i 0.21829 0.9074+0i 0.21218 0.83091+0i 0.057879 1+0i 0.00046238 0.0057923 0.033246 0.07121 0.11072 0.15157 0.1942 0.2812 0.42058 0.57544 0.68592 0.74806 0.80869 0.86581 0.92974 0.98521 0.99841
xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx fl_z_wage_rho = 0.65 fl_z_wage_sig = 0.35 xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx it_a_n = 750, it_z_n = 55 xxxxxxxxxxxxxxxxxxxxxxxxxxx Elapsed time is 18.188948 seconds. Elapsed time is 4.189742 seconds. mean sd coefofvar min max pYis0 pYls0 pYgr0 pYisMINY pYisMAXY p0_1 p1 p5 p10 p15 p20 p25 p35 p50 p65 p75 p80 p85 p90 p95 p99 p99_9 fl_cov_cl_mt_val fl_cor_cl_mt_val fl_cov_cl_mt_pol_a fl_cor_cl_mt_pol_a fl_cov_cl_mt_coh fl_cor_cl_mt_coh fl_cov_cl_mt_pol_c fl_cor_cl_mt_pol_c fracByP0_1 fracByP1 fracByP5 fracByP10 fracByP15 fracByP20 fracByP25 fracByP35 fracByP50 fracByP65 fracByP75 fracByP80 fracByP85 fracByP90 fracByP95 fracByP99 fracByP99_9 ____________ ___________ ___________ _______ ______ ________ _______ _______ ___________ __________ _______ _______ _______ _________ _______ ________ ________ _________ _______ ______ ______ ______ ______ ______ ______ ______ ______ ________________ ________________ __________________ __________________ ________________ ________________ __________________ __________________ __________ _________ ________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ ___________ cl_mt_val -3.7608e+117 6.1774e+133 -1.6426e+16 -2e+150 14.527 0 0.10347 0.89653 -6.1498e-34 5.5769e-28 -6.783 -3.363 -1.1185 -0.063842 0.56795 1.0705 1.418 2.0574 2.8305 3.5237 4.0485 4.3456 4.6965 5.1533 5.8051 7.0174 8.2458 3.8161e+267 1 8.6071e+118 6.9902e-16 6.6564e+118 5.2265e-16 4.5381e+117 1.9602e-16 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 cl_mt_pol_a 0.71106 1.9933 2.8032 -19.512 49.02 0.037772 0.38438 0.57785 1.7795e-34 3.4459e-27 -2.9019 -2.6021 -1.9397 -1.4509 -1.058 -0.76815 -0.5121 -0.1707 0.3665 1.0079 1.5576 2.0158 2.5655 3.2985 4.4896 7.33 10.995 8.6071e+118 6.9902e-16 3.9731 1 4.0107 0.97598 0.69794 0.93431 -0.0067864 -0.042477 -0.16541 -0.28805 -0.37094 -0.4415 -0.49032 -0.53035 -0.50341 -0.36342 -0.19781 -0.057361 0.11066 0.30528 0.56118 0.87852 0.98338 cl_mt_coh 0.81114 2.0617 2.5417 -19.968 50.122 0 0.36503 0.63497 -1.0923e-34 2.2284e-27 -3.056 -2.7045 -1.9793 -1.4672 -1.0447 -0.71952 -0.46155 -0.065316 0.45063 1.1496 1.7104 2.1777 2.716 3.4637 4.701 7.5758 11.314 6.6564e+118 5.2265e-16 4.0107 0.97598 4.2506 1 0.66874 0.8655 -0.0054826 -0.038384 -0.15083 -0.26034 -0.33276 -0.39229 -0.42306 -0.45549 -0.42136 -0.26557 -0.10525 0.027837 0.16531 0.35257 0.60079 0.88687 0.98448 cl_mt_pol_c 1.2689 0.37477 0.29534 1e-300 5.5434 0 0 1 -6.1498e-34 5.5769e-28 0.27837 0.51638 0.72142 0.81429 0.87257 0.91411 0.97957 1.0721 1.2486 1.4338 1.5437 1.5748 1.6629 1.7699 1.8982 2.1895 2.5032 4.5381e+117 1.9602e-16 0.69794 0.93431 0.66874 0.8655 0.14045 1 0.00017555 0.0033808 0.023559 0.05405 0.087504 0.12415 0.16 0.24159 0.38512 0.53834 0.65466 0.71425 0.77931 0.84765 0.9184 0.9818 0.99798
Simulate Model with Medium Persistence
close all for fl_z_wage_sig = ar_fl_z_wage_sig disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp(['fl_z_wage_rho = ' num2str(ar_fl_z_wage_rho(2))]); disp(['fl_z_wage_sig = ' num2str(fl_z_wage_sig)]); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp(''); disp(''); disp(''); disp(''); % Call Default Parameters <https://fanwangecon.github.io/CodeDynaAsset/m_az/paramfunc/html/ffs_abz_set_default_param.html ffs_abz_set_default_param> bl_input_override = true; it_param_set = 9; [param_map, support_map] = ffs_abz_set_default_param(it_param_set); % Simulation Accuracy param_map('fl_z_wage_rho') = ar_fl_z_wage_rho(2); param_map('fl_z_wage_sig') = fl_z_wage_sig; % Borrowing Parameters param_map('bl_default') = bl_default; param_map('fl_c_min') = fl_c_min; param_map('fl_b_bd') = fl_b_bd; % Interest Rates param_map('fl_r_save') = fl_r_save; % Display Parameters support_map('bl_display') = false; support_map('bl_display_final') = false; support_map('bl_time') = true; support_map('bl_profile') = false; for it_accuracy = 1:length(ar_it_a_n_hg) % Accuracy Regular param_map('it_a_n') = ar_it_a_n_hg(it_accuracy); it_z_n = param_map('it_z_n'); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp(['it_a_n = ' num2str(ar_it_a_n_hg(it_accuracy)) ', it_z_n = ' num2str(it_z_n)]); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); % Call Grid Generator <https://fanwangecon.github.io/CodeDynaAsset/m_az/paramfunc/html/ffs_abz_get_funcgrid.html ffs_abz_get_funcgrid> [armt_map, func_map] = ffs_abz_get_funcgrid(param_map, support_map, bl_input_override); % Call Dynamic Programming Problem <https://fanwangecon.github.io/CodeDynaAsset/m_az/solve/html/ff_abz_vf_vecsv.html ff_abz_vf_vecsv> result_map = ff_abz_vf_vecsv(param_map, support_map, armt_map, func_map); % Call Distribution CProgram result_map = ff_az_ds_vecsv(param_map, support_map, armt_map, func_map, result_map, bl_input_override); end % Snap snapnow; end % close all close all;
xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx fl_z_wage_rho = 0.8 fl_z_wage_sig = 0.05 xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx it_a_n = 750, it_z_n = 55 xxxxxxxxxxxxxxxxxxxxxxxxxxx Elapsed time is 17.084215 seconds. Elapsed time is 1.321009 seconds. mean sd coefofvar min max pYis0 pYls0 pYgr0 pYisMINY pYisMAXY p0_1 p1 p5 p10 p15 p20 p25 p35 p50 p65 p75 p80 p85 p90 p95 p99 p99_9 fl_cov_cl_mt_val fl_cor_cl_mt_val fl_cov_cl_mt_pol_a fl_cor_cl_mt_pol_a fl_cov_cl_mt_coh fl_cor_cl_mt_coh fl_cov_cl_mt_pol_c fl_cor_cl_mt_pol_c fracByP0_1 fracByP1 fracByP5 fracByP10 fracByP15 fracByP20 fracByP25 fracByP35 fracByP50 fracByP65 fracByP75 fracByP80 fracByP85 fracByP90 fracByP95 fracByP99 fracByP99_9 ____________ ___________ ___________ _______ ______ _______ __________ ________ __________ ___________ ________ ________ _________ ______ ______ ______ ______ ______ ______ ______ ______ ______ ______ ______ ______ ________ ________ ________________ ________________ __________________ __________________ ________________ ________________ __________________ __________________ __________ _________ _________ __________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ ___________ cl_mt_val -6.1389e+118 2.9704e+134 -4.8387e+15 -2e+150 14.162 0 7.3369e-16 1 5.2038e-33 -6.3096e-34 2.6673 2.9549 3.2493 3.433 3.433 3.5622 3.6149 3.6149 3.7944 3.9711 3.9711 3.9711 4.145 4.145 4.3156 4.4926 4.6452 8.8234e+268 1 7.8274e+119 6.0673e-14 1.0233e+120 6.9205e-14 7.654e+118 2.3536e-15 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 cl_mt_pol_a -0.010202 0.043431 -4.2573 -16.29 46.905 0.87741 0.092287 0.030306 2.6941e-34 -6.3096e-34 -0.40304 -0.22598 -0.062807 0 0 0 0 0 0 0 0 0 0 0 0 0.021778 0.021778 7.8274e+119 6.0673e-14 0.0018863 1 0.0019469 0.90049 0.0019284 0.40555 0.052229 0.41091 0.85735 1.0647 1.0647 1.0647 1.0647 1.0647 1.0647 1.0647 1.0647 1.0647 1.0647 1.0647 1.0647 1 1 cl_mt_coh 0.10003 0.049781 0.49765 -17.204 50.122 0 0.037906 0.96209 7.0161e-34 4.075e-34 -0.32803 -0.12585 0.031526 0.0544 0.0768 0.0992 0.0992 0.0992 0.1216 0.1216 0.1216 0.1216 0.1216 0.1216 0.1216 0.14381 0.14381 1.0233e+120 6.9205e-14 0.0019469 0.90049 0.0024781 1 0.0019109 0.35062 -0.0056725 -0.029062 -0.034601 -0.0074475 0.045187 0.27258 0.27258 0.27258 0.97182 0.97182 0.97182 0.97182 0.97182 0.97182 0.97182 1 1 cl_mt_pol_c 1.2791 0.10948 0.085593 1e-300 4.7533 0 0 1 5.2038e-33 -1.9308e-33 1.0226 1.0761 1.0976 1.1539 1.1539 1.2072 1.2131 1.2131 1.2752 1.3406 1.3406 1.3406 1.4094 1.4094 1.4816 1.558 1.6379 7.654e+118 2.3536e-15 0.0019284 0.40555 0.0019109 0.35062 0.011986 1 0.00079797 0.0085144 0.047775 0.16891 0.16891 0.18529 0.36199 0.36199 0.58437 0.78302 0.78302 0.78302 0.90791 0.90791 0.96361 0.99498 1
xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx fl_z_wage_rho = 0.8 fl_z_wage_sig = 0.2 xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx it_a_n = 750, it_z_n = 55 xxxxxxxxxxxxxxxxxxxxxxxxxxx Elapsed time is 16.680868 seconds. Elapsed time is 2.687246 seconds. mean sd coefofvar min max pYis0 pYls0 pYgr0 pYisMINY pYisMAXY p0_1 p1 p5 p10 p15 p20 p25 p35 p50 p65 p75 p80 p85 p90 p95 p99 p99_9 fl_cov_cl_mt_val fl_cor_cl_mt_val fl_cov_cl_mt_pol_a fl_cor_cl_mt_pol_a fl_cov_cl_mt_coh fl_cor_cl_mt_coh fl_cov_cl_mt_pol_c fl_cor_cl_mt_pol_c fracByP0_1 fracByP1 fracByP5 fracByP10 fracByP15 fracByP20 fracByP25 fracByP35 fracByP50 fracByP65 fracByP75 fracByP80 fracByP85 fracByP90 fracByP95 fracByP99 fracByP99_9 ___________ ____________________ __________________ _______ ______ ________ _______ _______ ___________ ___________ _______ _______ _______ ________ _______ ________ ________ ________ ______ _______ _______ _______ ______ ______ ______ ______ ______ ________________ ________________ __________________ ___________________ ________________ ___________________ __________________ ___________________ __________ _________ ________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ ___________ cl_mt_val 3.3332e+117 0+7.1797e+133i 0+2.154e+16i -2e+150 14.589 0 0.11682 0.88318 2.3358e-34 2.3217e-35 -8.9682 -4.4505 -1.5339 -0.29673 0.43541 0.96182 1.3696 2.0596 2.8477 3.5671 4.0176 4.3221 4.6078 5.0699 5.6464 6.7665 7.8111 -5.1549e+267 1+0i -6.6998e+118 0+6.3651e-16i -6.0541e+118 0+5.4882e-16i -4.0285e+117 0+1.5715e-16i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 cl_mt_pol_a 0.081285 1.466+0i 18.036+0i -19.512 48.287 0.064257 0.45228 0.48346 -5.8876e-37 2.3217e-35 -4.3368 -3.4694 -2.2551 -1.6216 -1.1949 -0.88168 -0.6828 -0.2645 0 0.3665 0.733 0.91625 1.2828 1.7409 2.6571 4.6729 7.2384 -6.6998e+118 0+6.3651e-16i 2.1493 1+0i 2.2245 0.98761+0i 0.47335 0.90433+0i -0.060103 -0.49472 -1.847 -3.0761 -3.9742 -4.5336 -4.9931 -5.5481 -5.7523 -5.3447 -4.7081 -4.3056 -3.5848 -2.6978 -1.3206 0.31282 0.90849 cl_mt_coh 0.16285 1.5364+0i 9.4345+0i -19.968 50.122 0 0.4165 0.5835 3.0262e-35 -1.1302e-34 -4.5737 -3.6391 -2.3307 -1.6541 -1.1868 -0.85778 -0.62606 -0.20357 0.1216 0.47303 0.82446 1.0562 1.43 1.8973 2.8319 4.888 7.4824 -6.0541e+118 0+5.4882e-16i 2.2245 0.98761+0i 2.3606 1+0i 0.46457 0.84691+0i -0.032346 -0.26114 -0.96714 -1.5982 -2.0488 -2.3204 -2.5379 -2.7713 -2.7677 -2.509 -2.1467 -1.8356 -1.4336 -0.95014 -0.21421 0.64359 0.94956 cl_mt_pol_c 1.2504 0.35703+0i 0.28553+0i 1e-300 4.9919 0 0 1 2.3358e-34 2.3217e-35 0.35143 0.5776 0.73959 0.82878 0.88998 0.93052 0.95572 1.0675 1.2061 1.3852 1.4914 1.5904 1.6307 1.7413 1.8623 2.107 2.3622 -4.0285e+117 0+1.5715e-16i 0.47335 0.90433+0i 0.46457 0.84691+0i 0.12747 1+0i 0.00022741 0.0039804 0.025515 0.057067 0.091538 0.12933 0.16599 0.24874 0.39305 0.54713 0.654 0.71841 0.78211 0.84821 0.92031 0.9822 0.99812
xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx fl_z_wage_rho = 0.8 fl_z_wage_sig = 0.35 xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx it_a_n = 750, it_z_n = 55 xxxxxxxxxxxxxxxxxxxxxxxxxxx Elapsed time is 16.470775 seconds. Elapsed time is 4.874312 seconds. mean sd coefofvar min max pYis0 pYls0 pYgr0 pYisMINY pYisMAXY p0_1 p1 p5 p10 p15 p20 p25 p35 p50 p65 p75 p80 p85 p90 p95 p99 p99_9 fl_cov_cl_mt_val fl_cor_cl_mt_val fl_cov_cl_mt_pol_a fl_cor_cl_mt_pol_a fl_cov_cl_mt_coh fl_cor_cl_mt_coh fl_cov_cl_mt_pol_c fl_cor_cl_mt_pol_c fracByP0_1 fracByP1 fracByP5 fracByP10 fracByP15 fracByP20 fracByP25 fracByP35 fracByP50 fracByP65 fracByP75 fracByP80 fracByP85 fracByP90 fracByP95 fracByP99 fracByP99_9 ____________ ___________ __________ _______ ______ ________ _______ _______ ___________ __________ _______ _______ _______ ________ ________ _________ ________ _______ ______ ______ ______ ______ ______ ______ ______ ______ ______ ________________ ________________ __________________ __________________ ________________ ________________ __________________ __________________ __________ __________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ ___________ cl_mt_val -1.7497e+118 1.4195e+134 -8.113e+15 -2e+150 15.047 0 0.27006 0.72994 3.6063e-34 4.4335e-09 -14.956 -8.5924 -4.3992 -2.711 -1.5953 -0.82337 -0.21871 0.89567 2.2899 3.5224 4.4621 5.0126 5.6398 6.425 7.5495 9.5522 11.364 2.015e+268 1 2.9938e+119 5.8815e-16 3.2761e+119 6.2923e-16 2.2333e+118 2.8523e-16 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 cl_mt_pol_a 2.3597 3.5859 1.5196 -19.512 49.02 0.024367 0.23151 0.74413 -1.6764e-34 1.9363e-08 -1.7923 -1.6216 -1.1095 -0.69389 -0.35267 -0.1707 0 0.45813 1.1911 2.3823 3.5734 4.398 5.4059 6.9635 9.6207 15.76 23.823 2.9938e+119 5.8815e-16 12.859 1 13.013 0.98939 1.7814 0.90065 -0.0016775 -0.0074733 -0.031507 -0.049329 -0.060583 -0.066998 -0.067741 -0.054397 -0.007762 0.10594 0.23247 0.31711 0.41775 0.54985 0.72011 0.91923 0.98881 cl_mt_coh 2.5074 3.6678 1.4628 -19.968 50.122 0 0.21305 0.78695 -7.3797e-35 1.2522e-08 -1.841 -1.6541 -1.0934 -0.64846 -0.29703 -0.065316 0.1216 0.56649 1.3366 2.5515 3.7665 4.5852 5.6356 7.202 9.9347 16.196 24.398 3.2761e+119 6.2923e-16 13.013 0.98939 13.453 1 1.7268 0.85354 -0.001593 -0.0071796 -0.029859 -0.045906 -0.055251 -0.059161 -0.057937 -0.044033 0.010611 0.12639 0.25228 0.33046 0.43449 0.55944 0.72915 0.92203 0.98906 cl_mt_pol_c 1.3165 0.55157 0.41896 1e-300 7.1222 0 0 1 3.6063e-34 4.4335e-09 0.14448 0.33659 0.52081 0.66016 0.71157 0.82124 0.86416 1.0545 1.2762 1.4979 1.6362 1.749 1.8754 2.0405 2.2977 2.8199 3.4554 2.2333e+118 2.8523e-16 1.7814 0.90065 1.7268 0.85354 0.30423 1 8.262e-05 0.0019908 0.01603 0.03881 0.064967 0.094572 0.12674 0.20272 0.33406 0.49081 0.61065 0.67433 0.74347 0.81792 0.9001 0.9765 0.99724
Simulate Model with High Persistence
close all for fl_z_wage_sig = ar_fl_z_wage_sig disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp(['fl_z_wage_rho = ' num2str(ar_fl_z_wage_rho(3))]); disp(['fl_z_wage_sig = ' num2str(fl_z_wage_sig)]); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp(''); disp(''); disp(''); disp(''); % Call Default Parameters <https://fanwangecon.github.io/CodeDynaAsset/m_az/paramfunc/html/ffs_abz_set_default_param.html ffs_abz_set_default_param> bl_input_override = true; it_param_set = 9; [param_map, support_map] = ffs_abz_set_default_param(it_param_set); % Simulation Accuracy param_map('fl_z_wage_rho') = ar_fl_z_wage_rho(3); param_map('fl_z_wage_sig') = fl_z_wage_sig; % Borrowing Parameters param_map('bl_default') = bl_default; param_map('fl_c_min') = fl_c_min; param_map('fl_b_bd') = fl_b_bd; % Interest Rates param_map('fl_r_save') = fl_r_save; % Display Parameters support_map('bl_display') = false; support_map('bl_display_final') = false; support_map('bl_time') = true; support_map('bl_profile') = false; for it_accuracy = 1:length(ar_it_a_n_hg) % Accuracy Regular param_map('it_a_n') = ar_it_a_n_hg(it_accuracy); it_z_n = param_map('it_z_n'); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); disp(['it_a_n = ' num2str(ar_it_a_n_hg(it_accuracy)) ', it_z_n = ' num2str(it_z_n)]); disp('xxxxxxxxxxxxxxxxxxxxxxxxxxx'); % Call Grid Generator <https://fanwangecon.github.io/CodeDynaAsset/m_az/paramfunc/html/ffs_abz_get_funcgrid.html ffs_abz_get_funcgrid> [armt_map, func_map] = ffs_abz_get_funcgrid(param_map, support_map, bl_input_override); % Call Dynamic Programming Problem <https://fanwangecon.github.io/CodeDynaAsset/m_az/solve/html/ff_abz_vf_vecsv.html ff_abz_vf_vecsv> result_map = ff_abz_vf_vecsv(param_map, support_map, armt_map, func_map); % Call Distribution CProgram result_map = ff_az_ds_vecsv(param_map, support_map, armt_map, func_map, result_map, bl_input_override); end % Snap snapnow; end % close all close all; clear all;
xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx fl_z_wage_rho = 0.95 fl_z_wage_sig = 0.05 xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx it_a_n = 750, it_z_n = 55 xxxxxxxxxxxxxxxxxxxxxxxxxxx Elapsed time is 15.437944 seconds. Elapsed time is 1.243355 seconds. mean sd coefofvar min max pYis0 pYls0 pYgr0 pYisMINY pYisMAXY p0_1 p1 p5 p10 p15 p20 p25 p35 p50 p65 p75 p80 p85 p90 p95 p99 p99_9 fl_cov_cl_mt_val fl_cor_cl_mt_val fl_cov_cl_mt_pol_a fl_cor_cl_mt_pol_a fl_cov_cl_mt_coh fl_cor_cl_mt_coh fl_cov_cl_mt_pol_c fl_cor_cl_mt_pol_c fracByP0_1 fracByP1 fracByP5 fracByP10 fracByP15 fracByP20 fracByP25 fracByP35 fracByP50 fracByP65 fracByP75 fracByP80 fracByP85 fracByP90 fracByP95 fracByP99 fracByP99_9 ____________ ___________ ___________ _______ ______ _______ __________ ___________ ___________ ___________ _________ _______ _______ ______ ______ ______ ______ ______ ______ ______ ______ ______ ______ ______ ______ ______ ______ ________________ ________________ __________________ __________________ ________________ ________________ __________________ __________________ ___________ _________ ________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ ___________ cl_mt_val -1.7133e+118 1.3087e+134 -7.6385e+15 -2e+150 14.828 0 0.0067917 0.99321 2.5038e-35 2.1531e-34 -0.70882 0.17784 1.04 1.9025 1.9025 1.9025 2.744 2.744 3.5598 4.3486 4.3486 5.1094 5.1094 5.1094 5.8391 6.5231 7.0893 1.7128e+268 1 2.5642e+119 2.5519e-13 2.7519e+119 1.1609e-13 2.107e+118 6.9575e-16 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 cl_mt_pol_a -0.00061387 0.0076776 -12.507 -19.239 46.912 0.99329 0.0067125 -2.6624e-31 -6.9734e-35 2.1531e-34 -0.088168 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2.5642e+119 2.5519e-13 5.8946e-05 1 6.2019e-05 0.446 0.00025565 0.1439 0.26167 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 cl_mt_coh 0.11048 0.018112 0.16394 -19.968 50.122 0 0.00098481 0.99902 -1.4985e-34 -1.4218e-34 0.0057421 0.0544 0.0768 0.0768 0.0992 0.0992 0.0992 0.0992 0.1216 0.1216 0.1216 0.1216 0.1216 0.1216 0.1216 0.1216 0.1216 2.7519e+119 1.1609e-13 6.2019e-05 0.446 0.00032805 1 0.00027055 0.064554 -0.00031171 0.0087674 0.062348 0.062348 0.29304 0.29304 0.29304 0.29304 1 1 1 1 1 1 1 1 1 cl_mt_pol_c 1.2799 0.23139 0.18078 1e-300 5.124 0 0 1 2.5038e-35 2.1531e-34 0.77097 0.85764 0.94413 1.0393 1.0393 1.0393 1.1441 1.1441 1.2595 1.3865 1.3865 1.5264 1.5264 1.5264 1.6803 1.8497 2.0363 2.107e+118 6.9575e-16 0.00025565 0.1439 0.00027055 0.064554 0.053543 1 0.00067635 0.019135 0.064199 0.16247 0.16247 0.16247 0.32637 0.32637 0.53355 0.73218 0.73218 0.87658 0.87658 0.87658 0.95617 0.98941 1
xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx fl_z_wage_rho = 0.95 fl_z_wage_sig = 0.2 xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx it_a_n = 750, it_z_n = 55 xxxxxxxxxxxxxxxxxxxxxxxxxxx Elapsed time is 18.027697 seconds. Elapsed time is 5.636764 seconds. mean sd coefofvar min max pYis0 pYls0 pYgr0 pYisMINY pYisMAXY p0_1 p1 p5 p10 p15 p20 p25 p35 p50 p65 p75 p80 p85 p90 p95 p99 p99_9 fl_cov_cl_mt_val fl_cor_cl_mt_val fl_cov_cl_mt_pol_a fl_cor_cl_mt_pol_a fl_cov_cl_mt_coh fl_cor_cl_mt_coh fl_cov_cl_mt_pol_c fl_cor_cl_mt_pol_c fracByP0_1 fracByP1 fracByP5 fracByP10 fracByP15 fracByP20 fracByP25 fracByP35 fracByP50 fracByP65 fracByP75 fracByP80 fracByP85 fracByP90 fracByP95 fracByP99 fracByP99_9 ___________ ___________ ___________ _______ ______ _______ ________ _______ ___________ __________ ________ ________ ________ _______ _______ _______ _______ ________ ________ ______ ______ ______ ______ ______ ______ ______ ______ ________________ ________________ __________________ __________________ ________________ ________________ __________________ __________________ ___________ __________ _________ _________ __________ __________ __________ __________ _________ _________ _________ _________ _________ _________ _________ _________ ___________ cl_mt_val -1.033e+118 1.1767e+134 -1.1392e+16 -2e+150 17.104 0 0.53883 0.46117 5.5748e-34 1.5338e-05 -44.356 -28.192 -15.414 -10.854 -9.8964 -7.3289 -5.3755 -4.1309 -0.74554 2.6906 4.7653 5.98 6.9816 8.841 10.709 13.934 15.868 1.3847e+268 1 1.5721e+119 2.6773e-16 1.8381e+119 3.0677e-16 1.3331e+118 1.2519e-16 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 cl_mt_pol_a 2.6553 4.9902 1.8793 -19.512 49.02 0.25421 0.088094 0.6577 -4.6901e-35 2.2406e-05 -1.3656 -0.86736 -0.3414 0 0 0 0 0.091625 0.64138 1.7409 3.0236 4.1231 5.5892 8.063 12.736 24.647 38.208 1.5721e+119 2.6773e-16 24.902 1 25.352 0.99772 3.8342 0.84906 -0.00052479 -0.0045643 -0.013638 -0.015673 -0.015673 -0.015673 -0.015673 -0.014633 0.0063101 0.067548 0.15248 0.22289 0.31004 0.4395 0.62602 0.88592 0.98437 cl_mt_coh 2.8168 5.0919 1.8077 -19.968 50.122 0 0.073802 0.9262 -2.2381e-35 1.4489e-05 -1.3737 -0.83538 -0.25223 0.0768 0.0992 0.1216 0.1216 0.19266 0.75341 1.8749 3.2057 4.3048 5.8225 8.3011 13.112 25.195 39.071 1.8381e+119 3.0677e-16 25.352 0.99772 25.927 1 3.7977 0.82419 -0.00049782 -0.0042814 -0.011878 -0.011961 -0.0096451 -0.0025486 -0.0025486 -0.0018475 0.022377 0.08709 0.17523 0.2415 0.33068 0.45381 0.63845 0.88868 0.98479 cl_mt_pol_c 1.3303 0.90494 0.68024 1e-300 7.7508 0 0 1 5.5748e-34 1.5338e-05 0.091433 0.19127 0.30479 0.4433 0.45951 0.58138 0.67484 0.87458 1.0974 1.4584 1.7747 1.9594 2.1393 2.5527 3.1582 4.4157 5.6435 1.3331e+118 1.2519e-16 3.8342 0.84906 3.7977 0.82419 0.81892 1 4.9471e-05 0.001049 0.0092447 0.022603 0.047394 0.059034 0.10966 0.13844 0.25001 0.39475 0.51316 0.58304 0.66092 0.74876 0.85335 0.96263 0.99556
xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx fl_z_wage_rho = 0.95 fl_z_wage_sig = 0.35 xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxx it_a_n = 750, it_z_n = 55 xxxxxxxxxxxxxxxxxxxxxxxxxxx Elapsed time is 19.547910 seconds. Elapsed time is 26.326025 seconds. mean sd coefofvar min max pYis0 pYls0 pYgr0 pYisMINY pYisMAXY p0_1 p1 p5 p10 p15 p20 p25 p35 p50 p65 p75 p80 p85 p90 p95 p99 p99_9 fl_cov_cl_mt_val fl_cor_cl_mt_val fl_cov_cl_mt_pol_a fl_cor_cl_mt_pol_a fl_cov_cl_mt_coh fl_cor_cl_mt_coh fl_cov_cl_mt_pol_c fl_cor_cl_mt_pol_c fracByP0_1 fracByP1 fracByP5 fracByP10 fracByP15 fracByP20 fracByP25 fracByP35 fracByP50 fracByP65 fracByP75 fracByP80 fracByP85 fracByP90 fracByP95 fracByP99 fracByP99_9 ____________ ___________ ___________ _______ ______ ________ ___________ _______ ___________ _________ ________ ________ _______ ________ _______ _______ _______ _______ _______ _______ ______ ______ ______ ______ ______ ______ ______ ________________ ________________ __________________ __________________ ________________ ________________ __________________ __________________ __________ __________ __________ __________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ _________ ___________ cl_mt_val -4.6635e+116 4.5896e+133 -9.8414e+16 -2e+150 18.504 0 0.6318 0.3682 1.6665e-34 0.0035502 -94.117 -73.379 -46.067 -35.689 -24.71 -21.946 -18.268 -11.506 -4.4329 0.76569 4.7287 6.4686 8.7022 10.968 13.751 17.131 18.504 2.1064e+267 1 1.0694e+118 2.0208e-17 1.6582e+118 3.0718e-17 6.8818e+116 8.2212e-18 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 cl_mt_pol_a 8.7267 11.531 1.3214 -19.512 49.02 0.084293 -4.5445e-33 0.91571 -1.2308e-35 0.011406 0 0 0 0.091625 0.45813 0.733 1.0079 1.9241 3.8483 7.4217 11.453 14.477 18.783 25.38 36.833 49.02 49.02 1.0694e+118 2.0208e-17 132.96 1 135.29 0.99757 17.778 0.84535 2.7571e-33 2.7571e-33 2.7571e-33 0.00030042 0.002085 0.0055743 0.010412 0.028097 0.074985 0.17067 0.27603 0.34916 0.44348 0.56906 0.74375 1 1 cl_mt_coh 9.0124 11.762 1.3051 -19.968 50.122 0 -6.6378e-18 1 -2.5222e-35 0.0073758 0.0544 0.0992 0.1216 0.21506 0.54409 0.84686 1.1496 2.0842 4.0468 7.6693 11.804 14.888 19.28 26.009 37.692 50.099 50.122 1.6582e+118 3.0718e-17 135.29 0.99757 138.34 1 17.671 0.82377 8.3611e-06 0.0003041 0.0010396 0.0016892 0.0033396 0.0072701 0.013424 0.032153 0.080227 0.17498 0.28188 0.35472 0.44849 0.57315 0.7463 0.95898 1 cl_mt_pol_c 1.4545 1.8239 1.2539 1e-300 18.114 0 0 1 1.6665e-34 0.0035502 0.020561 0.040284 0.15463 0.24919 0.30295 0.36324 0.42573 0.61553 0.91548 1.3809 1.7951 2.1171 2.5757 3.1762 4.3069 9.7256 18.114 6.8818e+116 8.2212e-18 17.778 0.84535 17.671 0.82377 3.3265 1 5.3287e-05 0.00037213 0.0050208 0.0084879 0.019577 0.028388 0.043102 0.078258 0.15639 0.27125 0.3802 0.44714 0.5273 0.62551 0.7517 0.93794 1