Risky + Safe Asset Dyna Prog Two-Step Solution (Loop)

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

Contents

function result_map = ff_wkz_vf(varargin)

FF_WKZ_VF solve infinite horizon exo shock + endo asset problem

This program solves the infinite horizon dynamic savings and risky capital asset problem with some ar1 shock. This is the two step solution version of ff_akz_vf. See ff_wkz_evf for details about the second stage.

@param param_map container parameter container

@param support_map container support container

@param armt_map container container with states, choices and shocks grids that are inputs for grid based solution algorithm

@param func_map container container with function handles for consumption cash-on-hand etc.

@return result_map container contains policy function matrix, value function matrix, iteration results, and policy function, value function and iteration results tables.

keys included in result_map:

@example

@include

@seealso

Default

it_param_set = 2;
bl_input_override = true;
[param_map, support_map] = ffs_akz_set_default_param(it_param_set);

% Note: param_map and support_map can be adjusted here or outside to override defaults
% param_map('it_w_n') = 50;
% param_map('it_z_n') = 15;

% get armt and func map
[armt_map, func_map] = ffs_akz_get_funcgrid(param_map, support_map, bl_input_override); % 1 for override
default_params = {param_map support_map armt_map func_map};

Parse Parameters 1

% if varargin only has param_map and support_map,
params_len = length(varargin);
[default_params{1:params_len}] = varargin{:};
param_map = [param_map; default_params{1}];
support_map = [support_map; default_params{2}];
if params_len >= 1 && params_len <= 2
    % If override param_map, re-generate armt and func if they are not
    % provided
    bl_input_override = true;
    [armt_map, func_map] = ffs_akz_get_funcgrid(param_map, support_map, bl_input_override);
else
    % Override all
    armt_map = [armt_map; default_params{3}];
    func_map = [func_map; default_params{4}];
end

% append function name
st_func_name = 'ff_wkz_vf';
support_map('st_profile_name_main') = [st_func_name support_map('st_profile_name_main')];
support_map('st_mat_name_main') = [st_func_name support_map('st_mat_name_main')];
support_map('st_img_name_main') = [st_func_name support_map('st_img_name_main')];

Parse Parameters 2

% armt_map
params_group = values(armt_map, {'ar_w', 'ar_z'});
[ ar_w, ar_z] = params_group{:};
params_group = values(armt_map, {'ar_a_meshk', 'ar_k_mesha', 'mt_coh_wkb', 'it_ameshk_n'});
[ar_a_meshk, ar_k_mesha, mt_coh_wkb, it_ameshk_n] = params_group{:};
% func_map
params_group = values(func_map, {'f_util_log', 'f_util_crra', 'f_cons', 'f_coh'});
[f_util_log, f_util_crra, f_cons, f_coh] = params_group{:};
% param_map
params_group = values(param_map, {'fl_r_save', 'fl_r_borr', 'fl_w',...
    'it_z_n', 'fl_crra', 'fl_beta', 'fl_c_min'});
[fl_r_save, fl_r_borr, fl_wage, it_z_n, fl_crra, fl_beta, fl_c_min] = params_group{:};
params_group = values(param_map, {'it_maxiter_val', 'fl_tol_val', 'fl_tol_pol', 'it_tol_pol_nochange'});
[it_maxiter_val, fl_tol_val, fl_tol_pol, it_tol_pol_nochange] = params_group{:};
% support_map
params_group = values(support_map, {'bl_profile', 'st_profile_path', ...
    'st_profile_prefix', 'st_profile_name_main', 'st_profile_suffix',...
    'bl_time', 'bl_display', 'it_display_every', 'bl_post'});
[bl_profile, st_profile_path, ...
    st_profile_prefix, st_profile_name_main, st_profile_suffix, ...
    bl_time, bl_display, it_display_every, bl_post] = params_group{:};

Initialize Output Matrixes

mt_val_cur = zeros(length(ar_a_meshk),length(ar_z));
mt_val = mt_val_cur - 1;
mt_pol_a = zeros(length(ar_a_meshk),length(ar_z));
mt_pol_a_cur = mt_pol_a - 1;
mt_pol_k = zeros(length(ar_a_meshk),length(ar_z));
mt_pol_k_cur = mt_pol_k - 1;

Initialize Convergence Conditions

bl_vfi_continue = true;
it_iter = 0;
ar_val_diff_norm = zeros([it_maxiter_val, 1]);
ar_pol_diff_norm = zeros([it_maxiter_val, 1]);
mt_pol_perc_change = zeros([it_maxiter_val, it_z_n]);

Iterate Value Function

Loop solution with 4 nested loops

  1. loop 1: over exogenous states
  2. loop 2: over endogenous states
  3. loop 3: over choices
  4. loop 4: add future utility, integration--loop over future shocks
% Start Profile
if (bl_profile)
    close all;
    profile off;
    profile on;
end

% Start Timer
if (bl_time)
    tic;
end

% Value Function Iteration
while bl_vfi_continue
    it_iter = it_iter + 1;

Solve Second Stage Problem k*(w,z)

This is the key difference between this function and ffs_akz_set_functions which solves the two stages jointly

    bl_input_override = true;
    [mt_ev_condi_z_max, ~, mt_ev_condi_z_max_kp, mt_ev_condi_z_max_bp] = ...
        ff_wkz_evf(mt_val_cur, param_map, support_map, armt_map, bl_input_override);

Solve First Stage Problem w*(z) given k*(w,z)

loop 1: over exogenous states

    for it_z_i = 1:length(ar_z)

        % Get 2nd Stage Arrays
        ar_ev_condi_z_max_z = mt_ev_condi_z_max(:, it_z_i);
        ar_w_kstar_z = mt_ev_condi_z_max_kp(:, it_z_i);
        ar_w_astar_z = mt_ev_condi_z_max_bp(:, it_z_i);

        % loop 2: over endogenous states
        for it_coh_j = 1:length(ar_a_meshk)
            % Get cash-on-hand which include k,b,z
            fl_coh = mt_coh_wkb(it_coh_j, it_z_i);

            % loop 3: over choices, only w vector
            % we choose w(z), know from ff_wkz_evf k*(w,z), b*=w-k*
            ar_val_cur = zeros(size(ar_w));
            for it_cohp_k = 1:length(ar_w)
                fl_w_kstar_z = ar_w_kstar_z(it_cohp_k);
                fl_w_astar_z = ar_w_astar_z(it_cohp_k);

                % consumption
                fl_c = f_cons(fl_coh, fl_w_astar_z, fl_w_kstar_z);

                % current utility
                if (fl_crra == 1)
                    ar_val_cur(it_cohp_k) = f_util_log(fl_c);
                    fl_u_neg_c = f_util_log(fl_c_min);
                else
                    ar_val_cur(it_cohp_k) = f_util_crra(fl_c);
                    fl_u_neg_c = f_util_crra(fl_c_min);
                end

                % loop 4: add future utility, integration already done in
                % ff_wkz_evf
                fl_ev_condi_z_max_z = ar_ev_condi_z_max_z(it_cohp_k);
                ar_val_cur(it_cohp_k) = ar_val_cur(it_cohp_k) + fl_beta*fl_ev_condi_z_max_z;

                % Replace if negative consumption
                if fl_c <= 0
                    ar_val_cur(it_cohp_k) = fl_u_neg_c;
                end

            end

            % maximization over loop 3 choices for loop 1+2 states
            it_max_lin_idx = find(ar_val_cur == max(ar_val_cur));
            mt_val(it_coh_j,it_z_i) = ar_val_cur(it_max_lin_idx(1));
            mt_pol_a(it_coh_j,it_z_i) = ar_w_astar_z(it_max_lin_idx(1));
            mt_pol_k(it_coh_j,it_z_i) = ar_w_kstar_z(it_max_lin_idx(1));

        end
    end

Check Tolerance and Continuation

    % Difference across iterations
    ar_val_diff_norm(it_iter) = norm(mt_val - mt_val_cur);
    ar_pol_diff_norm(it_iter) = norm(mt_pol_a - mt_pol_a_cur) + norm(mt_pol_k - mt_pol_k_cur);
    ar_pol_a_perc_change = sum((mt_pol_a ~= mt_pol_a_cur))/(it_ameshk_n);
    ar_pol_k_perc_change = sum((mt_pol_k ~= mt_pol_k_cur))/(it_ameshk_n);
    mt_pol_perc_change(it_iter, :) = mean([ar_pol_a_perc_change;ar_pol_k_perc_change]);

    % Update
    mt_val_cur = mt_val;
    mt_pol_a_cur = mt_pol_a;
    mt_pol_k_cur = mt_pol_k;

    % Print Iteration Results
    if (bl_display && (rem(it_iter, it_display_every)==0))
        fprintf('VAL it_iter:%d, fl_diff:%d, fl_diff_pol:%d\n', ...
            it_iter, ar_val_diff_norm(it_iter), ar_pol_diff_norm(it_iter));
        tb_valpol_iter = array2table([mean(mt_val_cur,1);...
                                      mean(mt_pol_a_cur,1); ...
                                      mean(mt_pol_k_cur,1); ...
                                      mt_val_cur(it_ameshk_n,:); ...
                                      mt_pol_a_cur(it_ameshk_n,:); ...
                                      mt_pol_k_cur(it_ameshk_n,:)]);
        tb_valpol_iter.Properties.VariableNames = strcat('z', string((1:size(mt_val_cur,2))));
        tb_valpol_iter.Properties.RowNames = {'mval', 'map', 'mak', 'Hval', 'Hap', 'Hak'};
        disp('mval = mean(mt_val_cur,1), average value over a')
        disp('map  = mean(mt_pol_a_cur,1), average choice over a')
        disp('mkp  = mean(mt_pol_k_cur,1), average choice over k')
        disp('Hval = mt_val_cur(it_ameshk_n,:), highest a state val')
        disp('Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice')
        disp('mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice')
        disp(tb_valpol_iter);
    end

    % Continuation Conditions:
    % 1. if value function convergence criteria reached
    % 2. if policy function variation over iterations is less than
    % threshold
    if (it_iter == (it_maxiter_val + 1))
        bl_vfi_continue = false;
    elseif ((it_iter == it_maxiter_val) || ...
            (ar_val_diff_norm(it_iter) < fl_tol_val) || ...
            (sum(ar_pol_diff_norm(max(1, it_iter-it_tol_pol_nochange):it_iter)) < fl_tol_pol))
        % Fix to max, run again to save results if needed
        it_iter_last = it_iter;
        it_iter = it_maxiter_val;
    end
VAL it_iter:5, fl_diff:1.134982e+02, fl_diff_pol:2.663606e+02
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    5.4955     5.516    5.5397    5.5666    5.5964    5.6295    5.6659    5.7058    5.7496    5.7973    5.8491    5.9048    5.9642    6.0264     6.089
    map     23.477    23.532    22.582    22.651    21.717    21.778    20.882     19.99    18.542    17.459     15.85    14.164    12.623    10.383    8.3802
    mak     2.0352    2.0352    3.0476    3.0492    4.0576     4.076     5.062    6.0584    7.6727    8.9964    10.872     12.74    14.564     17.14    19.573
    Hval     6.173    6.1848    6.1984     6.214    6.2314    6.2509    6.2731    6.2981    6.3256    6.3565    6.3908    6.4286    6.4698    6.5138    6.5591
    Hap     34.694    34.694    33.673    33.673    32.653    32.653    32.653    31.633    29.592    29.592    27.551    26.531     24.49    22.449    20.408
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    8.1633    9.1837    11.224    13.265    15.306    18.367    21.429

VAL it_iter:10, fl_diff:6.358055e+01, fl_diff_pol:8.816339e+01
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    8.0238    8.0756    8.1348    8.2008    8.2724    8.3503    8.4348    8.5261    8.6246    8.7308    8.8445    8.9656    9.0934    9.2255    9.3568
    map     26.437    26.511    25.581    25.681    24.778     24.89    24.038    23.198     21.95    20.736    19.181    17.616    16.083    13.841    11.818
    mak     2.0352    2.0368    3.0492    3.0492    4.0592    4.0784    5.0676    6.0696    7.5358    9.0324    10.911     12.76    14.724    17.417        20
    Hval    9.2934    9.3253    9.3617    9.4022    9.4488    9.4997    9.5548    9.6158    9.6835    9.7563    9.8376    9.9245    10.018    10.116    10.216
    Hap     38.776    38.776    37.755    37.755    37.755    37.755    36.735    36.735    35.714    34.694    32.653    31.633    29.592    27.551     25.51
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    21.429

VAL it_iter:15, fl_diff:3.879910e+01, fl_diff_pol:4.435572e+01
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    9.4982    9.5705    9.6528     9.744    9.8428    9.9498    10.065    10.189    10.323    10.467    10.621    10.784    10.957    11.136    11.313
    map     27.298    27.378    26.461    26.567    25.693    25.711    24.998    24.186    22.948    21.774     20.25    18.607    17.216    15.022    13.011
    mak     2.0368    2.0368    3.0492    3.0508      4.06    4.1865    5.0692    6.0712    7.5558    9.0388    10.901    12.888    14.761    17.486    20.117
    Hval    11.124    11.171    11.227     11.29    11.358    11.431    11.512    11.602    11.696    11.802    11.913    12.037    12.169    12.308     12.45
    Hap     39.796    39.796    39.796    39.796    38.776    37.755    38.776    37.755    35.714    35.714    34.694    32.653    31.633    29.592    27.551
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816     5.102     5.102    6.1224    8.1633    9.1837    11.224    13.265    15.306    18.367    21.429

VAL it_iter:20, fl_diff:2.470372e+01, fl_diff_pol:3.533853e+01
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    10.434    10.518    10.613    10.718    10.832    10.955    11.088    11.232    11.386    11.552     11.73    11.919     12.12    12.328    12.536
    map     27.672    27.766     26.85     26.97    26.097    26.146    25.432    24.635    23.381    22.238    20.682    19.135    17.751    15.559    13.539
    mak     2.0368    2.0368    3.0492    3.0516      4.06    4.1689      5.07     6.072     7.579    9.0396    10.964    12.883    14.787    17.512    20.191
    Hval    12.275    12.333    12.398     12.47    12.551    12.641    12.736    12.841    12.956     13.08    13.216    13.365    13.523     13.69    13.861
    Hap     40.816    40.816    39.796    39.796    39.796    39.796    38.776    38.776    36.735    36.735    34.694    33.673    32.653    30.612    28.571
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    8.1633    9.1837    11.224    13.265    15.306    18.367    21.429

VAL it_iter:25, fl_diff:1.620750e+01, fl_diff_pol:2.466611e+01
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    11.053    11.142    11.244    11.357    11.479    11.611    11.754    11.908    12.074    12.253    12.446    12.651    12.869    13.095    13.321
    map     27.854    27.946    27.034    27.152     26.28    26.247    25.626    24.855    23.748    22.491    20.936    19.439    18.003    15.862    13.824
    mak     2.0368    2.0368      3.05    3.0516      4.06    4.2569      5.07    6.0728    7.4686    9.0404    10.942    12.875    14.792    17.524    20.241
    Hval    13.008    13.071    13.141    13.224    13.314     13.41    13.516    13.635     13.76    13.902    14.052    14.216    14.392    14.578    14.767
    Hap     40.816    40.816    39.796    40.816    39.796    39.796    39.796    38.776    37.755    36.735    35.714    33.673    32.653    30.612    28.571
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    21.429

VAL it_iter:30, fl_diff:1.081928e+01, fl_diff_pol:1.970702e+01
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    11.471    11.564    11.669    11.786    11.912     12.05    12.198    12.358    12.531    12.718    12.918    13.133    13.361    13.598    13.836
    map     27.945    28.039    27.132    27.252    26.388    26.363    25.753    24.997    23.808    22.613    21.038    19.594    18.163    16.003    13.902
    mak     2.0368    2.0368    3.0508    3.0524    4.0616    4.2569    5.0708    6.0728    7.5326    9.0404    10.983    12.864    14.795    17.528    20.309
    Hval    13.492    13.557    13.636    13.723    13.818    13.919    14.035    14.159    14.293    14.441    14.601    14.773     14.96    15.158    15.357
    Hap     40.816    40.816    40.816    40.816    39.796    38.776    39.796    38.776    37.755    36.735    35.714    34.694    32.653    30.612    28.571
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816     5.102     5.102    6.1224    8.1633    9.1837    11.224    13.265    15.306    18.367    21.429

VAL it_iter:35, fl_diff:7.372941e+00, fl_diff_pol:1.684030e+01
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval     11.76    11.854    11.961     12.08    12.209    12.349      12.5    12.663     12.84     13.03    13.236    13.456     13.69    13.934    14.178
    map     27.988    28.084    27.177    27.305    26.445    26.536    25.825    25.065    23.865    22.685    21.132    19.657    18.271    16.098    13.962
    mak     2.0368    2.0376    3.0516    3.0524    4.0616    4.1441    5.0708    6.0728    7.5414    9.0412    10.984    12.877    14.796    17.536    20.349
    Hval    13.811    13.879    13.962    14.052    14.149    14.254    14.375    14.502    14.643    14.795    14.963    15.143    15.336    15.541    15.746
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    36.735    35.714    34.694    33.673    31.633    28.571
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    21.429

VAL it_iter:40, fl_diff:5.138657e+00, fl_diff_pol:4.699542e+00
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    11.964    12.058    12.166    12.286    12.416    12.557     12.71    12.875    13.053    13.246    13.454    13.677    13.914    14.161    14.409
    map     28.008    28.102      27.2    27.325    26.473    26.561    25.857    25.093    23.852    22.719    21.184    19.708    18.329    16.151    13.987
    mak     2.0368    2.0376    3.0516    3.0524    4.0616    4.1489    5.0708    6.0728    7.5774    9.0412    10.984    12.864    14.796    17.537    20.384
    Hval    14.029    14.099    14.183    14.275    14.373    14.481    14.603    14.732    14.876     15.03    15.203    15.388    15.586    15.796    16.003
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    36.735    35.714    34.694    33.673    31.633    28.571
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    21.429

VAL it_iter:45, fl_diff:3.641932e+00, fl_diff_pol:1.330242e+01
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    12.109    12.205    12.313    12.433    12.564    12.705    12.859    13.025    13.204    13.397    13.606    13.831     14.07    14.319    14.569
    map     28.015    28.111    27.211    27.338    26.483    26.575    25.871    25.107    23.859    22.737    21.202    19.727    18.352    16.182    13.982
    mak     2.0368    2.0376    3.0516    3.0524    4.0616    4.1489    5.0708    6.0728    7.5894    9.0412    10.984    12.864    14.796    17.539    20.423
    Hval    14.181    14.253    14.337     14.43    14.529    14.638    14.761     14.89    15.037    15.191    15.366    15.554    15.755    15.967    16.176
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    36.735    35.714    34.694    33.673    31.633    28.571
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    21.429

VAL it_iter:50, fl_diff:2.611019e+00, fl_diff_pol:1.428878e+01
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    12.215     12.31    12.419    12.539     12.67    12.812    12.965    13.132    13.311    13.505    13.715     13.94     14.18     14.43     14.68
    map     28.019    28.117    27.212    27.344    26.491    26.585    25.879    25.112    23.826    22.744    21.211    19.735    18.364    16.196        14
    mak     2.0368    2.0376    3.0516    3.0524    4.0616    4.1489    5.0708    6.0728    7.6287    9.0412    10.984    12.864    14.796     17.54    20.423
    Hval     14.29    14.363    14.447     14.54    14.639    14.749    14.872    15.002    15.149    15.304     15.48    15.669    15.871    16.085    16.294
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    36.735    35.714    34.694    33.673    31.633    28.571
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    21.429

VAL it_iter:55, fl_diff:1.886585e+00, fl_diff_pol:1.345811e+01
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    12.291    12.386    12.495    12.616    12.747    12.889    13.043    13.209    13.389    13.583    13.793    14.018    14.258    14.509     14.76
    map     28.021    28.118    27.216    27.347    26.496     26.59    25.881    25.115    23.831     22.75    21.216    19.739    18.371    16.206    13.978
    mak     2.0368    2.0376    3.0516    3.0524    4.0616    4.1481    5.0708    6.0728    7.6271    9.0412    10.984    12.864    14.796     17.54    20.457
    Hval    14.368    14.441    14.526    14.618    14.718    14.828    14.952    15.081    15.229    15.384     15.56     15.75    15.953    16.168    16.377
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    37.755    35.714    34.694    33.673    31.633    28.571
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    21.429

VAL it_iter:60, fl_diff:1.370342e+00, fl_diff_pol:1.767399e+00
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    12.347    12.442    12.551    12.672    12.803    12.945    13.099    13.265    13.445    13.639    13.849    14.075    14.315    14.566    14.817
    map     28.022    28.119    27.216    27.348    26.496    26.593    25.882    25.115    23.835    22.752    21.219    19.741    18.372    16.208    14.047
    mak     2.0376    2.0376    3.0516    3.0524    4.0616    4.1481    5.0708    6.0728    7.6271    9.0412    10.984    12.864    14.796     17.54     20.39
    Hval    14.425    14.498    14.582    14.675    14.775    14.885    15.009    15.138    15.286    15.442    15.618    15.808    16.012    16.226    16.436
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    37.755    35.714    34.694    33.673    31.633    29.592
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    20.408

VAL it_iter:65, fl_diff:9.988929e-01, fl_diff_pol:1.020408e+00
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    12.387    12.483    12.592    12.712    12.843    12.985    13.139    13.306    13.486     13.68     13.89    14.116    14.356    14.607    14.859
    map     28.023    28.119    27.216    27.348    26.496    26.595    25.882    25.115    23.835    22.753    21.221    19.741    18.372     16.21    14.047
    mak     2.0376    2.0376    3.0516    3.0524    4.0616    4.1481    5.0708    6.0728    7.6271    9.0412    10.984    12.864    14.796     17.54     20.39
    Hval    14.466    14.539    14.624    14.716    14.816    14.927     15.05     15.18    15.328    15.483     15.66     15.85    16.054    16.269    16.478
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    37.755    35.714    34.694    33.673    31.633    29.592
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    20.408

VAL it_iter:70, fl_diff:7.298583e-01, fl_diff_pol:1.020408e+00
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    12.417    12.513    12.621    12.742    12.873    13.015    13.169    13.336    13.516     13.71     13.92    14.146    14.386    14.637    14.889
    map     28.024    28.119    27.216    27.349    26.496    26.595    25.882    25.115    23.835    22.754    21.221    19.742    18.372    16.212     14.05
    mak     2.0376    2.0376    3.0516    3.0524    4.0616    4.1481    5.0708    6.0728    7.6271    9.0412    10.984    12.864    14.796     17.54     20.39
    Hval    14.496    14.569    14.654    14.746    14.846    14.957     15.08     15.21    15.358    15.514     15.69     15.88    16.084    16.299    16.509
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    37.755    35.714    34.694    33.673    31.633    29.592
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    20.408

VAL it_iter:75, fl_diff:5.341101e-01, fl_diff_pol:0
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    12.439    12.534    12.643    12.764    12.895    13.037    13.191    13.357    13.537    13.732    13.942    14.168    14.408    14.659    14.911
    map     28.024    28.119    27.216    27.349    26.496    26.595    25.882    25.115    23.835    22.755    21.221    19.742    18.372    16.212     14.05
    mak     2.0376    2.0376    3.0516    3.0524    4.0616    4.1481    5.0708    6.0728    7.6271    9.0412    10.984    12.864    14.796     17.54     20.39
    Hval    14.518    14.591    14.676    14.768    14.868    14.979    15.102    15.232     15.38    15.536    15.712    15.902    16.106    16.321    16.531
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    37.755    35.714    34.694    33.673    31.633    29.592
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    20.408

VAL it_iter:80, fl_diff:3.912520e-01, fl_diff_pol:0
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    12.455     12.55    12.659     12.78    12.911    13.053    13.207    13.373    13.553    13.748    13.958    14.184    14.424    14.675    14.927
    map     28.024    28.119    27.216    27.349    26.496    26.595    25.882    25.115    23.835    22.755    21.221    19.742    18.372    16.212     14.05
    mak     2.0376    2.0376    3.0516    3.0524    4.0616    4.1481    5.0708    6.0728    7.6271    9.0412    10.984    12.864    14.796     17.54     20.39
    Hval    14.534    14.607    14.692    14.784    14.884    14.995    15.118    15.248    15.396    15.552    15.728    15.918    16.122    16.337    16.547
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    37.755    35.714    34.694    33.673    31.633    29.592
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    20.408

VAL it_iter:85, fl_diff:2.867881e-01, fl_diff_pol:0
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    12.467    12.562    12.671    12.792    12.923    13.065    13.219    13.385    13.565     13.76     13.97    14.195    14.436    14.687    14.939
    map     28.024    28.119    27.216    27.349    26.496    26.595    25.882    25.115    23.835    22.755    21.221    19.742    18.372    16.212     14.05
    mak     2.0376    2.0376    3.0516    3.0524    4.0616    4.1481    5.0708    6.0728    7.6271    9.0412    10.984    12.864    14.796     17.54     20.39
    Hval    14.546    14.619    14.703    14.796    14.896    15.007     15.13     15.26    15.408    15.563     15.74     15.93    16.134    16.349    16.559
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    37.755    35.714    34.694    33.673    31.633    29.592
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    20.408

VAL it_iter:90, fl_diff:2.103033e-01, fl_diff_pol:0
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    12.475    12.571    12.679      12.8    12.931    13.073    13.227    13.394    13.574    13.768    13.979    14.204    14.445    14.695    14.947
    map     28.024    28.119    27.216    27.349    26.496    26.595    25.882    25.115    23.835    22.755    21.221    19.742    18.372    16.212     14.05
    mak     2.0376    2.0376    3.0516    3.0524    4.0616    4.1481    5.0708    6.0728    7.6271    9.0412    10.984    12.864    14.796     17.54     20.39
    Hval    14.554    14.627    14.712    14.805    14.904    15.015    15.139    15.268    15.416    15.572    15.749    15.939    16.143    16.358    16.568
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    37.755    35.714    34.694    33.673    31.633    29.592
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    20.408

VAL it_iter:95, fl_diff:1.542585e-01, fl_diff_pol:0
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    12.482    12.577    12.686    12.807    12.938     13.08    13.234      13.4     13.58    13.775    13.985     14.21    14.451    14.702    14.953
    map     28.024    28.119    27.216    27.349    26.496    26.595    25.882    25.115    23.835    22.755    21.221    19.742    18.372    16.212     14.05
    mak     2.0376    2.0376    3.0516    3.0524    4.0616    4.1481    5.0708    6.0728    7.6271    9.0412    10.984    12.864    14.796     17.54     20.39
    Hval    14.561    14.634    14.718    14.811    14.911    15.021    15.145    15.275    15.423    15.578    15.755    15.945    16.149    16.364    16.574
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    37.755    35.714    34.694    33.673    31.633    29.592
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    20.408

VAL it_iter:100, fl_diff:1.131695e-01, fl_diff_pol:0
mval = mean(mt_val_cur,1), average value over a
map  = mean(mt_pol_a_cur,1), average choice over a
mkp  = mean(mt_pol_k_cur,1), average choice over k
Hval = mt_val_cur(it_ameshk_n,:), highest a state val
Hap = mt_pol_a_cur(it_ameshk_n,:), highest a state choice
mak = mt_pol_k_cur(it_ameshk_n,:), highest k state choice
              z1        z2        z3        z4        z5        z6        z7        z8        z9       z10       z11       z12       z13       z14       z15  
            ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______    ______

    mval    12.486    12.582     12.69    12.811    12.942    13.084    13.238    13.405    13.585    13.779     13.99    14.215    14.456    14.706    14.958
    map     28.024    28.119    27.216    27.349    26.496    26.595    25.882    25.115    23.835    22.755    21.221    19.742    18.372    16.212     14.05
    mak     2.0376    2.0376    3.0516    3.0524    4.0616    4.1481    5.0708    6.0728    7.6271    9.0412    10.984    12.864    14.796     17.54     20.39
    Hval    14.565    14.638    14.723    14.816    14.915    15.026     15.15    15.279    15.427    15.583     15.76     15.95    16.154    16.369    16.579
    Hap     40.816    41.837    40.816    40.816    39.796    40.816    39.796    38.776    38.776    37.755    35.714    34.694    33.673    31.633    29.592
    Hak     2.0408    2.0408    3.0612    3.0612    4.0816    4.0816     5.102    6.1224    7.1429    9.1837    11.224    13.265    15.306    18.367    20.408

end

% End Timer
if (bl_time)
    toc;
end

% End Profile
if (bl_profile)
    profile off
    profile viewer
    st_file_name = [st_profile_prefix st_profile_name_main st_profile_suffix];
    profsave(profile('info'), strcat(st_profile_path, st_file_name));
end
Elapsed time is 114.535000 seconds.

Process Optimal Choices

result_map = containers.Map('KeyType','char', 'ValueType','any');
result_map('mt_val') = mt_val;

result_map('cl_mt_pol_coh') = {mt_coh_wkb, zeros(1)};
result_map('cl_mt_pol_a') = {mt_pol_a, zeros(1)};
result_map('cl_mt_pol_k') = {mt_pol_k, zeros(1)};
result_map('cl_mt_pol_c') = {f_cons(mt_coh_wkb, mt_pol_a, mt_pol_k), zeros(1)};
result_map('ar_st_pol_names') = ["cl_mt_pol_coh", "cl_mt_pol_a", "cl_mt_pol_k", "cl_mt_pol_c"];

if (bl_post)
    bl_input_override = true;
    result_map('ar_val_diff_norm') = ar_val_diff_norm(1:it_iter_last);
    result_map('ar_pol_diff_norm') = ar_pol_diff_norm(1:it_iter_last);
    result_map('mt_pol_perc_change') = mt_pol_perc_change(1:it_iter_last, :);
    result_map = ff_akz_vf_post(param_map, support_map, armt_map, func_map, result_map, bl_input_override);
end
Warning: Directory already exists. 
Warning: Directory already exists. 
Warning: Directory already exists. 
Warning: Directory already exists. 
valgap = norm(mt_val - mt_val_cur): value function difference across iterations
polgap = norm(mt_pol_a - mt_pol_a_cur): policy function difference across iterations
z1 = z1 perc change: (sum((mt_pol_a ~= mt_pol_a_cur))+sum((mt_pol_k ~= mt_pol_k_cur)))/(2*it_ameshk_n):percentage of state space points conditional on shock where the policy function is changing across iterations
                valgap     polgap        z1            z2            z3            z4            z5            z6            z7            z8            z9           z10           z11           z12           z13           z14           z15    
                _______    ______    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________

    iter=1       227.04    276.59             1             1             1             1             1             1             1             1             1             1             1             1             1             1             1
    iter=2       179.64    2790.4       0.99137       0.99137       0.98549       0.98549       0.97804       0.97804       0.96902       0.95804       0.94588       0.92039       0.88745       0.84824       0.80157          0.72       0.63529
    iter=3       151.16    920.98       0.65176       0.65098       0.52314       0.52275       0.61882       0.62196       0.71686       0.62275       0.89216       0.51608       0.69725       0.54235       0.63294        0.5949       0.62627
    iter=4        130.1    448.85       0.57725       0.57922       0.49608       0.49608       0.55608       0.57137        0.6251       0.52157       0.78431       0.48863       0.63765       0.50039       0.59569       0.50784       0.51647
    iter=5        113.5    266.36       0.53608       0.53569        0.4749       0.47529       0.55882       0.59333       0.58627       0.49137       0.71569       0.47176       0.56196       0.48588       0.51529       0.47686       0.48863
    iter=6       99.938     188.1       0.46471       0.46471       0.42588       0.42549       0.46627           0.5       0.54431       0.44627       0.73922       0.43647       0.54745       0.46471       0.47412       0.44902       0.45412
    iter=7       88.609    140.12       0.37294       0.37569       0.35176       0.35412       0.39608       0.43098       0.46471       0.37333       0.62392       0.37647       0.50431       0.40941       0.41451       0.40314       0.40824
    iter=8           79    120.94       0.27137       0.27412       0.26275        0.2651       0.28627       0.33098       0.32745       0.28588       0.73373       0.29059       0.42902       0.34549       0.33255        0.3298       0.33569
    iter=9       70.745    98.728        0.2051       0.20667       0.20078        0.2051       0.22235       0.26471          0.26       0.21529       0.72157       0.21961       0.36275        0.3102       0.25176       0.25294       0.27059
    iter=10      63.581    88.163       0.16118       0.16431       0.15412       0.15686        0.1698       0.21098       0.20275       0.17725       0.69804       0.17804       0.27686       0.32588       0.20471       0.21725       0.21922
    iter=11      57.334    71.859       0.13059       0.12863       0.12431       0.12431       0.13255       0.16235       0.15529       0.13137       0.61412       0.14353       0.25412       0.26549       0.16431       0.16392       0.17216
    iter=12      51.832    69.137      0.098431       0.10039      0.096863      0.096078       0.10902       0.12314       0.13569       0.12431       0.57569       0.11725       0.20902       0.21451       0.13451       0.14941       0.15333
    iter=13       46.97    54.289      0.086275      0.086275      0.083529      0.086275      0.093725       0.12431       0.11059      0.092157       0.43373      0.094902       0.16745       0.18824       0.11882       0.10941       0.12196
    iter=14      42.649     50.95      0.074118      0.074902      0.071765      0.072941      0.080784       0.13922      0.087059       0.07451       0.42902      0.082353       0.12784       0.11647      0.097255      0.099216       0.10039
    iter=15      38.799    44.356      0.056471      0.055686       0.05451      0.055294      0.064314       0.20196      0.073725      0.071373       0.29569      0.073725       0.15098      0.083137      0.072549      0.090588      0.094118
    iter=16      35.368    47.104      0.047451      0.049412      0.045882      0.047843      0.051765       0.23922      0.067451      0.065882       0.32314      0.060784       0.15451      0.094902      0.069804      0.072549      0.086275
    iter=17      32.283     40.36      0.042353      0.043529      0.043137      0.042353      0.047059       0.25137      0.058039      0.052941       0.22431      0.049804       0.12667      0.088627      0.060784      0.061176      0.061569
    iter=18      29.497    32.552      0.040784      0.040392      0.035686      0.041176      0.041569       0.13647       0.04902      0.042353       0.19529      0.042745        0.1302      0.088235      0.058039      0.051765      0.053333
    iter=19      26.977    34.912      0.034118      0.033725      0.038431      0.035686      0.039216       0.11608      0.041961      0.035686       0.21686      0.039608       0.11961      0.057647      0.056471      0.045098       0.04902
    iter=20      24.704    35.339      0.026275      0.031373      0.027451      0.030588       0.03098       0.19216      0.033333      0.027843       0.21529      0.034902      0.074902      0.040392          0.04      0.045098      0.063529
    iter=21      22.656    25.724      0.026275      0.024706      0.024706      0.024706      0.023529        0.1102      0.026275      0.027059       0.12745      0.034902      0.049804       0.04902      0.032157      0.043529      0.081961
    iter=22      20.807    24.188      0.019216      0.018431      0.020392      0.016863      0.021176       0.06902      0.025098      0.023922        0.1098      0.028627      0.083922      0.031765      0.027059      0.038431          0.08
    iter=23       19.13    28.633      0.018824      0.017647      0.016863      0.019608      0.018824       0.10941      0.021569      0.021176       0.15098      0.021569      0.089804      0.038824      0.025098       0.02902      0.073333
    iter=24      17.602    21.949      0.013725      0.016863      0.014902      0.014118      0.016863      0.087843      0.021961      0.018824      0.092549          0.02      0.087451       0.03451      0.021176      0.021961      0.065098
    iter=25      16.207    24.666      0.014902      0.013333      0.013725      0.014118      0.013725       0.10863       0.01451      0.017255       0.11686      0.019216      0.043529      0.031373      0.022745      0.021569      0.043137
    iter=26      14.932    25.232      0.011373      0.010588      0.012941      0.012549      0.016471      0.070588      0.016863      0.016078       0.12353      0.012941      0.036078      0.039216      0.016863      0.018431      0.061961
    iter=27      13.764    21.592       0.01098      0.013333     0.0090196      0.011373      0.012157      0.089412      0.015294      0.014902      0.055686      0.012549          0.04      0.019216      0.016863      0.015294      0.061961
    iter=28      12.694    21.639     0.0090196      0.010588      0.010196      0.010588      0.012941      0.087451      0.012941       0.01451      0.089412      0.011765      0.041176      0.014118      0.015686       0.01451      0.038039
    iter=29      11.714    23.752     0.0086275     0.0082353     0.0082353     0.0070588     0.0098039      0.072157      0.012941      0.011765       0.10549      0.012549      0.033333       0.01451      0.016078      0.011765      0.062745
    iter=30      10.819    19.707     0.0070588     0.0062745     0.0078431      0.007451     0.0070588      0.072941      0.010588      0.012549      0.013333     0.0098039      0.014118      0.011373       0.01451       0.01098      0.038431
    iter=31      10.001    20.629     0.0078431      0.005098     0.0047059     0.0082353     0.0058824      0.065098     0.0078431     0.0094118      0.081569     0.0082353      0.011765      0.011373       0.01098      0.012549      0.037255
    iter=32      9.2538    18.997     0.0039216     0.0058824     0.0058824     0.0035294      0.010196      0.027451      0.010588     0.0082353      0.070196     0.0094118     0.0094118      0.021569      0.013333      0.010196      0.037647
    iter=33      8.5706    18.069     0.0043137     0.0047059     0.0023529     0.0054902     0.0043137      0.044314     0.0062745     0.0054902      0.033725     0.0078431      0.010196     0.0070588      0.011765     0.0090196      0.061569
    iter=34      7.9453    10.089     0.0031373     0.0047059      0.007451     0.0062745     0.0035294     0.0070588     0.0062745     0.0047059      0.019608     0.0035294     0.0086275     0.0066667      0.010196     0.0098039     0.0094118
    iter=35      7.3729     16.84     0.0027451     0.0031373     0.0023529     0.0027451     0.0058824      0.032157     0.0070588     0.0054902       0.05451     0.0066667     0.0062745      0.007451     0.0070588     0.0090196      0.037647
    iter=36      6.8481    13.579     0.0039216     0.0027451     0.0019608     0.0015686     0.0043137      0.027451     0.0043137     0.0039216     0.0031373      0.005098     0.0094118     0.0039216     0.0090196      0.005098      0.034902
    iter=37      6.3662    14.429     0.0015686     0.0027451     0.0035294     0.0039216     0.0019608      0.020392     0.0062745     0.0023529      0.039216     0.0031373     0.0043137      0.020392     0.0047059      0.007451      0.034118
    iter=38       5.923    4.8384     0.0023529     0.0015686             0     0.0019608     0.0019608     0.0031373     0.0031373     0.0035294     0.0023529     0.0035294     0.0043137     0.0023529     0.0054902     0.0043137     0.0035294
    iter=39       5.515    13.871    0.00078431    0.00078431     0.0031373    0.00039216     0.0035294     0.0031373     0.0023529     0.0019608     0.0027451     0.0019608     0.0047059     0.0039216     0.0070588     0.0043137      0.036471
    iter=40      5.1387    4.6995     0.0015686     0.0011765     0.0023529     0.0019608     0.0031373     0.0015686     0.0011765     0.0019608    0.00039216     0.0031373     0.0027451     0.0039216     0.0027451      0.005098     0.0043137
    iter=41      4.7913    14.522    0.00078431     0.0023529     0.0019608     0.0023529     0.0019608     0.0031373     0.0015686    0.00078431          0.04     0.0031373     0.0031373     0.0015686     0.0035294     0.0039216      0.035294
    iter=42      4.4701    11.807     0.0023529     0.0011765    0.00078431    0.00078431    0.00078431     0.0011765     0.0031373     0.0023529      0.026667     0.0015686     0.0011765     0.0011765     0.0023529     0.0027451     0.0047059
    iter=43      4.1728    5.4465    0.00039216    0.00039216    0.00078431     0.0011765     0.0011765     0.0023529    0.00078431     0.0015686     0.0023529    0.00078431     0.0031373     0.0043137     0.0035294     0.0043137     0.0054902
    iter=44      3.8974    12.575    0.00039216    0.00078431     0.0011765     0.0015686    0.00078431    0.00078431     0.0015686    0.00039216     0.0019608     0.0019608     0.0011765     0.0031373     0.0015686     0.0023529      0.029804
    iter=45      3.6419    13.302             0             0    0.00078431    0.00039216    0.00039216             0    0.00039216     0.0015686    0.00039216     0.0011765    0.00039216    0.00078431    0.00078431     0.0031373      0.033333
    iter=46      3.4048    2.4995             0     0.0011765             0     0.0011765             0     0.0023529    0.00039216             0     0.0027451    0.00039216    0.00078431     0.0015686     0.0015686     0.0015686     0.0011765
    iter=47      3.1844    3.2577    0.00039216    0.00039216             0             0    0.00078431     0.0019608     0.0011765    0.00039216             0    0.00039216     0.0015686    0.00078431     0.0011765     0.0019608     0.0039216
    iter=48      2.9795    2.7878     0.0011765    0.00039216             0    0.00039216    0.00078431    0.00078431     0.0011765             0    0.00039216     0.0011765    0.00039216    0.00039216    0.00039216     0.0015686     0.0011765
    iter=49      2.7887    2.2817    0.00039216             0    0.00039216    0.00078431     0.0015686    0.00039216     0.0019608     0.0011765             0    0.00078431     0.0011765    0.00078431     0.0019608     0.0011765     0.0019608
    iter=50       2.611    14.289             0    0.00078431    0.00039216    0.00039216    0.00039216             0    0.00039216    0.00078431      0.038431    0.00078431    0.00039216     0.0011765    0.00078431    0.00078431    0.00078431
    iter=51      2.4454    3.7202    0.00039216             0             0             0    0.00078431    0.00039216    0.00039216    0.00039216    0.00039216    0.00039216    0.00039216             0     0.0015686     0.0011765     0.0019608
    iter=52       2.291    3.3021             0             0    0.00039216    0.00039216    0.00078431     0.0011765             0    0.00078431    0.00039216             0             0    0.00078431             0     0.0019608     0.0011765
    iter=53      2.1469    4.5043    0.00039216             0    0.00039216    0.00039216             0    0.00078431             0             0     0.0019608     0.0015686    0.00078431    0.00039216     0.0011765    0.00039216             0
    iter=54      2.0123    2.2817             0             0             0    0.00078431    0.00039216             0             0             0             0    0.00039216    0.00078431     0.0011765    0.00078431    0.00078431     0.0019608
    iter=55      1.8866    13.458             0    0.00039216    0.00078431             0     0.0011765             0    0.00039216    0.00039216             0    0.00039216    0.00078431             0             0    0.00078431      0.034118
    iter=56      1.7691    18.926    0.00039216             0             0             0             0             0             0             0    0.00078431    0.00078431    0.00039216             0             0    0.00078431      0.067451
    iter=57      1.6593    1.0204             0    0.00039216             0             0             0    0.00039216             0             0    0.00039216             0             0    0.00039216             0             0             0
    iter=58      1.5565    1.0204    0.00039216             0             0    0.00039216             0             0             0             0             0             0             0    0.00078431             0             0             0
    iter=59      1.4603    1.4431             0             0             0             0             0    0.00039216    0.00039216             0    0.00078431             0             0    0.00039216    0.00039216             0             0
    iter=60      1.3703    1.7674    0.00039216    0.00039216             0             0             0    0.00078431             0             0             0             0     0.0011765             0             0             0             0
    iter=61      1.2861    1.4431             0             0             0             0             0    0.00039216    0.00039216             0             0    0.00039216    0.00078431             0             0             0             0
    iter=62      1.2071    1.4431             0             0    0.00039216             0             0             0             0             0             0             0             0             0             0    0.00078431             0
    iter=63      1.1332    1.0204             0             0             0             0             0    0.00039216             0             0             0             0             0             0             0             0             0
    iter=64      1.0639    1.0204    0.00039216             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=65     0.99889    1.0204             0             0             0             0             0             0             0             0             0    0.00039216    0.00039216             0             0             0             0
    iter=66     0.93798    1.0204             0             0             0    0.00039216             0             0             0             0             0    0.00039216             0             0             0    0.00039216             0
    iter=67     0.88086    1.0204             0             0             0             0             0             0             0             0             0             0             0             0             0    0.00039216    0.00039216
    iter=68     0.82728    1.0204             0             0             0             0             0             0             0             0             0             0             0    0.00039216             0             0    0.00039216
    iter=69     0.77702    1.0204    0.00039216             0             0             0             0    0.00039216             0             0             0             0             0             0             0             0             0
    iter=70     0.72986    1.0204             0             0             0             0             0             0             0             0             0             0             0             0             0    0.00039216    0.00039216
    iter=71      0.6856         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=72     0.64407         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=73     0.60508         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=74     0.56848    1.0204             0             0             0             0             0             0             0             0             0    0.00039216             0             0             0             0             0
    iter=75     0.53411         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=76     0.50184         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=77     0.47154         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=78     0.44308         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=79     0.41636         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=80     0.39125         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=81     0.36767         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=82     0.34552         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=83     0.32471         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=84     0.30516         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=85     0.28679         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=86     0.26953         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=87     0.25331         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=88     0.23808         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=89     0.22376         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=90      0.2103         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=91     0.19766         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=92     0.18578         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=93     0.17461         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=94     0.16412         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=95     0.15426         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=96     0.14499         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=97     0.13628         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=98     0.12809         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=99      0.1204         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=100    0.11317         0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0

tb_val: V(a,z) value at each state space point
                                   z1_0_34741    z2_0_40076    z3_0_4623    z4_0_5333    z5_0_61519    z6_0_70966    z10_1_2567    z11_1_4496    z12_1_6723    z13_1_9291    z14_2_2253    z15_2_567
                                   __________    __________    _________    _________    __________    __________    __________    __________    __________    __________    __________    _________

    coh1:k=0,b=0                    -16.679       -16.679       -16.679      -16.679      -16.679       -16.679       -16.679       -16.679       -16.679       -16.679       -16.679       -16.679 
    coh2:k=1.02041,b=0               1.6951        1.9479        2.2271       2.5214       2.8265        3.1427        4.6343        5.0876        5.5573        6.0229        6.4636        6.8568 
    coh3:k=0,b=1.02041               2.2448        2.5822        2.9487       3.3322       3.7278        4.1353        6.0791        6.7832        7.4306        8.0338        8.5908        9.0897 
    coh4:k=2.04082,b=0               2.9914        3.2625        3.5574       3.8626       4.1718        4.4834        5.9294        6.3827        6.8524        7.3179        7.7586        8.1519 
    coh5:k=1.02041,b=1.02041         3.5061        3.8576        4.2358       4.6264       5.0226        5.4221        7.0898        7.5454        8.0444        8.6067        9.1318        9.6083 
    coh6:k=0,b=2.04082               3.4934        3.8679        4.2667       4.6754       5.0871        5.4996        7.1962        7.7648        8.3175        8.8438        9.3405        9.8339 
    coh7:k=3.06122,b=0               4.2237        4.4947        4.7897       5.0949       5.4041        5.7157        7.0339        7.4054        7.7937        8.1869         8.569         8.919 
    coh8:k=2.04082,b=1.02041         4.3374        4.6305        4.9498       5.3357       5.7345        6.1368        7.7599        8.1816        8.6138         9.051        9.5545        10.014 
    coh9:k=1.02041,b=2.04082          4.334        4.6334        4.9765       5.3827       5.7967        6.2118         7.864        8.2883        8.7762         9.281        9.7589        10.197 
    coh10:k=0,b=3.06122              4.3183         4.623        4.9609       5.3839       5.8119        6.2385        7.9162        8.3803        8.8955        9.3883         9.856        10.354 
    coh11:k=4.08163,b=0              4.9498        5.1987        5.4824       5.7892       6.1044        6.4219        7.7053        8.0447        8.3966        8.7558        9.1106        9.4413 
    coh12:k=3.06122,b=1.02041        5.0609        5.3315        5.6389       5.9718       6.3154        6.6635         8.246        8.6521        9.0659        9.4844        9.9108        10.356 
    coh13:k=2.04082,b=2.04082        5.0576        5.3343        5.6481       5.9875       6.3377        6.7307         8.348        8.7568         9.172        9.6484        10.111        10.536 
    coh14:k=1.02041,b=3.06122        5.0423        5.3242        5.6432       5.9879       6.3465        6.7567        8.3992        8.8106        9.2765        9.7536        10.207        10.653 
    coh15:k=0,b=4.08163              5.0203        5.3072        5.6313        5.981       6.3412        6.7626        8.4302        8.8525         9.346        9.8176        10.266         10.75 
    coh16:k=5.10204,b=0              5.5962         5.816        6.0667       6.3398       6.6306        6.9351        8.1927        8.5182        8.8533        9.1951        9.5337        9.8508 
    coh17:k=4.08163,b=1.02041        5.7049        5.9458        6.2199       6.5185       6.8372        7.1718        8.6475        9.0427        9.4435         9.848        10.247        10.655 
    coh18:k=3.06122,b=2.04082        5.7016        5.9486        6.2289       6.5339       6.8591        7.2002        8.7473        9.1454        9.5478        9.9687        10.418        10.832 
    coh19:k=2.04082,b=3.06122        5.6866        5.9387        6.2241       6.5343       6.8646        7.2107        8.7977        9.1984        9.6082        10.072        10.512         10.92 
    coh20:k=1.02041,b=4.08163        5.6651        5.9221        6.2124       6.5275       6.8625        7.2131        8.8281        9.2315        9.6762        10.135        10.571        11.016 
    coh21:k=0,b=5.10204              5.6387        5.9008        6.1961       6.5161       6.8557        7.2107        8.8479         9.254        9.7218        10.178        10.613        11.079 
    coh22:k=6.12245,b=0              6.1694         6.367        6.5911       6.8365       7.0995        7.3793        8.5954        8.9117        9.2353        9.5644        9.8905        10.205 
    coh23:k=5.10204,b=1.02041        6.2756         6.494        6.7409       7.0114       7.3017        7.6112        8.9965        9.3819        9.7719        10.165        10.552        10.922 
    coh24:k=4.08163,b=2.04082        6.2725        6.4967        6.7498       7.0266       7.3233        7.6391        9.0943        9.4828        9.8745        10.268        10.692        11.096 
    coh25:k=3.06122,b=3.06122        6.2577         6.487        6.7451        7.027       7.3287        7.6494        9.1437        9.5348        9.9281        10.355        10.785        11.181 
    coh26:k=2.04082,b=4.08163        6.2367        6.4707        6.7336       7.0203       7.3267        7.6518        9.1737        9.5674        9.9705        10.417        10.842        11.258 
    coh27:k=1.02041,b=5.10204        6.2109        6.4499        6.7176        7.009         7.32        7.6494        9.1931        9.5896        10.015         10.46        10.882        11.321 
    coh28:k=0,b=6.12245               6.181        6.4252        6.6981       6.9944       7.3099        7.6437         9.206        9.6053        10.047        10.491        10.918        11.366 
    coh29:k=7.14286,b=0              6.6824        6.8619        7.0655        7.289       7.5286        7.7849        8.9456        9.2589        9.5824        9.9115        10.238        10.543 
    coh30:k=6.12245,b=1.02041        6.7864        6.9862        7.2122       7.4603       7.7268        8.0122        9.3111        9.6849        10.064        10.447        10.824        11.182 
    coh31:k=5.10204,b=2.04082        6.7833        6.9889        7.2209       7.4752       7.7479        8.0396        9.4061         9.784        10.165        10.548        10.942        11.337 
    coh32:k=4.08163,b=3.06122        6.7689        6.9793        7.2162       7.4756       7.7533        8.0498        9.4546        9.8352        10.218        10.613        11.033         11.42 
    coh33:k=3.06122,b=4.08163        6.7482        6.9634         7.205        7.469       7.7512        8.0521         9.484        9.8673        10.252        10.674         11.09        11.482 
    coh34:k=2.04082,b=5.10204         6.723         6.943        7.1894        7.458       7.7447        8.0498        9.5031        9.8892        10.281        10.716        11.129        11.544 
    coh35:k=1.02041,b=6.12245        6.6938        6.9188        7.1702       7.4436       7.7348        8.0442        9.5158        9.9047        10.313        10.746        11.159        11.589 
    coh36:k=0,b=7.14286              6.6606        6.8912        7.1479       7.4263       7.7223         8.036         9.524        9.9158        10.336         10.77         11.19        11.623 
    coh37:k=8.16327,b=0              7.1458        7.3107        7.4979       7.7039       7.9242        8.1603        9.2863        9.5944        9.9086        10.227        10.543        10.839 
    coh38:k=7.14286,b=1.02041        7.2476        7.4324        7.6415       7.8717       8.1183        8.3831        9.6194        9.9597         10.33        10.703        11.072        11.422 
    coh39:k=6.12245,b=2.04082        7.2445         7.435          7.65       7.8863       8.1391        8.4101        9.6892        10.057        10.429        10.803        11.173         11.56 
    coh40:k=5.10204,b=3.06122        7.2304        7.4256        7.6455       7.8867       8.1443        8.4201        9.7368        10.107        10.481        10.856        11.263        11.642 
    coh41:k=4.08163,b=4.08163        7.2102          7.41        7.6345       7.8802       8.1424        8.4224        9.7657        10.139        10.515        10.912        11.319        11.694 
    coh42:k=3.06122,b=5.10204        7.1855          7.39        7.6191       7.8694       8.1359        8.4201        9.7846        10.161        10.538        10.953        11.358        11.753 
    coh43:k=2.04082,b=6.12245        7.1569        7.3663        7.6003       7.8553       8.1262        8.4146         9.797        10.176         10.56        10.983        11.387        11.797 
    coh44:k=1.02041,b=7.14286        7.1246        7.3393        7.5785       7.8384       8.1139        8.4065        9.8051        10.187        10.583        11.006        11.409        11.831 
    coh45:k=0,b=8.16327              7.0883         7.309        7.5538       7.8189       8.0994        8.3964        9.8098        10.195          10.6        11.024        11.436        11.858 
    coh46:k=9.18367,b=0              7.5679        7.7208        7.8945        8.086       8.2903        8.5093        9.5949        9.8939        10.199        10.508        10.815        11.103 
    coh47:k=8.16327,b=1.02041        7.6674        7.8399        8.0352       8.2504       8.4806        8.7277        9.8991        10.228        10.577        10.941        11.301        11.643 
    coh48:k=7.14286,b=2.04082        7.6645        7.8425        8.0436       8.2648        8.501        8.7542        9.9554        10.312        10.674        11.039          11.4        11.769 
    coh49:k=6.12245,b=3.06122        7.6506        7.8333        8.0392       8.2652       8.5062        8.7641        10.002        10.362        10.726        11.092        11.478         11.85 
    coh50:k=5.10204,b=4.08163        7.6309         7.818        8.0283       8.2588       8.5042        8.7664        10.031        10.393        10.759        11.134        11.533        11.901 
    coh1226:k=50,b=0                 14.953        14.999        15.052       15.111       15.177        15.247        15.595        15.701        15.818        15.945        16.077        16.208 
    coh1227:k=48.9796,b=1.02041      14.975        15.026        15.085        15.15       15.222        15.301        15.694        15.817         15.95        16.092        16.242        16.392 
    coh1228:k=47.9592,b=2.04082      14.975        15.027        15.087       15.154       15.228        15.308        15.717        15.842        15.978        16.124        16.282        16.439 
    coh1229:k=46.9388,b=3.06122      14.971        15.024        15.086       15.154       15.229        15.311        15.729        15.857        15.995        16.149        16.309        16.467 
    coh1230:k=45.9184,b=4.08163      14.967        15.021        15.083       15.153       15.229        15.312        15.737        15.866         16.01        16.165        16.327        16.488 
    coh1231:k=44.898,b=5.10204       14.961        15.016        15.079        15.15       15.227        15.311        15.742        15.873         16.02        16.177         16.34        16.506 
    coh1232:k=43.8776,b=6.12245      14.955        15.011        15.075       15.146       15.224         15.31        15.746        15.879        16.028        16.187        16.352        16.521 
    coh1233:k=42.8571,b=7.14286      14.948        15.005         15.07       15.142       15.221        15.307        15.748        15.884        16.035        16.194        16.363        16.532 
    coh1234:k=41.8367,b=8.16327      14.941        14.998        15.064       15.137       15.217        15.304         15.75        15.887        16.039          16.2        16.372        16.542 
    coh1235:k=40.8163,b=9.18367      14.934        14.992        15.058       15.132       15.213        15.301         15.75        15.889        16.043        16.205        16.379         16.55 
    coh1236:k=39.7959,b=10.2041      14.926        14.984        15.052       15.126       15.208        15.297         15.75        15.891        16.046        16.209        16.386        16.557 
    coh1237:k=38.7755,b=11.2245      14.917        14.977        15.045        15.12       15.203        15.293         15.75        15.892        16.048        16.213        16.391        16.563 
    coh1238:k=37.7551,b=12.2449      14.908        14.969        15.037       15.114       15.198        15.289        15.749        15.892         16.05        16.216        16.395        16.568 
    coh1239:k=36.7347,b=13.2653      14.899         14.96         15.03       15.108       15.192        15.284        15.748        15.892        16.051        16.219        16.399        16.572 
    coh1240:k=35.7143,b=14.2857      14.892        14.951        15.022       15.101       15.186        15.279        15.746        15.891        16.051        16.221        16.402        16.576 
    coh1241:k=34.6939,b=15.3061      14.884        14.943        15.014       15.093        15.18        15.274        15.744         15.89        16.051        16.222        16.404         16.58 
    coh1242:k=33.6735,b=16.3265      14.877        14.937        15.005       15.086       15.173        15.268        15.742        15.889        16.051        16.224        16.407        16.582 
    coh1243:k=32.6531,b=17.3469      14.869        14.929        14.998       15.078       15.166        15.262         15.74        15.887        16.051        16.224        16.408        16.585 
    coh1244:k=31.6327,b=18.3673      14.861        14.922        14.992       15.069       15.159        15.256        15.737        15.885         16.05        16.225         16.41        16.587 
    coh1245:k=30.6122,b=19.3878      14.853        14.914        14.985       15.063       15.152         15.25        15.734        15.882        16.049        16.225        16.411        16.589 
    coh1246:k=29.5918,b=20.4082      14.844        14.906        14.977       15.056       15.144        15.243        15.731         15.88        16.048        16.225        16.412        16.591 
    coh1247:k=28.5714,b=21.4286      14.835        14.898         14.97        15.05       15.136        15.236        15.727        15.877        16.046        16.224        16.413        16.592 
    coh1248:k=27.551,b=22.449        14.826        14.889        14.962       15.043        15.13        15.229        15.724        15.874        16.044        16.223        16.413        16.594 
    coh1249:k=26.5306,b=23.4694      14.816        14.881        14.954       15.036       15.124        15.222         15.72        15.872        16.043        16.222        16.413        16.595 
    coh1250:k=25.5102,b=24.4898      14.806        14.871        14.946       15.028       15.117        15.214        15.716        15.869         16.04        16.221        16.413        16.595 
    coh1251:k=24.4898,b=25.5102      14.796        14.862        14.937        15.02       15.111        15.208        15.712        15.866        16.038         16.22        16.413        16.596 
    coh1252:k=23.4694,b=26.5306      14.787        14.852        14.928       15.012       15.104        15.202        15.707        15.863        16.036        16.218        16.413        16.597 
    coh1253:k=22.449,b=27.551        14.779        14.842        14.919       15.004       15.096        15.196        15.703         15.86        16.033        16.216        16.412        16.597 
    coh1254:k=21.4286,b=28.5714      14.771        14.834         14.91       14.996       15.089        15.189        15.698        15.857         16.03        16.214        16.411        16.597 
    coh1255:k=20.4082,b=29.5918      14.762        14.826          14.9       14.987       15.081        15.183        15.693        15.853        16.027        16.212        16.411        16.597 
    coh1256:k=19.3878,b=30.6122      14.753        14.818        14.892       14.978       15.073        15.176        15.688         15.85        16.024         16.21         16.41        16.597 
    coh1257:k=18.3673,b=31.6327      14.744         14.81        14.885       14.968       15.065        15.169        15.683        15.846        16.021        16.208        16.409        16.597 
    coh1258:k=17.3469,b=32.6531      14.735        14.801        14.877       14.961       15.056        15.161        15.677        15.842        16.017        16.206        16.407        16.597 
    coh1259:k=16.3265,b=33.6735      14.725        14.792        14.869       14.953       15.048        15.154        15.672        15.838        16.014        16.204        16.406        16.597 
    coh1260:k=15.3061,b=34.6939      14.715        14.783         14.86       14.946       15.038        15.146        15.668        15.834         16.01        16.202        16.404        16.596 
    coh1261:k=14.2857,b=35.7143      14.705        14.774        14.852       14.938       15.031        15.138        15.664        15.829        16.006        16.199        16.403        16.596 
    coh1262:k=13.2653,b=36.7347      14.694        14.764        14.843        14.93       15.024         15.13        15.659        15.825        16.002        16.197        16.401        16.595 
    coh1263:k=12.2449,b=37.7551      14.683        14.754        14.834       14.922       15.017        15.121        15.654         15.82        15.998        16.194        16.399        16.594 
    coh1264:k=11.2245,b=38.7755      14.675        14.743        14.824       14.913       15.009        15.112        15.649        15.815        15.994        16.191        16.397        16.593 
    coh1265:k=10.2041,b=39.7959      14.666        14.733        14.814       14.905       15.002        15.106        15.644         15.81        15.991        16.189        16.395        16.592 
    coh1266:k=9.18367,b=40.8163      14.657        14.725        14.804       14.895       14.994        15.099        15.639        15.805        15.987        16.186        16.393        16.591 
    coh1267:k=8.16327,b=41.8367      14.648        14.716        14.794       14.886       14.986        15.092        15.633          15.8        15.984        16.183        16.391         16.59 
    coh1268:k=7.14286,b=42.8571      14.639        14.708        14.786       14.876       14.977        15.084        15.628        15.794         15.98         16.18        16.389        16.589 
    coh1269:k=6.12245,b=43.8776      14.629        14.699        14.778       14.866       14.968        15.077        15.622        15.789        15.976        16.176        16.386        16.588 
    coh1270:k=5.10204,b=44.898       14.619         14.69         14.77       14.857       14.959        15.069        15.616        15.784        15.972        16.173        16.384        16.587 
    coh1271:k=4.08163,b=45.9184      14.609         14.68        14.761        14.85        14.95        15.061         15.61         15.78        15.968         16.17        16.381        16.586 
    coh1272:k=3.06122,b=46.9388      14.598         14.67        14.752       14.842        14.94        15.053        15.604        15.775        15.964        16.166        16.378        16.584 
    coh1273:k=2.04082,b=47.9592      14.587         14.66        14.743       14.834       14.931        15.045        15.597         15.77        15.959        16.162        16.376        16.583 
    coh1274:k=1.02041,b=48.9796      14.576         14.65        14.734       14.825       14.924        15.036        15.591        15.765        15.955        16.158        16.373        16.581 
    coh1275:k=0,b=50                 14.566        14.639        14.724       14.817       14.916        15.027        15.584        15.761         15.95        16.155         16.37        16.579 

tb_pol_a: optimal safe savings choice for each state space point
                                   z1_0_34741    z2_0_40076    z3_0_4623    z4_0_5333    z5_0_61519    z6_0_70966    z10_1_2567    z11_1_4496    z12_1_6723    z13_1_9291    z14_2_2253    z15_2_567
                                   __________    __________    _________    _________    __________    __________    __________    __________    __________    __________    __________    _________

    coh1:k=0,b=0                          0             0            0            0             0             0             0             0             0             0             0            0  
    coh2:k=1.02041,b=0                    0             0            0            0             0             0             0             0             0             0             0            0  
    coh3:k=0,b=1.02041                    0             0            0            0             0             0             0             0             0             0             0            0  
    coh4:k=2.04082,b=0                    0             0            0            0             0             0             0             0             0             0             0            0  
    coh5:k=1.02041,b=1.02041              0             0            0            0             0             0             0             0             0             0             0            0  
    coh6:k=0,b=2.04082                    0             0            0            0             0             0             0             0             0             0             0            0  
    coh7:k=3.06122,b=0                    0             0            0            0             0             0             0             0             0             0             0            0  
    coh8:k=2.04082,b=1.02041              0             0            0            0             0             0             0             0             0             0             0            0  
    coh9:k=1.02041,b=2.04082              0             0            0            0             0             0             0             0             0             0             0            0  
    coh10:k=0,b=3.06122                   0             0            0            0             0             0             0             0             0             0             0            0  
    coh11:k=4.08163,b=0              1.0204        1.0204            0            0             0             0             0             0             0             0             0            0  
    coh12:k=3.06122,b=1.02041        1.0204        1.0204            0            0             0             0             0             0             0             0             0            0  
    coh13:k=2.04082,b=2.04082        1.0204        1.0204            0            0             0             0             0             0             0             0             0            0  
    coh14:k=1.02041,b=3.06122        1.0204        1.0204            0            0             0             0             0             0             0             0             0            0  
    coh15:k=0,b=4.08163              1.0204        1.0204            0            0             0             0             0             0             0             0             0            0  
    coh16:k=5.10204,b=0              2.0408        2.0408       1.0204       1.0204             0             0             0             0             0             0             0            0  
    coh17:k=4.08163,b=1.02041        2.0408        2.0408       1.0204       1.0204             0             0             0             0             0             0             0            0  
    coh18:k=3.06122,b=2.04082        2.0408        2.0408       1.0204       1.0204             0             0             0             0             0             0             0            0  
    coh19:k=2.04082,b=3.06122        2.0408        2.0408       1.0204       1.0204             0             0             0             0             0             0             0            0  
    coh20:k=1.02041,b=4.08163        2.0408        2.0408       1.0204       1.0204             0             0             0             0             0             0             0            0  
    coh21:k=0,b=5.10204              2.0408        2.0408       1.0204       1.0204             0             0             0             0             0             0             0            0  
    coh22:k=6.12245,b=0              3.0612        3.0612       2.0408       2.0408        1.0204             0             0             0             0             0             0            0  
    coh23:k=5.10204,b=1.02041        3.0612        3.0612       2.0408       2.0408        1.0204             0             0             0             0             0             0            0  
    coh24:k=4.08163,b=2.04082        3.0612        3.0612       2.0408       2.0408        1.0204             0             0             0             0             0             0            0  
    coh25:k=3.06122,b=3.06122        3.0612        3.0612       2.0408       2.0408        1.0204             0             0             0             0             0             0            0  
    coh26:k=2.04082,b=4.08163        3.0612        3.0612       2.0408       2.0408        1.0204             0             0             0             0             0             0            0  
    coh27:k=1.02041,b=5.10204        3.0612        3.0612       2.0408       2.0408        1.0204             0             0             0             0             0             0            0  
    coh28:k=0,b=6.12245              3.0612        3.0612       2.0408       2.0408        1.0204             0             0             0             0             0             0            0  
    coh29:k=7.14286,b=0              4.0816        4.0816       3.0612       3.0612        2.0408        1.0204             0             0             0             0             0            0  
    coh30:k=6.12245,b=1.02041        4.0816        4.0816       3.0612       3.0612        2.0408        1.0204             0             0             0             0             0            0  
    coh31:k=5.10204,b=2.04082        4.0816        4.0816       3.0612       3.0612        2.0408        1.0204             0             0             0             0             0            0  
    coh32:k=4.08163,b=3.06122        4.0816        4.0816       3.0612       3.0612        2.0408        1.0204             0             0             0             0             0            0  
    coh33:k=3.06122,b=4.08163        4.0816        4.0816       3.0612       3.0612        2.0408        1.0204             0             0             0             0             0            0  
    coh34:k=2.04082,b=5.10204        4.0816        4.0816       3.0612       3.0612        2.0408        1.0204             0             0             0             0             0            0  
    coh35:k=1.02041,b=6.12245        4.0816        4.0816       3.0612       3.0612        2.0408        1.0204             0             0             0             0             0            0  
    coh36:k=0,b=7.14286              4.0816        4.0816       3.0612       3.0612        2.0408        1.0204             0             0             0             0             0            0  
    coh37:k=8.16327,b=0               5.102         5.102       4.0816       4.0816        3.0612        2.0408             0             0             0             0             0            0  
    coh38:k=7.14286,b=1.02041         5.102         5.102       4.0816       4.0816        3.0612        2.0408             0             0             0             0             0            0  
    coh39:k=6.12245,b=2.04082         5.102         5.102       4.0816       4.0816        3.0612        2.0408             0             0             0             0             0            0  
    coh40:k=5.10204,b=3.06122         5.102         5.102       4.0816       4.0816        3.0612        2.0408             0             0             0             0             0            0  
    coh41:k=4.08163,b=4.08163         5.102         5.102       4.0816       4.0816        3.0612        2.0408             0             0             0             0             0            0  
    coh42:k=3.06122,b=5.10204         5.102         5.102       4.0816       4.0816        3.0612        2.0408             0             0             0             0             0            0  
    coh43:k=2.04082,b=6.12245         5.102         5.102       4.0816       4.0816        3.0612        2.0408             0             0             0             0             0            0  
    coh44:k=1.02041,b=7.14286         5.102         5.102       4.0816       4.0816        3.0612        2.0408             0             0             0             0             0            0  
    coh45:k=0,b=8.16327               5.102         5.102       4.0816       4.0816        3.0612        2.0408             0             0             0             0             0            0  
    coh46:k=9.18367,b=0              6.1224        6.1224        5.102        5.102        4.0816        4.0816             0             0             0             0             0            0  
    coh47:k=8.16327,b=1.02041        6.1224        6.1224        5.102        5.102        4.0816        4.0816             0             0             0             0             0            0  
    coh48:k=7.14286,b=2.04082        6.1224        6.1224        5.102        5.102        4.0816        4.0816             0             0             0             0             0            0  
    coh49:k=6.12245,b=3.06122        6.1224        6.1224        5.102        5.102        4.0816        4.0816             0             0             0             0             0            0  
    coh50:k=5.10204,b=4.08163        6.1224        6.1224        5.102        5.102        4.0816        4.0816             0             0             0             0             0            0  
    coh1226:k=50,b=0                 44.898        44.898       43.878       43.878        42.857        42.857        37.755        35.714        32.653        30.612        27.551        24.49  
    coh1227:k=48.9796,b=1.02041      44.898        44.898       43.878       43.878        42.857        42.857        38.776        36.735        34.694        32.653        29.592       26.531  
    coh1228:k=47.9592,b=2.04082      44.898        44.898       43.878       43.878        42.857        42.857        38.776        36.735        34.694        32.653        30.612       27.551  
    coh1229:k=46.9388,b=3.06122      44.898        44.898       43.878       43.878        42.857        42.857        38.776        36.735        34.694        33.673        30.612       27.551  
    coh1230:k=45.9184,b=4.08163      44.898        44.898       43.878       43.878        42.857        42.857        38.776        36.735        35.714        33.673        30.612       29.592  
    coh1231:k=44.898,b=5.10204       44.898        44.898       43.878       43.878        42.857        42.857        38.776        36.735        35.714        33.673        30.612       29.592  
    coh1232:k=43.8776,b=6.12245      44.898        44.898       43.878       43.878        42.857        42.857        38.776        37.755        35.714        33.673        31.633       29.592  
    coh1233:k=42.8571,b=7.14286      44.898        44.898       43.878       43.878        42.857        42.857        38.776        37.755        35.714        33.673        31.633       29.592  
    coh1234:k=41.8367,b=8.16327      44.898        44.898       43.878       43.878        42.857        42.857        38.776        37.755        35.714        33.673        31.633       29.592  
    coh1235:k=40.8163,b=9.18367      44.898        44.898       43.878       43.878        42.857        42.857        38.776        37.755        35.714        33.673        31.633       29.592  
    coh1236:k=39.7959,b=10.2041      44.898        44.898       43.878       43.878        42.857        42.857        38.776        37.755        35.714        33.673        31.633       29.592  
    coh1237:k=38.7755,b=11.2245      44.898        44.898       43.878       43.878        42.857        42.857        38.776        37.755        35.714        34.694        31.633       29.592  
    coh1238:k=37.7551,b=12.2449      44.898        44.898       43.878       43.878        42.857        42.857        38.776        37.755        35.714        34.694        31.633       29.592  
    coh1239:k=36.7347,b=13.2653      44.898        44.898       43.878       43.878        42.857        42.857        38.776        37.755        35.714        34.694        31.633       29.592  
    coh1240:k=35.7143,b=14.2857      43.878        44.898       43.878       43.878        42.857        42.857        38.776        37.755        35.714        34.694        31.633       29.592  
    coh1241:k=34.6939,b=15.3061      43.878        43.878       43.878       43.878        42.857        42.857        38.776        37.755        35.714        34.694        31.633       29.592  
    coh1242:k=33.6735,b=16.3265      43.878        43.878       43.878       43.878        42.857        42.857        38.776        37.755        35.714        34.694        31.633       29.592  
    coh1243:k=32.6531,b=17.3469      43.878        43.878       42.857       43.878        42.857        42.857        38.776        37.755        35.714        34.694        31.633       29.592  
    coh1244:k=31.6327,b=18.3673      43.878        43.878       42.857       43.878        42.857        42.857        38.776        37.755        35.714        34.694        31.633       29.592  
    coh1245:k=30.6122,b=19.3878      43.878        43.878       42.857       42.857        42.857        42.857        38.776        37.755        35.714        34.694        31.633       29.592  
    coh1246:k=29.5918,b=20.4082      43.878        43.878       42.857       42.857        42.857        42.857        38.776        37.755        35.714        34.694        31.633       29.592  
    coh1247:k=28.5714,b=21.4286      43.878        43.878       42.857       42.857        41.837        42.857        38.776        37.755        35.714        34.694        31.633       29.592  
    coh1248:k=27.551,b=22.449        43.878        43.878       42.857       42.857        41.837        42.857        38.776        36.735        35.714        34.694        31.633       29.592  
    coh1249:k=26.5306,b=23.4694      43.878        43.878       42.857       42.857        41.837        42.857        38.776        36.735        35.714        34.694        31.633       29.592  
    coh1250:k=25.5102,b=24.4898      43.878        43.878       42.857       42.857        41.837        42.857        38.776        36.735        35.714        34.694        31.633       29.592  
    coh1251:k=24.4898,b=25.5102      43.878        43.878       42.857       42.857        41.837        41.837        38.776        36.735        35.714        34.694        31.633       29.592  
    coh1252:k=23.4694,b=26.5306      42.857        43.878       42.857       42.857        41.837        41.837        38.776        36.735        35.714        34.694        31.633       29.592  
    coh1253:k=22.449,b=27.551        42.857        43.878       42.857       42.857        41.837        41.837        38.776        36.735        35.714        34.694        31.633       29.592  
    coh1254:k=21.4286,b=28.5714      42.857        42.857       42.857       42.857        41.837        41.837        38.776        36.735        35.714        34.694        31.633       29.592  
    coh1255:k=20.4082,b=29.5918      42.857        42.857       41.837       42.857        41.837        41.837        38.776        36.735        35.714        34.694        31.633       29.592  
    coh1256:k=19.3878,b=30.6122      42.857        42.857       41.837       42.857        41.837        41.837        38.776        36.735        35.714        33.673        31.633       29.592  
    coh1257:k=18.3673,b=31.6327      42.857        42.857       41.837       42.857        41.837        41.837        38.776        36.735        35.714        33.673        31.633       29.592  
    coh1258:k=17.3469,b=32.6531      42.857        42.857       41.837       41.837        41.837        41.837        38.776        36.735        35.714        33.673        31.633       29.592  
    coh1259:k=16.3265,b=33.6735      42.857        42.857       41.837       41.837        41.837        41.837        37.755        36.735        35.714        33.673        31.633       29.592  
    coh1260:k=15.3061,b=34.6939      42.857        42.857       41.837       41.837        41.837        41.837        37.755        36.735        35.714        33.673        31.633       29.592  
    coh1261:k=14.2857,b=35.7143      42.857        42.857       41.837       41.837        40.816        41.837        37.755        36.735        35.714        33.673        31.633       29.592  
    coh1262:k=13.2653,b=36.7347      42.857        42.857       41.837       41.837        40.816        41.837        37.755        36.735        35.714        33.673        31.633       29.592  
    coh1263:k=12.2449,b=37.7551      42.857        42.857       41.837       41.837        40.816        41.837        37.755        36.735        35.714        33.673        31.633       29.592  
    coh1264:k=11.2245,b=38.7755      41.837        42.857       41.837       41.837        40.816        40.816        37.755        36.735        34.694        33.673        31.633       29.592  
    coh1265:k=10.2041,b=39.7959      41.837        41.837       41.837       41.837        40.816        40.816        37.755        36.735        34.694        33.673        31.633       29.592  
    coh1266:k=9.18367,b=40.8163      41.837        41.837       41.837       41.837        40.816        40.816        37.755        36.735        34.694        33.673        31.633       29.592  
    coh1267:k=8.16327,b=41.8367      41.837        41.837       40.816       41.837        40.816        40.816        37.755        36.735        34.694        33.673        31.633       29.592  
    coh1268:k=7.14286,b=42.8571      41.837        41.837       40.816       41.837        40.816        40.816        37.755        36.735        34.694        33.673        31.633       29.592  
    coh1269:k=6.12245,b=43.8776      41.837        41.837       40.816       41.837        40.816        40.816        37.755        36.735        34.694        33.673        31.633       29.592  
    coh1270:k=5.10204,b=44.898       41.837        41.837       40.816       40.816        40.816        40.816        37.755        35.714        34.694        33.673        31.633       29.592  
    coh1271:k=4.08163,b=45.9184      41.837        41.837       40.816       40.816        40.816        40.816        37.755        35.714        34.694        33.673        31.633       29.592  
    coh1272:k=3.06122,b=46.9388      41.837        41.837       40.816       40.816        40.816        40.816        37.755        35.714        34.694        33.673        31.633       29.592  
    coh1273:k=2.04082,b=47.9592      41.837        41.837       40.816       40.816        39.796        40.816        37.755        35.714        34.694        33.673        31.633       29.592  
    coh1274:k=1.02041,b=48.9796      41.837        41.837       40.816       40.816        39.796        40.816        37.755        35.714        34.694        33.673        31.633       29.592  
    coh1275:k=0,b=50                 40.816        41.837       40.816       40.816        39.796        40.816        37.755        35.714        34.694        33.673        31.633       29.592  

tb_pol_k: optimal risky investment choice for each state space point
                                   z1_0_34741    z2_0_40076    z3_0_4623    z4_0_5333    z5_0_61519    z6_0_70966    z10_1_2567    z11_1_4496    z12_1_6723    z13_1_9291    z14_2_2253    z15_2_567
                                   __________    __________    _________    _________    __________    __________    __________    __________    __________    __________    __________    _________

    coh1:k=0,b=0                          0             0            0            0             0             0             0             0             0             0             0            0  
    coh2:k=1.02041,b=0               1.0204        1.0204       1.0204       1.0204        1.0204        1.0204        1.0204        1.0204        1.0204        1.0204        1.0204       1.0204  
    coh3:k=0,b=1.02041               1.0204        1.0204       1.0204       1.0204        1.0204        1.0204        2.0408        2.0408        2.0408        2.0408        2.0408       2.0408  
    coh4:k=2.04082,b=0               2.0408        2.0408       2.0408       2.0408        2.0408        2.0408        1.0204        1.0204        1.0204        1.0204        1.0204       1.0204  
    coh5:k=1.02041,b=1.02041         2.0408        2.0408       2.0408       2.0408        2.0408        2.0408        2.0408        2.0408        3.0612        3.0612        3.0612       3.0612  
    coh6:k=0,b=2.04082               2.0408        2.0408       2.0408       2.0408        2.0408        2.0408        2.0408        3.0612        3.0612        3.0612        3.0612       4.0816  
    coh7:k=3.06122,b=0               2.0408        2.0408       2.0408       2.0408        2.0408        2.0408        2.0408        2.0408        2.0408        2.0408        2.0408       2.0408  
    coh8:k=2.04082,b=1.02041         2.0408        2.0408       2.0408       3.0612        3.0612        3.0612        3.0612        3.0612        3.0612        3.0612        4.0816       4.0816  
    coh9:k=1.02041,b=2.04082         2.0408        2.0408       3.0612       3.0612        3.0612        3.0612        3.0612        3.0612        4.0816        4.0816        4.0816       4.0816  
    coh10:k=0,b=3.06122              2.0408        2.0408       3.0612       3.0612        3.0612        3.0612        3.0612        4.0816        4.0816        4.0816        4.0816        5.102  
    coh11:k=4.08163,b=0              2.0408        2.0408       3.0612       3.0612        3.0612        3.0612        3.0612        3.0612        3.0612        3.0612        3.0612       3.0612  
    coh12:k=3.06122,b=1.02041        2.0408        2.0408       3.0612       3.0612        3.0612        3.0612        4.0816        4.0816        4.0816        4.0816         5.102        5.102  
    coh13:k=2.04082,b=2.04082        2.0408        2.0408       3.0612       3.0612        3.0612        4.0816        4.0816        4.0816        4.0816         5.102         5.102        5.102  
    coh14:k=1.02041,b=3.06122        2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        4.0816        4.0816         5.102         5.102         5.102       6.1224  
    coh15:k=0,b=4.08163              2.0408        2.0408       3.0612       3.0612        3.0612        4.0816        4.0816         5.102         5.102         5.102         5.102       6.1224  
    coh16:k=5.10204,b=0              2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        4.0816        4.0816        4.0816        4.0816        4.0816       4.0816  
    coh17:k=4.08163,b=1.02041        2.0408        2.0408       3.0612       3.0612        4.0816        4.0816         5.102         5.102         5.102         5.102         5.102       6.1224  
    coh18:k=3.06122,b=2.04082        2.0408        2.0408       3.0612       3.0612        4.0816        4.0816         5.102         5.102         5.102        6.1224        6.1224       6.1224  
    coh19:k=2.04082,b=3.06122        2.0408        2.0408       3.0612       3.0612        4.0816        4.0816         5.102         5.102        6.1224        6.1224        6.1224       7.1429  
    coh20:k=1.02041,b=4.08163        2.0408        2.0408       3.0612       3.0612        4.0816        4.0816         5.102         5.102        6.1224        6.1224        6.1224       7.1429  
    coh21:k=0,b=5.10204              2.0408        2.0408       3.0612       3.0612        4.0816        4.0816         5.102         5.102        6.1224        6.1224        7.1429       7.1429  
    coh22:k=6.12245,b=0              2.0408        2.0408       3.0612       3.0612        4.0816         5.102         5.102         5.102         5.102         5.102         5.102       4.0816  
    coh23:k=5.10204,b=1.02041        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        6.1224        6.1224        6.1224        6.1224        6.1224       7.1429  
    coh24:k=4.08163,b=2.04082        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        6.1224        6.1224        6.1224        6.1224        7.1429       7.1429  
    coh25:k=3.06122,b=3.06122        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        6.1224        6.1224        6.1224        7.1429        7.1429       7.1429  
    coh26:k=2.04082,b=4.08163        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        6.1224        6.1224        7.1429        7.1429        7.1429       8.1633  
    coh27:k=1.02041,b=5.10204        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        6.1224        6.1224        7.1429        7.1429        7.1429       8.1633  
    coh28:k=0,b=6.12245              2.0408        2.0408       3.0612       3.0612        4.0816         5.102        6.1224        6.1224        7.1429        7.1429        8.1633       8.1633  
    coh29:k=7.14286,b=0              2.0408        2.0408       3.0612       3.0612        4.0816         5.102        6.1224         5.102         5.102         5.102         5.102        5.102  
    coh30:k=6.12245,b=1.02041        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        6.1224        7.1429        7.1429        7.1429        7.1429       7.1429  
    coh31:k=5.10204,b=2.04082        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        7.1429        7.1429        7.1429        7.1429        8.1633       8.1633  
    coh32:k=4.08163,b=3.06122        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        7.1429        7.1429        7.1429        8.1633        8.1633       8.1633  
    coh33:k=3.06122,b=4.08163        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        7.1429        7.1429        7.1429        8.1633        8.1633       9.1837  
    coh34:k=2.04082,b=5.10204        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        7.1429        7.1429        8.1633        8.1633        8.1633       9.1837  
    coh35:k=1.02041,b=6.12245        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        7.1429        7.1429        8.1633        8.1633        8.1633       9.1837  
    coh36:k=0,b=7.14286              2.0408        2.0408       3.0612       3.0612        4.0816         5.102        7.1429        7.1429        8.1633        8.1633        9.1837       9.1837  
    coh37:k=8.16327,b=0              2.0408        2.0408       3.0612       3.0612        4.0816         5.102        6.1224        6.1224        6.1224        6.1224        6.1224       6.1224  
    coh38:k=7.14286,b=1.02041        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        7.1429        8.1633        8.1633        8.1633        8.1633       8.1633  
    coh39:k=6.12245,b=2.04082        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        8.1633        8.1633        8.1633        8.1633        9.1837       9.1837  
    coh40:k=5.10204,b=3.06122        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        8.1633        8.1633        8.1633        8.1633        9.1837       9.1837  
    coh41:k=4.08163,b=4.08163        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        8.1633        8.1633        8.1633        9.1837        9.1837       9.1837  
    coh42:k=3.06122,b=5.10204        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        8.1633        8.1633        8.1633        9.1837        9.1837       10.204  
    coh43:k=2.04082,b=6.12245        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        8.1633        8.1633        9.1837        9.1837        9.1837       10.204  
    coh44:k=1.02041,b=7.14286        2.0408        2.0408       3.0612       3.0612        4.0816         5.102        8.1633        8.1633        9.1837        9.1837        9.1837       10.204  
    coh45:k=0,b=8.16327              2.0408        2.0408       3.0612       3.0612        4.0816         5.102        8.1633        8.1633        9.1837        9.1837        10.204       10.204  
    coh46:k=9.18367,b=0              2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        7.1429        7.1429        7.1429        7.1429        7.1429       7.1429  
    coh47:k=8.16327,b=1.02041        2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        8.1633        8.1633        9.1837        9.1837        9.1837       9.1837  
    coh48:k=7.14286,b=2.04082        2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        9.1837        9.1837        9.1837        9.1837       10.204  
    coh49:k=6.12245,b=3.06122        2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        9.1837        9.1837        9.1837        10.204       10.204  
    coh50:k=5.10204,b=4.08163        2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        9.1837        9.1837        10.204        10.204       10.204  
    coh1226:k=50,b=0                 2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       21.429  
    coh1227:k=48.9796,b=1.02041      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       21.429  
    coh1228:k=47.9592,b=2.04082      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       21.429  
    coh1229:k=46.9388,b=3.06122      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       21.429  
    coh1230:k=45.9184,b=4.08163      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1231:k=44.898,b=5.10204       2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1232:k=43.8776,b=6.12245      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1233:k=42.8571,b=7.14286      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1234:k=41.8367,b=8.16327      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1235:k=40.8163,b=9.18367      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1236:k=39.7959,b=10.2041      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1237:k=38.7755,b=11.2245      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1238:k=37.7551,b=12.2449      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1239:k=36.7347,b=13.2653      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1240:k=35.7143,b=14.2857      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1241:k=34.6939,b=15.3061      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1242:k=33.6735,b=16.3265      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1243:k=32.6531,b=17.3469      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1244:k=31.6327,b=18.3673      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1245:k=30.6122,b=19.3878      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1246:k=29.5918,b=20.4082      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1247:k=28.5714,b=21.4286      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1248:k=27.551,b=22.449        2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1249:k=26.5306,b=23.4694      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1250:k=25.5102,b=24.4898      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1251:k=24.4898,b=25.5102      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1252:k=23.4694,b=26.5306      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1253:k=22.449,b=27.551        2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1254:k=21.4286,b=28.5714      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1255:k=20.4082,b=29.5918      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1256:k=19.3878,b=30.6122      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1257:k=18.3673,b=31.6327      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1258:k=17.3469,b=32.6531      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1259:k=16.3265,b=33.6735      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1260:k=15.3061,b=34.6939      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1261:k=14.2857,b=35.7143      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1262:k=13.2653,b=36.7347      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1263:k=12.2449,b=37.7551      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1264:k=11.2245,b=38.7755      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1265:k=10.2041,b=39.7959      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1266:k=9.18367,b=40.8163      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1267:k=8.16327,b=41.8367      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1268:k=7.14286,b=42.8571      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1269:k=6.12245,b=43.8776      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1270:k=5.10204,b=44.898       2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1271:k=4.08163,b=45.9184      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1272:k=3.06122,b=46.9388      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1273:k=2.04082,b=47.9592      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1274:k=1.02041,b=48.9796      2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  
    coh1275:k=0,b=50                 2.0408        2.0408       3.0612       3.0612        4.0816        4.0816        9.1837        11.224        13.265        15.306        18.367       20.408  

tb_pol_w: risky + safe investment choices (first stage choice, choose within risky vs safe)
                                   z1_0_34741    z2_0_40076    z3_0_4623    z4_0_5333    z5_0_61519    z6_0_70966    z10_1_2567    z11_1_4496    z12_1_6723    z13_1_9291    z14_2_2253    z15_2_567
                                   __________    __________    _________    _________    __________    __________    __________    __________    __________    __________    __________    _________

    coh1:k=0,b=0                          0             0            0            0             0             0             0             0             0             0             0            0  
    coh2:k=1.02041,b=0               1.0204        1.0204       1.0204       1.0204        1.0204        1.0204        1.0204        1.0204        1.0204        1.0204        1.0204       1.0204  
    coh3:k=0,b=1.02041               1.0204        1.0204       1.0204       1.0204        1.0204        1.0204        2.0408        2.0408        2.0408        2.0408        2.0408       2.0408  
    coh4:k=2.04082,b=0               2.0408        2.0408       2.0408       2.0408        2.0408        2.0408        1.0204        1.0204        1.0204        1.0204        1.0204       1.0204  
    coh5:k=1.02041,b=1.02041         2.0408        2.0408       2.0408       2.0408        2.0408        2.0408        2.0408        2.0408        3.0612        3.0612        3.0612       3.0612  
    coh6:k=0,b=2.04082               2.0408        2.0408       2.0408       2.0408        2.0408        2.0408        2.0408        3.0612        3.0612        3.0612        3.0612       4.0816  
    coh7:k=3.06122,b=0               2.0408        2.0408       2.0408       2.0408        2.0408        2.0408        2.0408        2.0408        2.0408        2.0408        2.0408       2.0408  
    coh8:k=2.04082,b=1.02041         2.0408        2.0408       2.0408       3.0612        3.0612        3.0612        3.0612        3.0612        3.0612        3.0612        4.0816       4.0816  
    coh9:k=1.02041,b=2.04082         2.0408        2.0408       3.0612       3.0612        3.0612        3.0612        3.0612        3.0612        4.0816        4.0816        4.0816       4.0816  
    coh10:k=0,b=3.06122              2.0408        2.0408       3.0612       3.0612        3.0612        3.0612        3.0612        4.0816        4.0816        4.0816        4.0816        5.102  
    coh11:k=4.08163,b=0              3.0612        3.0612       3.0612       3.0612        3.0612        3.0612        3.0612        3.0612        3.0612        3.0612        3.0612       3.0612  
    coh12:k=3.06122,b=1.02041        3.0612        3.0612       3.0612       3.0612        3.0612        3.0612        4.0816        4.0816        4.0816        4.0816         5.102        5.102  
    coh13:k=2.04082,b=2.04082        3.0612        3.0612       3.0612       3.0612        3.0612        4.0816        4.0816        4.0816        4.0816         5.102         5.102        5.102  
    coh14:k=1.02041,b=3.06122        3.0612        3.0612       3.0612       3.0612        4.0816        4.0816        4.0816        4.0816         5.102         5.102         5.102       6.1224  
    coh15:k=0,b=4.08163              3.0612        3.0612       3.0612       3.0612        3.0612        4.0816        4.0816         5.102         5.102         5.102         5.102       6.1224  
    coh16:k=5.10204,b=0              4.0816        4.0816       4.0816       4.0816        4.0816        4.0816        4.0816        4.0816        4.0816        4.0816        4.0816       4.0816  
    coh17:k=4.08163,b=1.02041        4.0816        4.0816       4.0816       4.0816        4.0816        4.0816         5.102         5.102         5.102         5.102         5.102       6.1224  
    coh18:k=3.06122,b=2.04082        4.0816        4.0816       4.0816       4.0816        4.0816        4.0816         5.102         5.102         5.102        6.1224        6.1224       6.1224  
    coh19:k=2.04082,b=3.06122        4.0816        4.0816       4.0816       4.0816        4.0816        4.0816         5.102         5.102        6.1224        6.1224        6.1224       7.1429  
    coh20:k=1.02041,b=4.08163        4.0816        4.0816       4.0816       4.0816        4.0816        4.0816         5.102         5.102        6.1224        6.1224        6.1224       7.1429  
    coh21:k=0,b=5.10204              4.0816        4.0816       4.0816       4.0816        4.0816        4.0816         5.102         5.102        6.1224        6.1224        7.1429       7.1429  
    coh22:k=6.12245,b=0               5.102         5.102        5.102        5.102         5.102         5.102         5.102         5.102         5.102         5.102         5.102       4.0816  
    coh23:k=5.10204,b=1.02041         5.102         5.102        5.102        5.102         5.102         5.102        6.1224        6.1224        6.1224        6.1224        6.1224       7.1429  
    coh24:k=4.08163,b=2.04082         5.102         5.102        5.102        5.102         5.102         5.102        6.1224        6.1224        6.1224        6.1224        7.1429       7.1429  
    coh25:k=3.06122,b=3.06122         5.102         5.102        5.102        5.102         5.102         5.102        6.1224        6.1224        6.1224        7.1429        7.1429       7.1429  
    coh26:k=2.04082,b=4.08163         5.102         5.102        5.102        5.102         5.102         5.102        6.1224        6.1224        7.1429        7.1429        7.1429       8.1633  
    coh27:k=1.02041,b=5.10204         5.102         5.102        5.102        5.102         5.102         5.102        6.1224        6.1224        7.1429        7.1429        7.1429       8.1633  
    coh28:k=0,b=6.12245               5.102         5.102        5.102        5.102         5.102         5.102        6.1224        6.1224        7.1429        7.1429        8.1633       8.1633  
    coh29:k=7.14286,b=0              6.1224        6.1224       6.1224       6.1224        6.1224        6.1224        6.1224         5.102         5.102         5.102         5.102        5.102  
    coh30:k=6.12245,b=1.02041        6.1224        6.1224       6.1224       6.1224        6.1224        6.1224        6.1224        7.1429        7.1429        7.1429        7.1429       7.1429  
    coh31:k=5.10204,b=2.04082        6.1224        6.1224       6.1224       6.1224        6.1224        6.1224        7.1429        7.1429        7.1429        7.1429        8.1633       8.1633  
    coh32:k=4.08163,b=3.06122        6.1224        6.1224       6.1224       6.1224        6.1224        6.1224        7.1429        7.1429        7.1429        8.1633        8.1633       8.1633  
    coh33:k=3.06122,b=4.08163        6.1224        6.1224       6.1224       6.1224        6.1224        6.1224        7.1429        7.1429        7.1429        8.1633        8.1633       9.1837  
    coh34:k=2.04082,b=5.10204        6.1224        6.1224       6.1224       6.1224        6.1224        6.1224        7.1429        7.1429        8.1633        8.1633        8.1633       9.1837  
    coh35:k=1.02041,b=6.12245        6.1224        6.1224       6.1224       6.1224        6.1224        6.1224        7.1429        7.1429        8.1633        8.1633        8.1633       9.1837  
    coh36:k=0,b=7.14286              6.1224        6.1224       6.1224       6.1224        6.1224        6.1224        7.1429        7.1429        8.1633        8.1633        9.1837       9.1837  
    coh37:k=8.16327,b=0              7.1429        7.1429       7.1429       7.1429        7.1429        7.1429        6.1224        6.1224        6.1224        6.1224        6.1224       6.1224  
    coh38:k=7.14286,b=1.02041        7.1429        7.1429       7.1429       7.1429        7.1429        7.1429        7.1429        8.1633        8.1633        8.1633        8.1633       8.1633  
    coh39:k=6.12245,b=2.04082        7.1429        7.1429       7.1429       7.1429        7.1429        7.1429        8.1633        8.1633        8.1633        8.1633        9.1837       9.1837  
    coh40:k=5.10204,b=3.06122        7.1429        7.1429       7.1429       7.1429        7.1429        7.1429        8.1633        8.1633        8.1633        8.1633        9.1837       9.1837  
    coh41:k=4.08163,b=4.08163        7.1429        7.1429       7.1429       7.1429        7.1429        7.1429        8.1633        8.1633        8.1633        9.1837        9.1837       9.1837  
    coh42:k=3.06122,b=5.10204        7.1429        7.1429       7.1429       7.1429        7.1429        7.1429        8.1633        8.1633        8.1633        9.1837        9.1837       10.204  
    coh43:k=2.04082,b=6.12245        7.1429        7.1429       7.1429       7.1429        7.1429        7.1429        8.1633        8.1633        9.1837        9.1837        9.1837       10.204  
    coh44:k=1.02041,b=7.14286        7.1429        7.1429       7.1429       7.1429        7.1429        7.1429        8.1633        8.1633        9.1837        9.1837        9.1837       10.204  
    coh45:k=0,b=8.16327              7.1429        7.1429       7.1429       7.1429        7.1429        7.1429        8.1633        8.1633        9.1837        9.1837        10.204       10.204  
    coh46:k=9.18367,b=0              8.1633        8.1633       8.1633       8.1633        8.1633        8.1633        7.1429        7.1429        7.1429        7.1429        7.1429       7.1429  
    coh47:k=8.16327,b=1.02041        8.1633        8.1633       8.1633       8.1633        8.1633        8.1633        8.1633        8.1633        9.1837        9.1837        9.1837       9.1837  
    coh48:k=7.14286,b=2.04082        8.1633        8.1633       8.1633       8.1633        8.1633        8.1633        9.1837        9.1837        9.1837        9.1837        9.1837       10.204  
    coh49:k=6.12245,b=3.06122        8.1633        8.1633       8.1633       8.1633        8.1633        8.1633        9.1837        9.1837        9.1837        9.1837        10.204       10.204  
    coh50:k=5.10204,b=4.08163        8.1633        8.1633       8.1633       8.1633        8.1633        8.1633        9.1837        9.1837        9.1837        10.204        10.204       10.204  
    coh1226:k=50,b=0                 46.939        46.939       46.939       46.939        46.939        46.939        46.939        46.939        45.918        45.918        45.918       45.918  
    coh1227:k=48.9796,b=1.02041      46.939        46.939       46.939       46.939        46.939        46.939        47.959        47.959        47.959        47.959        47.959       47.959  
    coh1228:k=47.9592,b=2.04082      46.939        46.939       46.939       46.939        46.939        46.939        47.959        47.959        47.959        47.959         48.98        48.98  
    coh1229:k=46.9388,b=3.06122      46.939        46.939       46.939       46.939        46.939        46.939        47.959        47.959        47.959         48.98         48.98        48.98  
    coh1230:k=45.9184,b=4.08163      46.939        46.939       46.939       46.939        46.939        46.939        47.959        47.959         48.98         48.98         48.98           50  
    coh1231:k=44.898,b=5.10204       46.939        46.939       46.939       46.939        46.939        46.939        47.959        47.959         48.98         48.98         48.98           50  
    coh1232:k=43.8776,b=6.12245      46.939        46.939       46.939       46.939        46.939        46.939        47.959         48.98         48.98         48.98            50           50  
    coh1233:k=42.8571,b=7.14286      46.939        46.939       46.939       46.939        46.939        46.939        47.959         48.98         48.98         48.98            50           50  
    coh1234:k=41.8367,b=8.16327      46.939        46.939       46.939       46.939        46.939        46.939        47.959         48.98         48.98         48.98            50           50  
    coh1235:k=40.8163,b=9.18367      46.939        46.939       46.939       46.939        46.939        46.939        47.959         48.98         48.98         48.98            50           50  
    coh1236:k=39.7959,b=10.2041      46.939        46.939       46.939       46.939        46.939        46.939        47.959         48.98         48.98         48.98            50           50  
    coh1237:k=38.7755,b=11.2245      46.939        46.939       46.939       46.939        46.939        46.939        47.959         48.98         48.98            50            50           50  
    coh1238:k=37.7551,b=12.2449      46.939        46.939       46.939       46.939        46.939        46.939        47.959         48.98         48.98            50            50           50  
    coh1239:k=36.7347,b=13.2653      46.939        46.939       46.939       46.939        46.939        46.939        47.959         48.98         48.98            50            50           50  
    coh1240:k=35.7143,b=14.2857      45.918        46.939       46.939       46.939        46.939        46.939        47.959         48.98         48.98            50            50           50  
    coh1241:k=34.6939,b=15.3061      45.918        45.918       46.939       46.939        46.939        46.939        47.959         48.98         48.98            50            50           50  
    coh1242:k=33.6735,b=16.3265      45.918        45.918       46.939       46.939        46.939        46.939        47.959         48.98         48.98            50            50           50  
    coh1243:k=32.6531,b=17.3469      45.918        45.918       45.918       46.939        46.939        46.939        47.959         48.98         48.98            50            50           50  
    coh1244:k=31.6327,b=18.3673      45.918        45.918       45.918       46.939        46.939        46.939        47.959         48.98         48.98            50            50           50  
    coh1245:k=30.6122,b=19.3878      45.918        45.918       45.918       45.918        46.939        46.939        47.959         48.98         48.98            50            50           50  
    coh1246:k=29.5918,b=20.4082      45.918        45.918       45.918       45.918        46.939        46.939        47.959         48.98         48.98            50            50           50  
    coh1247:k=28.5714,b=21.4286      45.918        45.918       45.918       45.918        45.918        46.939        47.959         48.98         48.98            50            50           50  
    coh1248:k=27.551,b=22.449        45.918        45.918       45.918       45.918        45.918        46.939        47.959        47.959         48.98            50            50           50  
    coh1249:k=26.5306,b=23.4694      45.918        45.918       45.918       45.918        45.918        46.939        47.959        47.959         48.98            50            50           50  
    coh1250:k=25.5102,b=24.4898      45.918        45.918       45.918       45.918        45.918        46.939        47.959        47.959         48.98            50            50           50  
    coh1251:k=24.4898,b=25.5102      45.918        45.918       45.918       45.918        45.918        45.918        47.959        47.959         48.98            50            50           50  
    coh1252:k=23.4694,b=26.5306      44.898        45.918       45.918       45.918        45.918        45.918        47.959        47.959         48.98            50            50           50  
    coh1253:k=22.449,b=27.551        44.898        45.918       45.918       45.918        45.918        45.918        47.959        47.959         48.98            50            50           50  
    coh1254:k=21.4286,b=28.5714      44.898        44.898       45.918       45.918        45.918        45.918        47.959        47.959         48.98            50            50           50  
    coh1255:k=20.4082,b=29.5918      44.898        44.898       44.898       45.918        45.918        45.918        47.959        47.959         48.98            50            50           50  
    coh1256:k=19.3878,b=30.6122      44.898        44.898       44.898       45.918        45.918        45.918        47.959        47.959         48.98         48.98            50           50  
    coh1257:k=18.3673,b=31.6327      44.898        44.898       44.898       45.918        45.918        45.918        47.959        47.959         48.98         48.98            50           50  
    coh1258:k=17.3469,b=32.6531      44.898        44.898       44.898       44.898        45.918        45.918        47.959        47.959         48.98         48.98            50           50  
    coh1259:k=16.3265,b=33.6735      44.898        44.898       44.898       44.898        45.918        45.918        46.939        47.959         48.98         48.98            50           50  
    coh1260:k=15.3061,b=34.6939      44.898        44.898       44.898       44.898        45.918        45.918        46.939        47.959         48.98         48.98            50           50  
    coh1261:k=14.2857,b=35.7143      44.898        44.898       44.898       44.898        44.898        45.918        46.939        47.959         48.98         48.98            50           50  
    coh1262:k=13.2653,b=36.7347      44.898        44.898       44.898       44.898        44.898        45.918        46.939        47.959         48.98         48.98            50           50  
    coh1263:k=12.2449,b=37.7551      44.898        44.898       44.898       44.898        44.898        45.918        46.939        47.959         48.98         48.98            50           50  
    coh1264:k=11.2245,b=38.7755      43.878        44.898       44.898       44.898        44.898        44.898        46.939        47.959        47.959         48.98            50           50  
    coh1265:k=10.2041,b=39.7959      43.878        43.878       44.898       44.898        44.898        44.898        46.939        47.959        47.959         48.98            50           50  
    coh1266:k=9.18367,b=40.8163      43.878        43.878       44.898       44.898        44.898        44.898        46.939        47.959        47.959         48.98            50           50  
    coh1267:k=8.16327,b=41.8367      43.878        43.878       43.878       44.898        44.898        44.898        46.939        47.959        47.959         48.98            50           50  
    coh1268:k=7.14286,b=42.8571      43.878        43.878       43.878       44.898        44.898        44.898        46.939        47.959        47.959         48.98            50           50  
    coh1269:k=6.12245,b=43.8776      43.878        43.878       43.878       44.898        44.898        44.898        46.939        47.959        47.959         48.98            50           50  
    coh1270:k=5.10204,b=44.898       43.878        43.878       43.878       43.878        44.898        44.898        46.939        46.939        47.959         48.98            50           50  
    coh1271:k=4.08163,b=45.9184      43.878        43.878       43.878       43.878        44.898        44.898        46.939        46.939        47.959         48.98            50           50  
    coh1272:k=3.06122,b=46.9388      43.878        43.878       43.878       43.878        44.898        44.898        46.939        46.939        47.959         48.98            50           50  
    coh1273:k=2.04082,b=47.9592      43.878        43.878       43.878       43.878        43.878        44.898        46.939        46.939        47.959         48.98            50           50  
    coh1274:k=1.02041,b=48.9796      43.878        43.878       43.878       43.878        43.878        44.898        46.939        46.939        47.959         48.98            50           50  
    coh1275:k=0,b=50                 42.857        43.878       43.878       43.878        43.878        44.898        46.939        46.939        47.959         48.98            50           50  

end
ans = 

  Map with properties:

        Count: 14
      KeyType: char
    ValueType: any