Contents

ipwkbz Stock vs Risky Capital Investment

close all

Simulate Risky Asset with CRS and full depreciation, IID shocks

CRRA utility with Constant return to scale risky asset and safe asset at the parameters below, constant share of asset in riksy asset, at the parameters below, shock does not have additional effects except through coh so plotting conditional on cash-on-hand, there is no difference across shocks, hence solve at only it_z_n == 3. Note we have to have many grid points for it_w_n to get the exact solution. For pratical purposes, this does not need to be as high.

it_param_set = 4;
[param_map, support_map] = ffs_ipwkbz_set_default_param(it_param_set);

% Simulation Accuracy
param_map('it_w_perc_n') = 500;
param_map('it_ak_perc_n') = param_map('it_w_perc_n');
param_map('it_z_n') = 15;

param_map('fl_coh_interp_grid_gap') = 0.01;
param_map('fl_w_interp_grid_gap') = 0.01;
param_map('it_c_interp_grid_gap') = 10^-4;

% Iteration Limits
% param_map('it_maxiter_val') = 2;

% Turn on 2nd stage graphs
support_map('bl_graph_evf') = true;
support_map('bl_display_evf') = false;

% Production Function Parameters
% note shock is log normal
param_map('fl_Amean') = 1.0265;
param_map('fl_alpha') = 1;
param_map('fl_delta') = 1;
param_map('fl_r') = 0.03;
param_map('fl_w') = 0.20;

% Shock Parameter, iid shocks
param_map('fl_z_rho') = 0;
param_map('fl_z_sig') = 0.05;

% Display Parameters
support_map('bl_display') = false;
support_map('bl_display_final') = true;
support_map('bl_time') = true;
% support_map('bl_profile') = false;


% Call Program
ff_ipwkbz_vf_vecsv(param_map, support_map);
Warning:
"C:\Users\fan\CodeDynaAsset\m_ipwkbz\test\ff_ipwkbz_vf_vecsv\test_prod\html"
not found in path. 
Warning: Ignoring extra legend entries. 
Elapsed time is 651.191191 seconds.
Warning: Directory already exists. 
valgap = norm(mt_val - mt_val_cur)
polgap = norm(mt_pol_a - mt_pol_a_cur)
z1 = z1 perc change: sum((mt_pol_a ~= mt_pol_a_cur))/(it_ameshk_n)
                  valgap      polgap         z1            z2            z3            z4            z5            z6            z7            z8            z9           z10           z11           z12           z13           z14           z15    
                __________    _______    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________

    iter=1          628.62     7044.4             1             1             1             1             1             1             1             1             1             1             1             1             1             1             1
    iter=2          1296.2     8289.9       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904
    iter=3          1118.6     2772.2       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904       0.99904
    iter=4          973.78     1388.9       0.99506       0.99506       0.99506       0.99506       0.99506       0.99506       0.99506       0.99506       0.99506       0.99506       0.99506       0.99506       0.99506       0.99506       0.99506
    iter=5          854.92     828.95       0.99121       0.99121       0.99121       0.99121       0.99121       0.99121       0.99121       0.99121       0.99121       0.99121       0.99121       0.99121       0.99121       0.99121       0.99121
    iter=6          752.58     557.73       0.98796       0.98796       0.98796       0.98796       0.98796       0.98796       0.98796       0.98796       0.98796       0.98796       0.98796       0.98796       0.98796       0.98796       0.98796
    iter=7          672.49     410.25       0.98387       0.98387       0.98387       0.98387       0.98387       0.98387       0.98387       0.98387       0.98387       0.98387       0.98387       0.98387       0.98387       0.98387       0.98387
    iter=8          603.02     298.84       0.98038       0.98038       0.98038       0.98038       0.98038       0.98038       0.98038       0.98038       0.98038       0.98038       0.98038       0.98038       0.98038       0.98038       0.98038
    iter=9          546.59     252.06       0.97628       0.97628       0.97628       0.97628       0.97628       0.97628       0.97628       0.97628       0.97628       0.97628       0.97628       0.97628       0.97628       0.97628       0.97628
    iter=10         497.44     215.06       0.94486       0.94486       0.94486       0.94486       0.94486       0.94486       0.94486       0.94486       0.94486       0.94486       0.94486       0.94486       0.94486       0.94486       0.94486
    iter=11         455.82     226.26       0.96942       0.96942       0.96942       0.96942       0.96942       0.96942       0.96942       0.96942       0.96942       0.96942       0.96942       0.96942       0.96942       0.96942       0.96942
    iter=12         419.17     207.75       0.96737       0.96737       0.96737       0.96737       0.96737       0.96737       0.96737       0.96737       0.96737       0.96737       0.96737       0.96737       0.96737       0.96737       0.96737
    iter=13          387.2     229.72       0.90573       0.90573       0.90573       0.90573       0.90573       0.90573       0.90573       0.90573       0.90573       0.90573       0.90573       0.90573       0.90573       0.90573       0.90573
    iter=14         358.71     269.16       0.96316       0.96316       0.96316       0.96316       0.96316       0.96316       0.96316       0.96316       0.96316       0.96316       0.96316       0.96316       0.96316       0.96316       0.96316
    iter=15         333.37     256.64       0.89201       0.89201       0.89201       0.89201       0.89201       0.89201       0.89201       0.89201       0.89201       0.89201       0.89201       0.89201       0.89201       0.89201       0.89201
    iter=16         310.61     277.22       0.96111       0.96111       0.96111       0.96111       0.96111       0.96111       0.96111       0.96111       0.96111       0.96111       0.96111       0.96111       0.96111       0.96111       0.96111
    iter=17         290.11     295.71       0.96232       0.96232       0.96232       0.96232       0.96232       0.96232       0.96232       0.96232       0.96232       0.96232       0.96232       0.96232       0.96232       0.96232       0.96232
    iter=18         271.56     314.78        0.8784        0.8784        0.8784        0.8784        0.8784        0.8784        0.8784        0.8784        0.8784        0.8784        0.8784        0.8784        0.8784        0.8784        0.8784
    iter=19         254.71     320.61       0.87455       0.87455       0.87455       0.87455       0.87455       0.87455       0.87455       0.87455       0.87455       0.87455       0.87455       0.87455       0.87455       0.87455       0.87455
    iter=20         239.34     340.85       0.95762       0.95762       0.95762       0.95762       0.95762       0.95762       0.95762       0.95762       0.95762       0.95762       0.95762       0.95762       0.95762       0.95762       0.95762
    iter=21         225.28     361.84       0.89791       0.89791       0.89791       0.89791       0.89791       0.89791       0.89791       0.89791       0.89791       0.89791       0.89791       0.89791       0.89791       0.89791       0.89791
    iter=22         212.39     400.29        0.9256        0.9256        0.9256        0.9256        0.9256        0.9256        0.9256        0.9256        0.9256        0.9256        0.9256        0.9256        0.9256        0.9256        0.9256
    iter=23         200.53     394.06       0.83879       0.83879       0.83879       0.83879       0.83879       0.83879       0.83879       0.83879       0.83879       0.83879       0.83879       0.83879       0.83879       0.83879       0.83879
    iter=24         189.59     416.38       0.82747       0.82747       0.82747       0.82747       0.82747       0.82747       0.82747       0.82747       0.82747       0.82747       0.82747       0.82747       0.82747       0.82747       0.82747
    iter=25         179.48     481.23        0.8837        0.8837        0.8837        0.8837        0.8837        0.8837        0.8837        0.8837        0.8837        0.8837        0.8837        0.8837        0.8837        0.8837        0.8837
    iter=26         170.11     446.93        0.8117        0.8117        0.8117        0.8117        0.8117        0.8117        0.8117        0.8117        0.8117        0.8117        0.8117        0.8117        0.8117        0.8117        0.8117
    iter=27         161.42     476.71         0.853         0.853         0.853         0.853         0.853         0.853         0.853         0.853         0.853         0.853         0.853         0.853         0.853         0.853         0.853
    iter=28         153.33     490.36       0.83482       0.83482       0.83482       0.83482       0.83482       0.83482       0.83482       0.83482       0.83482       0.83482       0.83482       0.83482       0.83482       0.83482       0.83482
    iter=29          145.8     492.46       0.77257       0.77257       0.77257       0.77257       0.77257       0.77257       0.77257       0.77257       0.77257       0.77257       0.77257       0.77257       0.77257       0.77257       0.77257
    iter=30         138.77     536.51       0.81447       0.81447       0.81447       0.81447       0.81447       0.81447       0.81447       0.81447       0.81447       0.81447       0.81447       0.81447       0.81447       0.81447       0.81447
    iter=31         132.21     513.19       0.77968       0.77968       0.77968       0.77968       0.77968       0.77968       0.77968       0.77968       0.77968       0.77968       0.77968       0.77968       0.77968       0.77968       0.77968
    iter=32         126.06     524.07       0.69709       0.69709       0.69709       0.69709       0.69709       0.69709       0.69709       0.69709       0.69709       0.69709       0.69709       0.69709       0.69709       0.69709       0.69709
    iter=33          120.3     540.88       0.71382       0.71382       0.71382       0.71382       0.71382       0.71382       0.71382       0.71382       0.71382       0.71382       0.71382       0.71382       0.71382       0.71382       0.71382
    iter=34          114.9     570.03       0.78184       0.78184       0.78184       0.78184       0.78184       0.78184       0.78184       0.78184       0.78184       0.78184       0.78184       0.78184       0.78184       0.78184       0.78184
    iter=35         109.82     544.95       0.67229       0.67229       0.67229       0.67229       0.67229       0.67229       0.67229       0.67229       0.67229       0.67229       0.67229       0.67229       0.67229       0.67229       0.67229
    iter=36         105.04     556.28       0.68721       0.68721       0.68721       0.68721       0.68721       0.68721       0.68721       0.68721       0.68721       0.68721       0.68721       0.68721       0.68721       0.68721       0.68721
    iter=37         100.54     557.83       0.71755       0.71755       0.71755       0.71755       0.71755       0.71755       0.71755       0.71755       0.71755       0.71755       0.71755       0.71755       0.71755       0.71755       0.71755
    iter=38         96.294     581.88       0.75295       0.75295       0.75295       0.75295       0.75295       0.75295       0.75295       0.75295       0.75295       0.75295       0.75295       0.75295       0.75295       0.75295       0.75295
    iter=39         92.285     547.73       0.68481       0.68481       0.68481       0.68481       0.68481       0.68481       0.68481       0.68481       0.68481       0.68481       0.68481       0.68481       0.68481       0.68481       0.68481
    iter=40         88.496     545.24       0.67493       0.67493       0.67493       0.67493       0.67493       0.67493       0.67493       0.67493       0.67493       0.67493       0.67493       0.67493       0.67493       0.67493       0.67493
    iter=41         84.913     540.53       0.68023       0.68023       0.68023       0.68023       0.68023       0.68023       0.68023       0.68023       0.68023       0.68023       0.68023       0.68023       0.68023       0.68023       0.68023
    iter=42         81.515     531.71       0.65122       0.65122       0.65122       0.65122       0.65122       0.65122       0.65122       0.65122       0.65122       0.65122       0.65122       0.65122       0.65122       0.65122       0.65122
    iter=43          78.29     573.65       0.78534       0.78534       0.78534       0.78534       0.78534       0.78534       0.78534       0.78534       0.78534       0.78534       0.78534       0.78534       0.78534       0.78534       0.78534
    iter=44         75.226     513.03       0.65627       0.65627       0.65627       0.65627       0.65627       0.65627       0.65627       0.65627       0.65627       0.65627       0.65627       0.65627       0.65627       0.65627       0.65627
    iter=45         72.251     511.01       0.67638       0.67638       0.67638       0.67638       0.67638       0.67638       0.67638       0.67638       0.67638       0.67638       0.67638       0.67638       0.67638       0.67638       0.67638
    iter=46         69.249      503.4       0.65603       0.65603       0.65603       0.65603       0.65603       0.65603       0.65603       0.65603       0.65603       0.65603       0.65603       0.65603       0.65603       0.65603       0.65603
    iter=47         66.205     487.01       0.65796       0.65796       0.65796       0.65796       0.65796       0.65796       0.65796       0.65796       0.65796       0.65796       0.65796       0.65796       0.65796       0.65796       0.65796
    iter=48         63.153     473.46       0.64772       0.64772       0.64772       0.64772       0.64772       0.64772       0.64772       0.64772       0.64772       0.64772       0.64772       0.64772       0.64772       0.64772       0.64772
    iter=49          60.13     463.59       0.66446       0.66446       0.66446       0.66446       0.66446       0.66446       0.66446       0.66446       0.66446       0.66446       0.66446       0.66446       0.66446       0.66446       0.66446
    iter=50         57.164     451.94       0.66025       0.66025       0.66025       0.66025       0.66025       0.66025       0.66025       0.66025       0.66025       0.66025       0.66025       0.66025       0.66025       0.66025       0.66025
    iter=201     0.0062558    0.86677    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158
    iter=202      0.005869    0.71489    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237
    iter=203      0.005502    0.45069    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079
    iter=204     0.0051545    0.73791    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237    0.00072237
    iter=205     0.0048262    0.17413    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079
    iter=206     0.0045166    0.78512    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158
    iter=207     0.0042255     0.4529    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079
    iter=208     0.0039525    0.87986    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158
    iter=209     0.0036971     1.0646    0.00096316    0.00096316    0.00096316    0.00096316    0.00096316    0.00096316    0.00096316    0.00096316    0.00096316    0.00096316    0.00096316    0.00096316    0.00096316    0.00096316    0.00096316
    iter=210     0.0034586          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=211     0.0032365    0.46888    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079
    iter=212     0.0030301    0.50298    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079
    iter=213     0.0028384          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=214     0.0026607    0.36356    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079
    iter=215      0.002496          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=216     0.0023435          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=217      0.002202    0.46986    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079
    iter=218     0.0020708          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=219     0.0019489    0.51398    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079
    iter=220     0.0018355      0.437    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079
    iter=221     0.0017298          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=222      0.001631    0.68505    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158
    iter=223     0.0015385    0.20826    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079
    iter=224     0.0014517          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=225       0.00137          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=226      0.001293          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=227     0.0012203          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=228     0.0011515     0.2407    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079    0.00024079
    iter=229     0.0010863          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=230     0.0010245          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=231    0.00096571          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=232    0.00090988          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=233    0.00085681    0.78495    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158
    iter=234    0.00080635          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=235    0.00075838          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=236     0.0007128          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=237    0.00066952    0.79418    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158    0.00048158
    iter=238    0.00062848          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=239     0.0005896          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=240    0.00055283          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=241    0.00051811          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=242    0.00048539          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=243    0.00045461          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=244    0.00042571          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=245    0.00039863          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=246     0.0003733          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=247    0.00034965          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=248    0.00032761          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=249    0.00030708          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=250    0.00028798          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0

tb_val
                             z1_0_85962    z2_0_87824    z3_0_89726    z4_0_9167    z5_0_93655    z6_0_95684    z10_1_0425    z11_1_0651    z12_1_0881    z13_1_1117    z14_1_1358    z15_1_1604
                             __________    __________    __________    _________    __________    __________    __________    __________    __________    __________    __________    __________

    coh1:k=-20,b=0            -77.899       -77.899       -77.899       -77.899      -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899  
    coh2:k=-19.99,b=0         -77.899       -77.899       -77.899       -77.899      -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899  
    coh3:k=-19.98,b=0         -77.899       -77.899       -77.899       -77.899      -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899  
    coh4:k=-19.97,b=0         -77.899       -77.899       -77.899       -77.899      -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899  
    coh5:k=-19.96,b=0         -77.899       -77.899       -77.899       -77.899      -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899  
    coh6:k=-19.95,b=0         -77.899       -77.899       -77.899       -77.899      -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899  
    coh7:k=-19.94,b=0         -77.899       -77.899       -77.899       -77.899      -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899  
    coh8:k=-19.93,b=0         -77.899       -77.899       -77.899       -77.899      -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899  
    coh9:k=-19.92,b=0         -77.899       -77.899       -77.899       -77.899      -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899       -77.899  
    coh10:k=-19.91,b=0        -77.894       -77.894       -77.894       -77.894      -77.894       -77.894       -77.894       -77.894       -77.894       -77.894       -77.894       -77.894  
    coh11:k=-19.9,b=0         -77.552       -77.552       -77.552       -77.552      -77.552       -77.552       -77.552       -77.552       -77.552       -77.552       -77.552       -77.552  
    coh12:k=-19.89,b=0        -77.257       -77.257       -77.257       -77.257      -77.257       -77.257       -77.257       -77.257       -77.257       -77.257       -77.257       -77.257  
    coh13:k=-19.88,b=0        -77.001       -77.001       -77.001       -77.001      -77.001       -77.001       -77.001       -77.001       -77.001       -77.001       -77.001       -77.001  
    coh14:k=-19.87,b=0        -76.774       -76.774       -76.774       -76.774      -76.774       -76.774       -76.774       -76.774       -76.774       -76.774       -76.774       -76.774  
    coh15:k=-19.86,b=0        -76.572       -76.572       -76.572       -76.572      -76.572       -76.572       -76.572       -76.572       -76.572       -76.572       -76.572       -76.572  
    coh16:k=-19.85,b=0        -76.391       -76.391       -76.391       -76.391      -76.391       -76.391       -76.391       -76.391       -76.391       -76.391       -76.391       -76.391  
    coh17:k=-19.84,b=0        -76.227       -76.227       -76.227       -76.227      -76.227       -76.227       -76.227       -76.227       -76.227       -76.227       -76.227       -76.227  
    coh18:k=-19.83,b=0        -76.077       -76.077       -76.077       -76.077      -76.077       -76.077       -76.077       -76.077       -76.077       -76.077       -76.077       -76.077  
    coh19:k=-19.82,b=0        -75.941       -75.941       -75.941       -75.941      -75.941       -75.941       -75.941       -75.941       -75.941       -75.941       -75.941       -75.941  
    coh20:k=-19.81,b=0        -75.815       -75.815       -75.815       -75.815      -75.815       -75.815       -75.815       -75.815       -75.815       -75.815       -75.815       -75.815  
    coh21:k=-19.8,b=0         -75.699       -75.699       -75.699       -75.699      -75.699       -75.699       -75.699       -75.699       -75.699       -75.699       -75.699       -75.699  
    coh22:k=-19.79,b=0        -75.591       -75.591       -75.591       -75.591      -75.591       -75.591       -75.591       -75.591       -75.591       -75.591       -75.591       -75.591  
    coh23:k=-19.78,b=0         -75.49        -75.49        -75.49        -75.49       -75.49        -75.49        -75.49        -75.49        -75.49        -75.49        -75.49        -75.49  
    coh24:k=-19.77,b=0        -75.397       -75.397       -75.397       -75.397      -75.397       -75.397       -75.397       -75.397       -75.397       -75.397       -75.397       -75.397  
    coh25:k=-19.76,b=0        -75.309       -75.309       -75.309       -75.309      -75.309       -75.309       -75.309       -75.309       -75.309       -75.309       -75.309       -75.309  
    coh26:k=-19.75,b=0        -75.226       -75.226       -75.226       -75.226      -75.226       -75.226       -75.226       -75.226       -75.226       -75.226       -75.226       -75.226  
    coh27:k=-19.7399,b=0      -75.149       -75.149       -75.149       -75.149      -75.149       -75.149       -75.149       -75.149       -75.149       -75.149       -75.149       -75.149  
    coh28:k=-19.7299,b=0      -75.075       -75.075       -75.075       -75.075      -75.075       -75.075       -75.075       -75.075       -75.075       -75.075       -75.075       -75.075  
    coh29:k=-19.7199,b=0      -75.006       -75.006       -75.006       -75.006      -75.006       -75.006       -75.006       -75.006       -75.006       -75.006       -75.006       -75.006  
    coh30:k=-19.7099,b=0       -74.94        -74.94        -74.94        -74.94       -74.94        -74.94        -74.94        -74.94        -74.94        -74.94        -74.94        -74.94  
    coh31:k=-19.6999,b=0      -74.878       -74.878       -74.878       -74.878      -74.878       -74.878       -74.878       -74.878       -74.878       -74.878       -74.878       -74.878  
    coh32:k=-19.6899,b=0      -74.818       -74.818       -74.818       -74.818      -74.818       -74.818       -74.818       -74.818       -74.818       -74.818       -74.818       -74.818  
    coh33:k=-19.6799,b=0      -74.762       -74.762       -74.762       -74.762      -74.762       -74.762       -74.762       -74.762       -74.762       -74.762       -74.762       -74.762  
    coh34:k=-19.6699,b=0      -74.708       -74.708       -74.708       -74.708      -74.708       -74.708       -74.708       -74.708       -74.708       -74.708       -74.708       -74.708  
    coh35:k=-19.6599,b=0      -74.656       -74.656       -74.656       -74.656      -74.656       -74.656       -74.656       -74.656       -74.656       -74.656       -74.656       -74.656  
    coh36:k=-19.6499,b=0      -74.607       -74.607       -74.607       -74.607      -74.607       -74.607       -74.607       -74.607       -74.607       -74.607       -74.607       -74.607  
    coh37:k=-19.6399,b=0      -74.559       -74.559       -74.559       -74.559      -74.559       -74.559       -74.559       -74.559       -74.559       -74.559       -74.559       -74.559  
    coh38:k=-19.6299,b=0      -74.514       -74.514       -74.514       -74.514      -74.514       -74.514       -74.514       -74.514       -74.514       -74.514       -74.514       -74.514  
    coh39:k=-19.6199,b=0       -74.47        -74.47        -74.47        -74.47       -74.47        -74.47        -74.47        -74.47        -74.47        -74.47        -74.47        -74.47  
    coh40:k=-19.6099,b=0      -74.429       -74.429       -74.429       -74.429      -74.429       -74.429       -74.429       -74.429       -74.429       -74.429       -74.429       -74.429  
    coh41:k=-19.5999,b=0      -74.388       -74.388       -74.388       -74.388      -74.388       -74.388       -74.388       -74.388       -74.388       -74.388       -74.388       -74.388  
    coh42:k=-19.5899,b=0      -74.349       -74.349       -74.349       -74.349      -74.349       -74.349       -74.349       -74.349       -74.349       -74.349       -74.349       -74.349  
    coh43:k=-19.5799,b=0      -74.312       -74.312       -74.312       -74.312      -74.312       -74.312       -74.312       -74.312       -74.312       -74.312       -74.312       -74.312  
    coh44:k=-19.5699,b=0      -74.276       -74.276       -74.276       -74.276      -74.276       -74.276       -74.276       -74.276       -74.276       -74.276       -74.276       -74.276  
    coh45:k=-19.5599,b=0      -74.241       -74.241       -74.241       -74.241      -74.241       -74.241       -74.241       -74.241       -74.241       -74.241       -74.241       -74.241  
    coh46:k=-19.5499,b=0      -74.207       -74.207       -74.207       -74.207      -74.207       -74.207       -74.207       -74.207       -74.207       -74.207       -74.207       -74.207  
    coh47:k=-19.5399,b=0      -74.175       -74.175       -74.175       -74.175      -74.175       -74.175       -74.175       -74.175       -74.175       -74.175       -74.175       -74.175  
    coh48:k=-19.5299,b=0      -74.143       -74.143       -74.143       -74.143      -74.143       -74.143       -74.143       -74.143       -74.143       -74.143       -74.143       -74.143  
    coh49:k=-19.5199,b=0      -74.113       -74.113       -74.113       -74.113      -74.113       -74.113       -74.113       -74.113       -74.113       -74.113       -74.113       -74.113  
    coh50:k=-19.5099,b=0      -74.083       -74.083       -74.083       -74.083      -74.083       -74.083       -74.083       -74.083       -74.083       -74.083       -74.083       -74.083  
    coh8257:k=62.5765,b=0      10.103        10.103        10.103        10.103       10.103        10.103        10.103        10.103        10.103        10.103        10.103        10.103  
    coh8258:k=62.5865,b=0      10.103        10.103        10.103        10.103       10.103        10.103        10.103        10.103        10.103        10.103        10.103        10.103  
    coh8259:k=62.5965,b=0      10.103        10.103        10.103        10.103       10.103        10.103        10.103        10.103        10.103        10.103        10.103        10.103  
    coh8260:k=62.6065,b=0      10.103        10.103        10.103        10.103       10.103        10.103        10.103        10.103        10.103        10.103        10.103        10.103  
    coh8261:k=62.6165,b=0      10.103        10.103        10.103        10.103       10.103        10.103        10.103        10.103        10.103        10.103        10.103        10.103  
    coh8262:k=62.6265,b=0      10.104        10.104        10.104        10.104       10.104        10.104        10.104        10.104        10.104        10.104        10.104        10.104  
    coh8263:k=62.6365,b=0      10.105        10.105        10.105        10.105       10.105        10.105        10.105        10.105        10.105        10.105        10.105        10.105  
    coh8264:k=62.6465,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8265:k=62.6565,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8266:k=62.6665,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8267:k=62.6765,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8268:k=62.6865,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8269:k=62.6965,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8270:k=62.7065,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8271:k=62.7165,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8272:k=62.7265,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8273:k=62.7365,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8274:k=62.7465,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8275:k=62.7565,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8276:k=62.7665,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8277:k=62.7765,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8278:k=62.7865,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8279:k=62.7965,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8280:k=62.8065,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8281:k=62.8165,b=0      10.107        10.107        10.107        10.107       10.107        10.107        10.107        10.107        10.107        10.107        10.107        10.107  
    coh8282:k=62.8265,b=0      10.109        10.109        10.109        10.109       10.109        10.109        10.109        10.109        10.109        10.109        10.109        10.109  
    coh8283:k=62.8365,b=0       10.11         10.11         10.11         10.11        10.11         10.11         10.11         10.11         10.11         10.11         10.11         10.11  
    coh8284:k=62.8465,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8285:k=62.8565,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8286:k=62.8665,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8287:k=62.8765,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8288:k=62.8865,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8289:k=62.8965,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8290:k=62.9065,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8291:k=62.9165,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8292:k=62.9265,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8293:k=62.9365,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8294:k=62.9465,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8295:k=62.9565,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8296:k=62.9666,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8297:k=62.9766,b=0      10.111        10.111        10.111        10.111       10.111        10.111        10.111        10.111        10.111        10.111        10.111        10.111  
    coh8298:k=62.9866,b=0      10.112        10.112        10.112        10.112       10.112        10.112        10.112        10.112        10.112        10.112        10.112        10.112  
    coh8299:k=62.9966,b=0      10.112        10.112        10.112        10.112       10.112        10.112        10.112        10.112        10.112        10.112        10.112        10.112  
    coh8300:k=63.0066,b=0      10.112        10.112        10.112        10.112       10.112        10.112        10.112        10.112        10.112        10.112        10.112        10.112  
    coh8301:k=63.0166,b=0      10.112        10.112        10.112        10.112       10.112        10.112        10.112        10.112        10.112        10.112        10.112        10.112  
    coh8302:k=63.0266,b=0      10.114        10.114        10.114        10.114       10.114        10.114        10.114        10.114        10.114        10.114        10.114        10.114  
    coh8303:k=63.0366,b=0      10.115        10.115        10.115        10.115       10.115        10.115        10.115        10.115        10.115        10.115        10.115        10.115  
    coh8304:k=63.0466,b=0      10.115        10.115        10.115        10.115       10.115        10.115        10.115        10.115        10.115        10.115        10.115        10.115  
    coh8305:k=63.0566,b=0      10.115        10.115        10.115        10.115       10.115        10.115        10.115        10.115        10.115        10.115        10.115        10.115  
    coh8306:k=63.0666,b=0      10.115        10.115        10.115        10.115       10.115        10.115        10.115        10.115        10.115        10.115        10.115        10.115  

tb_pol_a
                             z1_0_85962    z2_0_87824    z3_0_89726    z4_0_9167    z5_0_93655    z6_0_95684    z10_1_0425    z11_1_0651    z12_1_0881    z13_1_1117    z14_1_1358    z15_1_1604
                             __________    __________    __________    _________    __________    __________    __________    __________    __________    __________    __________    __________

    coh1:k=-20,b=0              -9.14         -9.14         -9.14        -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14  
    coh2:k=-19.99,b=0           -9.14         -9.14         -9.14        -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14  
    coh3:k=-19.98,b=0           -9.14         -9.14         -9.14        -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14  
    coh4:k=-19.97,b=0           -9.14         -9.14         -9.14        -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14  
    coh5:k=-19.96,b=0           -9.14         -9.14         -9.14        -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14  
    coh6:k=-19.95,b=0           -9.14         -9.14         -9.14        -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14  
    coh7:k=-19.94,b=0           -9.14         -9.14         -9.14        -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14  
    coh8:k=-19.93,b=0           -9.14         -9.14         -9.14        -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14  
    coh9:k=-19.92,b=0           -9.14         -9.14         -9.14        -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14         -9.14  
    coh10:k=-19.91,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh11:k=-19.9,b=0             -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh12:k=-19.89,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh13:k=-19.88,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh14:k=-19.87,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh15:k=-19.86,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh16:k=-19.85,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh17:k=-19.84,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh18:k=-19.83,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh19:k=-19.82,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh20:k=-19.81,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh21:k=-19.8,b=0             -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh22:k=-19.79,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh23:k=-19.78,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh24:k=-19.77,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh25:k=-19.76,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh26:k=-19.75,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh27:k=-19.7399,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh28:k=-19.7299,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh29:k=-19.7199,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh30:k=-19.7099,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh31:k=-19.6999,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh32:k=-19.6899,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh33:k=-19.6799,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh34:k=-19.6699,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh35:k=-19.6599,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh36:k=-19.6499,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh37:k=-19.6399,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh38:k=-19.6299,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh39:k=-19.6199,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh40:k=-19.6099,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh41:k=-19.5999,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh42:k=-19.5899,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh43:k=-19.5799,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh44:k=-19.5699,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh45:k=-19.5599,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh46:k=-19.5499,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh47:k=-19.5399,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh48:k=-19.5299,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh49:k=-19.5199,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh50:k=-19.5099,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh8257:k=62.5765,b=0      34.217        34.217        34.217       34.217        34.217        34.217        34.217        34.217        34.217        34.217        34.217        34.217  
    coh8258:k=62.5865,b=0      34.226        34.226        34.226       34.226        34.226        34.226        34.226        34.226        34.226        34.226        34.226        34.226  
    coh8259:k=62.5965,b=0      34.234        34.234        34.234       34.234        34.234        34.234        34.234        34.234        34.234        34.234        34.234        34.234  
    coh8260:k=62.6065,b=0      34.243        34.243        34.243       34.243        34.243        34.243        34.243        34.243        34.243        34.243        34.243        34.243  
    coh8261:k=62.6165,b=0      34.251        34.251        34.251       34.251        34.251        34.251        34.251        34.251        34.251        34.251        34.251        34.251  
    coh8262:k=62.6265,b=0      34.023        34.023        34.023       34.023        34.023        34.023        34.023        34.023        34.023        34.023        34.023        34.023  
    coh8263:k=62.6365,b=0      34.105        34.105        34.105       34.105        34.105        34.105        34.105        34.105        34.105        34.105        34.105        34.105  
    coh8264:k=62.6465,b=0      34.112        34.112        34.112       34.112        34.112        34.112        34.112        34.112        34.112        34.112        34.112        34.112  
    coh8265:k=62.6565,b=0       34.12         34.12         34.12        34.12         34.12         34.12         34.12         34.12         34.12         34.12         34.12         34.12  
    coh8266:k=62.6665,b=0      34.129        34.129        34.129       34.129        34.129        34.129        34.129        34.129        34.129        34.129        34.129        34.129  
    coh8267:k=62.6765,b=0      34.137        34.137        34.137       34.137        34.137        34.137        34.137        34.137        34.137        34.137        34.137        34.137  
    coh8268:k=62.6865,b=0      34.145        34.145        34.145       34.145        34.145        34.145        34.145        34.145        34.145        34.145        34.145        34.145  
    coh8269:k=62.6965,b=0      34.154        34.154        34.154       34.154        34.154        34.154        34.154        34.154        34.154        34.154        34.154        34.154  
    coh8270:k=62.7065,b=0      34.162        34.162        34.162       34.162        34.162        34.162        34.162        34.162        34.162        34.162        34.162        34.162  
    coh8271:k=62.7165,b=0      34.171        34.171        34.171       34.171        34.171        34.171        34.171        34.171        34.171        34.171        34.171        34.171  
    coh8272:k=62.7265,b=0      34.179        34.179        34.179       34.179        34.179        34.179        34.179        34.179        34.179        34.179        34.179        34.179  
    coh8273:k=62.7365,b=0      34.188        34.188        34.188       34.188        34.188        34.188        34.188        34.188        34.188        34.188        34.188        34.188  
    coh8274:k=62.7465,b=0      34.196        34.196        34.196       34.196        34.196        34.196        34.196        34.196        34.196        34.196        34.196        34.196  
    coh8275:k=62.7565,b=0      34.205        34.205        34.205       34.205        34.205        34.205        34.205        34.205        34.205        34.205        34.205        34.205  
    coh8276:k=62.7665,b=0      34.213        34.213        34.213       34.213        34.213        34.213        34.213        34.213        34.213        34.213        34.213        34.213  
    coh8277:k=62.7765,b=0      34.222        34.222        34.222       34.222        34.222        34.222        34.222        34.222        34.222        34.222        34.222        34.222  
    coh8278:k=62.7865,b=0       34.23         34.23         34.23        34.23         34.23         34.23         34.23         34.23         34.23         34.23         34.23         34.23  
    coh8279:k=62.7965,b=0      34.239        34.239        34.239       34.239        34.239        34.239        34.239        34.239        34.239        34.239        34.239        34.239  
    coh8280:k=62.8065,b=0      34.247        34.247        34.247       34.247        34.247        34.247        34.247        34.247        34.247        34.247        34.247        34.247  
    coh8281:k=62.8165,b=0      34.256        34.256        34.256       34.256        34.256        34.256        34.256        34.256        34.256        34.256        34.256        34.256  
    coh8282:k=62.8265,b=0      34.079        34.079        34.079       34.079        34.079        34.079        34.079        34.079        34.079        34.079        34.079        34.079  
    coh8283:k=62.8365,b=0      34.108        34.108        34.108       34.108        34.108        34.108        34.108        34.108        34.108        34.108        34.108        34.108  
    coh8284:k=62.8465,b=0      34.115        34.115        34.115       34.115        34.115        34.115        34.115        34.115        34.115        34.115        34.115        34.115  
    coh8285:k=62.8565,b=0      34.124        34.124        34.124       34.124        34.124        34.124        34.124        34.124        34.124        34.124        34.124        34.124  
    coh8286:k=62.8665,b=0      34.132        34.132        34.132       34.132        34.132        34.132        34.132        34.132        34.132        34.132        34.132        34.132  
    coh8287:k=62.8765,b=0      34.141        34.141        34.141       34.141        34.141        34.141        34.141        34.141        34.141        34.141        34.141        34.141  
    coh8288:k=62.8865,b=0      34.149        34.149        34.149       34.149        34.149        34.149        34.149        34.149        34.149        34.149        34.149        34.149  
    coh8289:k=62.8965,b=0      34.158        34.158        34.158       34.158        34.158        34.158        34.158        34.158        34.158        34.158        34.158        34.158  
    coh8290:k=62.9065,b=0      34.166        34.166        34.166       34.166        34.166        34.166        34.166        34.166        34.166        34.166        34.166        34.166  
    coh8291:k=62.9165,b=0      34.174        34.174        34.174       34.174        34.174        34.174        34.174        34.174        34.174        34.174        34.174        34.174  
    coh8292:k=62.9265,b=0      34.183        34.183        34.183       34.183        34.183        34.183        34.183        34.183        34.183        34.183        34.183        34.183  
    coh8293:k=62.9365,b=0      34.191        34.191        34.191       34.191        34.191        34.191        34.191        34.191        34.191        34.191        34.191        34.191  
    coh8294:k=62.9465,b=0        34.2          34.2          34.2         34.2          34.2          34.2          34.2          34.2          34.2          34.2          34.2          34.2  
    coh8295:k=62.9565,b=0      34.208        34.208        34.208       34.208        34.208        34.208        34.208        34.208        34.208        34.208        34.208        34.208  
    coh8296:k=62.9666,b=0      34.217        34.217        34.217       34.217        34.217        34.217        34.217        34.217        34.217        34.217        34.217        34.217  
    coh8297:k=62.9766,b=0      34.225        34.225        34.225       34.225        34.225        34.225        34.225        34.225        34.225        34.225        34.225        34.225  
    coh8298:k=62.9866,b=0      34.234        34.234        34.234       34.234        34.234        34.234        34.234        34.234        34.234        34.234        34.234        34.234  
    coh8299:k=62.9966,b=0      34.242        34.242        34.242       34.242        34.242        34.242        34.242        34.242        34.242        34.242        34.242        34.242  
    coh8300:k=63.0066,b=0      34.251        34.251        34.251       34.251        34.251        34.251        34.251        34.251        34.251        34.251        34.251        34.251  
    coh8301:k=63.0166,b=0      33.998        33.998        33.998       33.998        33.998        33.998        33.998        33.998        33.998        33.998        33.998        33.998  
    coh8302:k=63.0266,b=0      34.103        34.103        34.103       34.103        34.103        34.103        34.103        34.103        34.103        34.103        34.103        34.103  
    coh8303:k=63.0366,b=0       34.11         34.11         34.11        34.11         34.11         34.11         34.11         34.11         34.11         34.11         34.11         34.11  
    coh8304:k=63.0466,b=0      34.118        34.118        34.118       34.118        34.118        34.118        34.118        34.118        34.118        34.118        34.118        34.118  
    coh8305:k=63.0566,b=0      34.127        34.127        34.127       34.127        34.127        34.127        34.127        34.127        34.127        34.127        34.127        34.127  
    coh8306:k=63.0666,b=0      34.135        34.135        34.135       34.135        34.135        34.135        34.135        34.135        34.135        34.135        34.135        34.135  

tb_pol_k
                             z1_0_85962    z2_0_87824    z3_0_89726    z4_0_9167     z5_0_93655    z6_0_95684    z10_1_0425    z11_1_0651    z12_1_0881    z13_1_1117    z14_1_1358    z15_1_1604
                             __________    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________    __________

    coh1:k=-20,b=0                 9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14
    coh2:k=-19.99,b=0              9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14
    coh3:k=-19.98,b=0              9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14
    coh4:k=-19.97,b=0              9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14
    coh5:k=-19.96,b=0              9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14
    coh6:k=-19.95,b=0              9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14
    coh7:k=-19.94,b=0              9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14
    coh8:k=-19.93,b=0              9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14
    coh9:k=-19.92,b=0              9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14          9.14
    coh10:k=-19.91,b=0       9.0018e-08    9.0018e-08    9.0018e-08    9.0018e-08    9.0018e-08    9.0018e-08    9.0018e-08    9.0018e-08    9.0018e-08    9.0018e-08    9.0018e-08    9.0018e-08
    coh11:k=-19.9,b=0        1.0002e-07    1.0002e-07    1.0002e-07    1.0002e-07    1.0002e-07    1.0002e-07    1.0002e-07    1.0002e-07    1.0002e-07    1.0002e-07    1.0002e-07    1.0002e-07
    coh12:k=-19.89,b=0       1.1002e-07    1.1002e-07    1.1002e-07    1.1002e-07    1.1002e-07    1.1002e-07    1.1002e-07    1.1002e-07    1.1002e-07    1.1002e-07    1.1002e-07    1.1002e-07
    coh13:k=-19.88,b=0       1.2002e-07    1.2002e-07    1.2002e-07    1.2002e-07    1.2002e-07    1.2002e-07    1.2002e-07    1.2002e-07    1.2002e-07    1.2002e-07    1.2002e-07    1.2002e-07
    coh14:k=-19.87,b=0       1.3003e-07    1.3003e-07    1.3003e-07    1.3003e-07    1.3003e-07    1.3003e-07    1.3003e-07    1.3003e-07    1.3003e-07    1.3003e-07    1.3003e-07    1.3003e-07
    coh15:k=-19.86,b=0       1.4003e-07    1.4003e-07    1.4003e-07    1.4003e-07    1.4003e-07    1.4003e-07    1.4003e-07    1.4003e-07    1.4003e-07    1.4003e-07    1.4003e-07    1.4003e-07
    coh16:k=-19.85,b=0       1.5003e-07    1.5003e-07    1.5003e-07    1.5003e-07    1.5003e-07    1.5003e-07    1.5003e-07    1.5003e-07    1.5003e-07    1.5003e-07    1.5003e-07    1.5003e-07
    coh17:k=-19.84,b=0       1.6003e-07    1.6003e-07    1.6003e-07    1.6003e-07    1.6003e-07    1.6003e-07    1.6003e-07    1.6003e-07    1.6003e-07    1.6003e-07    1.6003e-07    1.6003e-07
    coh18:k=-19.83,b=0       1.7003e-07    1.7003e-07    1.7003e-07    1.7003e-07    1.7003e-07    1.7003e-07    1.7003e-07    1.7003e-07    1.7003e-07    1.7003e-07    1.7003e-07    1.7003e-07
    coh19:k=-19.82,b=0       1.8004e-07    1.8004e-07    1.8004e-07    1.8004e-07    1.8004e-07    1.8004e-07    1.8004e-07    1.8004e-07    1.8004e-07    1.8004e-07    1.8004e-07    1.8004e-07
    coh20:k=-19.81,b=0       1.9004e-07    1.9004e-07    1.9004e-07    1.9004e-07    1.9004e-07    1.9004e-07    1.9004e-07    1.9004e-07    1.9004e-07    1.9004e-07    1.9004e-07    1.9004e-07
    coh21:k=-19.8,b=0        2.0004e-07    2.0004e-07    2.0004e-07    2.0004e-07    2.0004e-07    2.0004e-07    2.0004e-07    2.0004e-07    2.0004e-07    2.0004e-07    2.0004e-07    2.0004e-07
    coh22:k=-19.79,b=0       2.1004e-07    2.1004e-07    2.1004e-07    2.1004e-07    2.1004e-07    2.1004e-07    2.1004e-07    2.1004e-07    2.1004e-07    2.1004e-07    2.1004e-07    2.1004e-07
    coh23:k=-19.78,b=0       2.2004e-07    2.2004e-07    2.2004e-07    2.2004e-07    2.2004e-07    2.2004e-07    2.2004e-07    2.2004e-07    2.2004e-07    2.2004e-07    2.2004e-07    2.2004e-07
    coh24:k=-19.77,b=0       2.3005e-07    2.3005e-07    2.3005e-07    2.3005e-07    2.3005e-07    2.3005e-07    2.3005e-07    2.3005e-07    2.3005e-07    2.3005e-07    2.3005e-07    2.3005e-07
    coh25:k=-19.76,b=0       2.4005e-07    2.4005e-07    2.4005e-07    2.4005e-07    2.4005e-07    2.4005e-07    2.4005e-07    2.4005e-07    2.4005e-07    2.4005e-07    2.4005e-07    2.4005e-07
    coh26:k=-19.75,b=0       2.5005e-07    2.5005e-07    2.5005e-07    2.5005e-07    2.5005e-07    2.5005e-07    2.5005e-07    2.5005e-07    2.5005e-07    2.5005e-07    2.5005e-07    2.5005e-07
    coh27:k=-19.7399,b=0     2.6005e-07    2.6005e-07    2.6005e-07    2.6005e-07    2.6005e-07    2.6005e-07    2.6005e-07    2.6005e-07    2.6005e-07    2.6005e-07    2.6005e-07    2.6005e-07
    coh28:k=-19.7299,b=0     2.7005e-07    2.7005e-07    2.7005e-07    2.7005e-07    2.7005e-07    2.7005e-07    2.7005e-07    2.7005e-07    2.7005e-07    2.7005e-07    2.7005e-07    2.7005e-07
    coh29:k=-19.7199,b=0     2.8006e-07    2.8006e-07    2.8006e-07    2.8006e-07    2.8006e-07    2.8006e-07    2.8006e-07    2.8006e-07    2.8006e-07    2.8006e-07    2.8006e-07    2.8006e-07
    coh30:k=-19.7099,b=0     2.9006e-07    2.9006e-07    2.9006e-07    2.9006e-07    2.9006e-07    2.9006e-07    2.9006e-07    2.9006e-07    2.9006e-07    2.9006e-07    2.9006e-07    2.9006e-07
    coh31:k=-19.6999,b=0     3.0006e-07    3.0006e-07    3.0006e-07    3.0006e-07    3.0006e-07    3.0006e-07    3.0006e-07    3.0006e-07    3.0006e-07    3.0006e-07    3.0006e-07    3.0006e-07
    coh32:k=-19.6899,b=0     3.1006e-07    3.1006e-07    3.1006e-07    3.1006e-07    3.1006e-07    3.1006e-07    3.1006e-07    3.1006e-07    3.1006e-07    3.1006e-07    3.1006e-07    3.1006e-07
    coh33:k=-19.6799,b=0     3.2006e-07    3.2006e-07    3.2006e-07    3.2006e-07    3.2006e-07    3.2006e-07    3.2006e-07    3.2006e-07    3.2006e-07    3.2006e-07    3.2006e-07    3.2006e-07
    coh34:k=-19.6699,b=0     3.3007e-07    3.3007e-07    3.3007e-07    3.3007e-07    3.3007e-07    3.3007e-07    3.3007e-07    3.3007e-07    3.3007e-07    3.3007e-07    3.3007e-07    3.3007e-07
    coh35:k=-19.6599,b=0     3.4007e-07    3.4007e-07    3.4007e-07    3.4007e-07    3.4007e-07    3.4007e-07    3.4007e-07    3.4007e-07    3.4007e-07    3.4007e-07    3.4007e-07    3.4007e-07
    coh36:k=-19.6499,b=0     3.5007e-07    3.5007e-07    3.5007e-07    3.5007e-07    3.5007e-07    3.5007e-07    3.5007e-07    3.5007e-07    3.5007e-07    3.5007e-07    3.5007e-07    3.5007e-07
    coh37:k=-19.6399,b=0     3.6007e-07    3.6007e-07    3.6007e-07    3.6007e-07    3.6007e-07    3.6007e-07    3.6007e-07    3.6007e-07    3.6007e-07    3.6007e-07    3.6007e-07    3.6007e-07
    coh38:k=-19.6299,b=0     3.7007e-07    3.7007e-07    3.7007e-07    3.7007e-07    3.7007e-07    3.7007e-07    3.7007e-07    3.7007e-07    3.7007e-07    3.7007e-07    3.7007e-07    3.7007e-07
    coh39:k=-19.6199,b=0     3.8008e-07    3.8008e-07    3.8008e-07    3.8008e-07    3.8008e-07    3.8008e-07    3.8008e-07    3.8008e-07    3.8008e-07    3.8008e-07    3.8008e-07    3.8008e-07
    coh40:k=-19.6099,b=0     3.9008e-07    3.9008e-07    3.9008e-07    3.9008e-07    3.9008e-07    3.9008e-07    3.9008e-07    3.9008e-07    3.9008e-07    3.9008e-07    3.9008e-07    3.9008e-07
    coh41:k=-19.5999,b=0     4.0008e-07    4.0008e-07    4.0008e-07    4.0008e-07    4.0008e-07    4.0008e-07    4.0008e-07    4.0008e-07    4.0008e-07    4.0008e-07    4.0008e-07    4.0008e-07
    coh42:k=-19.5899,b=0     4.1008e-07    4.1008e-07    4.1008e-07    4.1008e-07    4.1008e-07    4.1008e-07    4.1008e-07    4.1008e-07    4.1008e-07    4.1008e-07    4.1008e-07    4.1008e-07
    coh43:k=-19.5799,b=0     4.2008e-07    4.2008e-07    4.2008e-07    4.2008e-07    4.2008e-07    4.2008e-07    4.2008e-07    4.2008e-07    4.2008e-07    4.2008e-07    4.2008e-07    4.2008e-07
    coh44:k=-19.5699,b=0     4.3009e-07    4.3009e-07    4.3009e-07    4.3009e-07    4.3009e-07    4.3009e-07    4.3009e-07    4.3009e-07    4.3009e-07    4.3009e-07    4.3009e-07    4.3009e-07
    coh45:k=-19.5599,b=0     4.4009e-07    4.4009e-07    4.4009e-07    4.4009e-07    4.4009e-07    4.4009e-07    4.4009e-07    4.4009e-07    4.4009e-07    4.4009e-07    4.4009e-07    4.4009e-07
    coh46:k=-19.5499,b=0     4.5009e-07    4.5009e-07    4.5009e-07    4.5009e-07    4.5009e-07    4.5009e-07    4.5009e-07    4.5009e-07    4.5009e-07    4.5009e-07    4.5009e-07    4.5009e-07
    coh47:k=-19.5399,b=0     4.6009e-07    4.6009e-07    4.6009e-07    4.6009e-07    4.6009e-07    4.6009e-07    4.6009e-07    4.6009e-07    4.6009e-07    4.6009e-07    4.6009e-07    4.6009e-07
    coh48:k=-19.5299,b=0     4.7009e-07    4.7009e-07    4.7009e-07    4.7009e-07    4.7009e-07    4.7009e-07    4.7009e-07    4.7009e-07    4.7009e-07    4.7009e-07    4.7009e-07    4.7009e-07
    coh49:k=-19.5199,b=0      4.801e-07     4.801e-07     4.801e-07     4.801e-07     4.801e-07     4.801e-07     4.801e-07     4.801e-07     4.801e-07     4.801e-07     4.801e-07     4.801e-07
    coh50:k=-19.5099,b=0      4.901e-07     4.901e-07     4.901e-07     4.901e-07     4.901e-07     4.901e-07     4.901e-07     4.901e-07     4.901e-07     4.901e-07     4.901e-07     4.901e-07
    coh8257:k=62.5765,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8258:k=62.5865,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8259:k=62.5965,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8260:k=62.6065,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8261:k=62.6165,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8262:k=62.6265,b=0        15.962        15.962        15.962        15.962        15.962        15.962        15.962        15.962        15.962        15.962        15.962        15.962
    coh8263:k=62.6365,b=0        15.888        15.888        15.888        15.888        15.888        15.888        15.888        15.888        15.888        15.888        15.888        15.888
    coh8264:k=62.6465,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8265:k=62.6565,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8266:k=62.6665,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8267:k=62.6765,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8268:k=62.6865,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8269:k=62.6965,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8270:k=62.7065,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8271:k=62.7165,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8272:k=62.7265,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8273:k=62.7365,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8274:k=62.7465,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8275:k=62.7565,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8276:k=62.7665,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8277:k=62.7765,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8278:k=62.7865,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8279:k=62.7965,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8280:k=62.8065,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8281:k=62.8165,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8282:k=62.8265,b=0         15.91         15.91         15.91         15.91         15.91         15.91         15.91         15.91         15.91         15.91         15.91         15.91
    coh8283:k=62.8365,b=0        15.889        15.889        15.889        15.889        15.889        15.889        15.889        15.889        15.889        15.889        15.889        15.889
    coh8284:k=62.8465,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8285:k=62.8565,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8286:k=62.8665,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8287:k=62.8765,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8288:k=62.8865,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8289:k=62.8965,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8290:k=62.9065,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8291:k=62.9165,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8292:k=62.9265,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8293:k=62.9365,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8294:k=62.9465,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8295:k=62.9565,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8296:k=62.9666,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8297:k=62.9766,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8298:k=62.9866,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8299:k=62.9966,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8300:k=63.0066,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8301:k=63.0166,b=0        15.985        15.985        15.985        15.985        15.985        15.985        15.985        15.985        15.985        15.985        15.985        15.985
    coh8302:k=63.0266,b=0        15.888        15.888        15.888        15.888        15.888        15.888        15.888        15.888        15.888        15.888        15.888        15.888
    coh8303:k=63.0366,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8304:k=63.0466,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8305:k=63.0566,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89
    coh8306:k=63.0666,b=0         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89         15.89

tb_pol_w
                             z1_0_85962    z2_0_87824    z3_0_89726    z4_0_9167    z5_0_93655    z6_0_95684    z10_1_0425    z11_1_0651    z12_1_0881    z13_1_1117    z14_1_1358    z15_1_1604
                             __________    __________    __________    _________    __________    __________    __________    __________    __________    __________    __________    __________

    coh1:k=-20,b=0                  0             0             0            0             0             0             0             0             0             0             0             0  
    coh2:k=-19.99,b=0               0             0             0            0             0             0             0             0             0             0             0             0  
    coh3:k=-19.98,b=0               0             0             0            0             0             0             0             0             0             0             0             0  
    coh4:k=-19.97,b=0               0             0             0            0             0             0             0             0             0             0             0             0  
    coh5:k=-19.96,b=0               0             0             0            0             0             0             0             0             0             0             0             0  
    coh6:k=-19.95,b=0               0             0             0            0             0             0             0             0             0             0             0             0  
    coh7:k=-19.94,b=0               0             0             0            0             0             0             0             0             0             0             0             0  
    coh8:k=-19.93,b=0               0             0             0            0             0             0             0             0             0             0             0             0  
    coh9:k=-19.92,b=0               0             0             0            0             0             0             0             0             0             0             0             0  
    coh10:k=-19.91,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh11:k=-19.9,b=0             -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh12:k=-19.89,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh13:k=-19.88,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh14:k=-19.87,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh15:k=-19.86,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh16:k=-19.85,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh17:k=-19.84,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh18:k=-19.83,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh19:k=-19.82,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh20:k=-19.81,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh21:k=-19.8,b=0             -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh22:k=-19.79,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh23:k=-19.78,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh24:k=-19.77,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh25:k=-19.76,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh26:k=-19.75,b=0            -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh27:k=-19.7399,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh28:k=-19.7299,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh29:k=-19.7199,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh30:k=-19.7099,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh31:k=-19.6999,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh32:k=-19.6899,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh33:k=-19.6799,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh34:k=-19.6699,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh35:k=-19.6599,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh36:k=-19.6499,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh37:k=-19.6399,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh38:k=-19.6299,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh39:k=-19.6199,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh40:k=-19.6099,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh41:k=-19.5999,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh42:k=-19.5899,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh43:k=-19.5799,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh44:k=-19.5699,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh45:k=-19.5599,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh46:k=-19.5499,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh47:k=-19.5399,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh48:k=-19.5299,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh49:k=-19.5199,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh50:k=-19.5099,b=0          -20           -20           -20          -20           -20           -20           -20           -20           -20           -20           -20           -20  
    coh8257:k=62.5765,b=0      50.107        50.107        50.107       50.107        50.107        50.107        50.107        50.107        50.107        50.107        50.107        50.107  
    coh8258:k=62.5865,b=0      50.116        50.116        50.116       50.116        50.116        50.116        50.116        50.116        50.116        50.116        50.116        50.116  
    coh8259:k=62.5965,b=0      50.124        50.124        50.124       50.124        50.124        50.124        50.124        50.124        50.124        50.124        50.124        50.124  
    coh8260:k=62.6065,b=0      50.133        50.133        50.133       50.133        50.133        50.133        50.133        50.133        50.133        50.133        50.133        50.133  
    coh8261:k=62.6165,b=0      50.141        50.141        50.141       50.141        50.141        50.141        50.141        50.141        50.141        50.141        50.141        50.141  
    coh8262:k=62.6265,b=0      49.985        49.985        49.985       49.985        49.985        49.985        49.985        49.985        49.985        49.985        49.985        49.985  
    coh8263:k=62.6365,b=0      49.993        49.993        49.993       49.993        49.993        49.993        49.993        49.993        49.993        49.993        49.993        49.993  
    coh8264:k=62.6465,b=0      50.002        50.002        50.002       50.002        50.002        50.002        50.002        50.002        50.002        50.002        50.002        50.002  
    coh8265:k=62.6565,b=0       50.01         50.01         50.01        50.01         50.01         50.01         50.01         50.01         50.01         50.01         50.01         50.01  
    coh8266:k=62.6665,b=0      50.019        50.019        50.019       50.019        50.019        50.019        50.019        50.019        50.019        50.019        50.019        50.019  
    coh8267:k=62.6765,b=0      50.027        50.027        50.027       50.027        50.027        50.027        50.027        50.027        50.027        50.027        50.027        50.027  
    coh8268:k=62.6865,b=0      50.035        50.035        50.035       50.035        50.035        50.035        50.035        50.035        50.035        50.035        50.035        50.035  
    coh8269:k=62.6965,b=0      50.044        50.044        50.044       50.044        50.044        50.044        50.044        50.044        50.044        50.044        50.044        50.044  
    coh8270:k=62.7065,b=0      50.052        50.052        50.052       50.052        50.052        50.052        50.052        50.052        50.052        50.052        50.052        50.052  
    coh8271:k=62.7165,b=0      50.061        50.061        50.061       50.061        50.061        50.061        50.061        50.061        50.061        50.061        50.061        50.061  
    coh8272:k=62.7265,b=0      50.069        50.069        50.069       50.069        50.069        50.069        50.069        50.069        50.069        50.069        50.069        50.069  
    coh8273:k=62.7365,b=0      50.078        50.078        50.078       50.078        50.078        50.078        50.078        50.078        50.078        50.078        50.078        50.078  
    coh8274:k=62.7465,b=0      50.086        50.086        50.086       50.086        50.086        50.086        50.086        50.086        50.086        50.086        50.086        50.086  
    coh8275:k=62.7565,b=0      50.095        50.095        50.095       50.095        50.095        50.095        50.095        50.095        50.095        50.095        50.095        50.095  
    coh8276:k=62.7665,b=0      50.103        50.103        50.103       50.103        50.103        50.103        50.103        50.103        50.103        50.103        50.103        50.103  
    coh8277:k=62.7765,b=0      50.112        50.112        50.112       50.112        50.112        50.112        50.112        50.112        50.112        50.112        50.112        50.112  
    coh8278:k=62.7865,b=0       50.12         50.12         50.12        50.12         50.12         50.12         50.12         50.12         50.12         50.12         50.12         50.12  
    coh8279:k=62.7965,b=0      50.129        50.129        50.129       50.129        50.129        50.129        50.129        50.129        50.129        50.129        50.129        50.129  
    coh8280:k=62.8065,b=0      50.137        50.137        50.137       50.137        50.137        50.137        50.137        50.137        50.137        50.137        50.137        50.137  
    coh8281:k=62.8165,b=0      50.146        50.146        50.146       50.146        50.146        50.146        50.146        50.146        50.146        50.146        50.146        50.146  
    coh8282:k=62.8265,b=0      49.988        49.988        49.988       49.988        49.988        49.988        49.988        49.988        49.988        49.988        49.988        49.988  
    coh8283:k=62.8365,b=0      49.997        49.997        49.997       49.997        49.997        49.997        49.997        49.997        49.997        49.997        49.997        49.997  
    coh8284:k=62.8465,b=0      50.005        50.005        50.005       50.005        50.005        50.005        50.005        50.005        50.005        50.005        50.005        50.005  
    coh8285:k=62.8565,b=0      50.014        50.014        50.014       50.014        50.014        50.014        50.014        50.014        50.014        50.014        50.014        50.014  
    coh8286:k=62.8665,b=0      50.022        50.022        50.022       50.022        50.022        50.022        50.022        50.022        50.022        50.022        50.022        50.022  
    coh8287:k=62.8765,b=0      50.031        50.031        50.031       50.031        50.031        50.031        50.031        50.031        50.031        50.031        50.031        50.031  
    coh8288:k=62.8865,b=0      50.039        50.039        50.039       50.039        50.039        50.039        50.039        50.039        50.039        50.039        50.039        50.039  
    coh8289:k=62.8965,b=0      50.048        50.048        50.048       50.048        50.048        50.048        50.048        50.048        50.048        50.048        50.048        50.048  
    coh8290:k=62.9065,b=0      50.056        50.056        50.056       50.056        50.056        50.056        50.056        50.056        50.056        50.056        50.056        50.056  
    coh8291:k=62.9165,b=0      50.064        50.064        50.064       50.064        50.064        50.064        50.064        50.064        50.064        50.064        50.064        50.064  
    coh8292:k=62.9265,b=0      50.073        50.073        50.073       50.073        50.073        50.073        50.073        50.073        50.073        50.073        50.073        50.073  
    coh8293:k=62.9365,b=0      50.081        50.081        50.081       50.081        50.081        50.081        50.081        50.081        50.081        50.081        50.081        50.081  
    coh8294:k=62.9465,b=0       50.09         50.09         50.09        50.09         50.09         50.09         50.09         50.09         50.09         50.09         50.09         50.09  
    coh8295:k=62.9565,b=0      50.098        50.098        50.098       50.098        50.098        50.098        50.098        50.098        50.098        50.098        50.098        50.098  
    coh8296:k=62.9666,b=0      50.107        50.107        50.107       50.107        50.107        50.107        50.107        50.107        50.107        50.107        50.107        50.107  
    coh8297:k=62.9766,b=0      50.115        50.115        50.115       50.115        50.115        50.115        50.115        50.115        50.115        50.115        50.115        50.115  
    coh8298:k=62.9866,b=0      50.124        50.124        50.124       50.124        50.124        50.124        50.124        50.124        50.124        50.124        50.124        50.124  
    coh8299:k=62.9966,b=0      50.132        50.132        50.132       50.132        50.132        50.132        50.132        50.132        50.132        50.132        50.132        50.132  
    coh8300:k=63.0066,b=0      50.141        50.141        50.141       50.141        50.141        50.141        50.141        50.141        50.141        50.141        50.141        50.141  
    coh8301:k=63.0166,b=0      49.983        49.983        49.983       49.983        49.983        49.983        49.983        49.983        49.983        49.983        49.983        49.983  
    coh8302:k=63.0266,b=0      49.991        49.991        49.991       49.991        49.991        49.991        49.991        49.991        49.991        49.991        49.991        49.991  
    coh8303:k=63.0366,b=0          50            50            50           50            50            50            50            50            50            50            50            50  
    coh8304:k=63.0466,b=0      50.008        50.008        50.008       50.008        50.008        50.008        50.008        50.008        50.008        50.008        50.008        50.008  
    coh8305:k=63.0566,b=0      50.017        50.017        50.017       50.017        50.017        50.017        50.017        50.017        50.017        50.017        50.017        50.017  
    coh8306:k=63.0666,b=0      50.025        50.025        50.025       50.025        50.025        50.025        50.025        50.025        50.025        50.025        50.025        50.025  

Simulate Risky Asset with CRS and full depreciation, Persistent Shocks

Now Shocks matter conditional on coh(z) still.

it_param_set = 4;
[param_map, support_map] = ffs_ipwkbz_set_default_param(it_param_set);

% Simulation Accuracy
param_map('it_w_perc_n') = 500;
param_map('it_ak_perc_n') = param_map('it_w_perc_n');
param_map('it_z_n') = 15;

param_map('fl_coh_interp_grid_gap') = 0.01;
param_map('fl_w_interp_grid_gap') = 0.01;
param_map('it_c_interp_grid_gap') = 10^-4;

% Turn on 2nd stage graphs
support_map('bl_graph_evf') = true;
support_map('bl_display_evf') = false;

% Production Function Parameters
% note shock is log normal
param_map('fl_Amean') = 1.0265;
param_map('fl_alpha') = 1;
param_map('fl_delta') = 1;
param_map('fl_r') = 0.03;
param_map('fl_w') = 0.20;

% Shock Parameter, iid shocks
param_map('fl_z_rho') = 0.05;
param_map('fl_z_sig') = 0.05;

% Display Parameters
support_map('bl_display') = false;
support_map('bl_display_final') = false;
support_map('bl_time') = true;
% support_map('bl_profile') = false;

% Call Program
ff_ipwkbz_vf_vecsv(param_map, support_map);
Warning: Ignoring extra legend entries. 
Elapsed time is 423.821679 seconds.
Warning: Directory already exists.