Test Shock Persistence and Variance Default Low Cmin (Save + Borr Distribution)

back to Fan's Dynamic Assets Repository Table of Content.

@seealso

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