Contents

IPWKZ 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_ipwkz_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_ipwkz_vf_vecsv(param_map, support_map);
Warning: Ignoring extra legend entries. 
Elapsed time is 273.519580 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        465.46      606.7             1             1             1             1             1             1             1             1             1             1             1             1             1             1             1
    iter=2        355.22     4995.9             1             1             1             1             1             1             1             1             1             1             1             1             1             1             1
    iter=3        292.33     1669.6        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995
    iter=4        248.17     836.83        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995        0.9995
    iter=5        215.31     499.22       0.99882       0.99882       0.99882       0.99882       0.99882       0.99882       0.99882       0.99882       0.99882       0.99882       0.99882       0.99882       0.99882       0.99882       0.99882
    iter=6        190.18     332.05       0.99832       0.99832       0.99832       0.99832       0.99832       0.99832       0.99832       0.99832       0.99832       0.99832       0.99832       0.99832       0.99832       0.99832       0.99832
    iter=7        170.67      233.4       0.99798       0.99798       0.99798       0.99798       0.99798       0.99798       0.99798       0.99798       0.99798       0.99798       0.99798       0.99798       0.99798       0.99798       0.99798
    iter=8        155.34      180.3       0.99714       0.99714       0.99714       0.99714       0.99714       0.99714       0.99714       0.99714       0.99714       0.99714       0.99714       0.99714       0.99714       0.99714       0.99714
    iter=9        143.14     133.51       0.99631       0.99631       0.99631       0.99631       0.99631       0.99631       0.99631       0.99631       0.99631       0.99631       0.99631       0.99631       0.99631       0.99631       0.99631
    iter=10       133.31      108.3       0.99555       0.99555       0.99555       0.99555       0.99555       0.99555       0.99555       0.99555       0.99555       0.99555       0.99555       0.99555       0.99555       0.99555       0.99555
    iter=11       125.29     85.513       0.95432       0.95432       0.95432       0.95432       0.95432       0.95432       0.95432       0.95432       0.95432       0.95432       0.95432       0.95432       0.95432       0.95432       0.95432
    iter=12       118.63     74.263       0.99261       0.99261       0.99261       0.99261       0.99261       0.99261       0.99261       0.99261       0.99261       0.99261       0.99261       0.99261       0.99261       0.99261       0.99261
    iter=13       112.98     61.697       0.93366       0.93366       0.93366       0.93366       0.93366       0.93366       0.93366       0.93366       0.93366       0.93366       0.93366       0.93366       0.93366       0.93366       0.93366
    iter=14        108.1     55.945       0.99009       0.99009       0.99009       0.99009       0.99009       0.99009       0.99009       0.99009       0.99009       0.99009       0.99009       0.99009       0.99009       0.99009       0.99009
    iter=15       103.77     49.547       0.95583       0.95583       0.95583       0.95583       0.95583       0.95583       0.95583       0.95583       0.95583       0.95583       0.95583       0.95583       0.95583       0.95583       0.95583
    iter=16       99.858     46.646       0.94995       0.94995       0.94995       0.94995       0.94995       0.94995       0.94995       0.94995       0.94995       0.94995       0.94995       0.94995       0.94995       0.94995       0.94995
    iter=17        96.25     42.088       0.90427       0.90427       0.90427       0.90427       0.90427       0.90427       0.90427       0.90427       0.90427       0.90427       0.90427       0.90427       0.90427       0.90427       0.90427
    iter=18        92.87     43.999         0.978         0.978         0.978         0.978         0.978         0.978         0.978         0.978         0.978         0.978         0.978         0.978         0.978         0.978         0.978
    iter=19       89.661     52.898       0.93433       0.93433       0.93433       0.93433       0.93433       0.93433       0.93433       0.93433       0.93433       0.93433       0.93433       0.93433       0.93433       0.93433       0.93433
    iter=20       86.585     63.056       0.93567       0.93567       0.93567       0.93567       0.93567       0.93567       0.93567       0.93567       0.93567       0.93567       0.93567       0.93567       0.93567       0.93567       0.93567
    iter=21       83.612     43.727       0.90091       0.90091       0.90091       0.90091       0.90091       0.90091       0.90091       0.90091       0.90091       0.90091       0.90091       0.90091       0.90091       0.90091       0.90091
    iter=22       80.723     57.013       0.86765       0.86765       0.86765       0.86765       0.86765       0.86765       0.86765       0.86765       0.86765       0.86765       0.86765       0.86765       0.86765       0.86765       0.86765
    iter=23       77.908     126.51        0.9041        0.9041        0.9041        0.9041        0.9041        0.9041        0.9041        0.9041        0.9041        0.9041        0.9041        0.9041        0.9041        0.9041        0.9041
    iter=24       75.157     38.843       0.78032       0.78032       0.78032       0.78032       0.78032       0.78032       0.78032       0.78032       0.78032       0.78032       0.78032       0.78032       0.78032       0.78032       0.78032
    iter=25       72.466      59.07       0.83994       0.83994       0.83994       0.83994       0.83994       0.83994       0.83994       0.83994       0.83994       0.83994       0.83994       0.83994       0.83994       0.83994       0.83994
    iter=26       69.833     150.14       0.91048       0.91048       0.91048       0.91048       0.91048       0.91048       0.91048       0.91048       0.91048       0.91048       0.91048       0.91048       0.91048       0.91048       0.91048
    iter=27       67.257     48.473       0.78116       0.78116       0.78116       0.78116       0.78116       0.78116       0.78116       0.78116       0.78116       0.78116       0.78116       0.78116       0.78116       0.78116       0.78116
    iter=28       64.737     56.406       0.74756       0.74756       0.74756       0.74756       0.74756       0.74756       0.74756       0.74756       0.74756       0.74756       0.74756       0.74756       0.74756       0.74756       0.74756
    iter=29       62.275      50.54       0.73396       0.73396       0.73396       0.73396       0.73396       0.73396       0.73396       0.73396       0.73396       0.73396       0.73396       0.73396       0.73396       0.73396       0.73396
    iter=30       59.872     79.692       0.75579       0.75579       0.75579       0.75579       0.75579       0.75579       0.75579       0.75579       0.75579       0.75579       0.75579       0.75579       0.75579       0.75579       0.75579
    iter=31       57.527        140       0.78435       0.78435       0.78435       0.78435       0.78435       0.78435       0.78435       0.78435       0.78435       0.78435       0.78435       0.78435       0.78435       0.78435       0.78435
    iter=32       55.242     53.339       0.68408       0.68408       0.68408       0.68408       0.68408       0.68408       0.68408       0.68408       0.68408       0.68408       0.68408       0.68408       0.68408       0.68408       0.68408
    iter=33       53.019      32.92       0.66443       0.66443       0.66443       0.66443       0.66443       0.66443       0.66443       0.66443       0.66443       0.66443       0.66443       0.66443       0.66443       0.66443       0.66443
    iter=34       50.857     108.42       0.69399       0.69399       0.69399       0.69399       0.69399       0.69399       0.69399       0.69399       0.69399       0.69399       0.69399       0.69399       0.69399       0.69399       0.69399
    iter=35       48.759     132.85       0.76218       0.76218       0.76218       0.76218       0.76218       0.76218       0.76218       0.76218       0.76218       0.76218       0.76218       0.76218       0.76218       0.76218       0.76218
    iter=36       46.724     36.334        0.6876        0.6876        0.6876        0.6876        0.6876        0.6876        0.6876        0.6876        0.6876        0.6876        0.6876        0.6876        0.6876        0.6876        0.6876
    iter=37        44.75     39.503          0.63          0.63          0.63          0.63          0.63          0.63          0.63          0.63          0.63          0.63          0.63          0.63          0.63          0.63          0.63
    iter=38        42.84     112.74       0.70037       0.70037       0.70037       0.70037       0.70037       0.70037       0.70037       0.70037       0.70037       0.70037       0.70037       0.70037       0.70037       0.70037       0.70037
    iter=39       40.994     61.555       0.63823       0.63823       0.63823       0.63823       0.63823       0.63823       0.63823       0.63823       0.63823       0.63823       0.63823       0.63823       0.63823       0.63823       0.63823
    iter=40       39.209     28.349       0.59237       0.59237       0.59237       0.59237       0.59237       0.59237       0.59237       0.59237       0.59237       0.59237       0.59237       0.59237       0.59237       0.59237       0.59237
    iter=41       37.485     39.845       0.58263       0.58263       0.58263       0.58263       0.58263       0.58263       0.58263       0.58263       0.58263       0.58263       0.58263       0.58263       0.58263       0.58263       0.58263
    iter=42       35.822     56.801       0.62529       0.62529       0.62529       0.62529       0.62529       0.62529       0.62529       0.62529       0.62529       0.62529       0.62529       0.62529       0.62529       0.62529       0.62529
    iter=43        34.22     32.535       0.52385       0.52385       0.52385       0.52385       0.52385       0.52385       0.52385       0.52385       0.52385       0.52385       0.52385       0.52385       0.52385       0.52385       0.52385
    iter=44       32.677     29.328       0.50218       0.50218       0.50218       0.50218       0.50218       0.50218       0.50218       0.50218       0.50218       0.50218       0.50218       0.50218       0.50218       0.50218       0.50218
    iter=45       31.192     165.52       0.61891       0.61891       0.61891       0.61891       0.61891       0.61891       0.61891       0.61891       0.61891       0.61891       0.61891       0.61891       0.61891       0.61891       0.61891
    iter=46       29.762     64.827       0.53695       0.53695       0.53695       0.53695       0.53695       0.53695       0.53695       0.53695       0.53695       0.53695       0.53695       0.53695       0.53695       0.53695       0.53695
    iter=47       28.389     30.908       0.51058       0.51058       0.51058       0.51058       0.51058       0.51058       0.51058       0.51058       0.51058       0.51058       0.51058       0.51058       0.51058       0.51058       0.51058
    iter=48       27.069     30.623       0.48354       0.48354       0.48354       0.48354       0.48354       0.48354       0.48354       0.48354       0.48354       0.48354       0.48354       0.48354       0.48354       0.48354       0.48354
    iter=49       25.802     25.537       0.45432       0.45432       0.45432       0.45432       0.45432       0.45432       0.45432       0.45432       0.45432       0.45432       0.45432       0.45432       0.45432       0.45432       0.45432
    iter=50       24.587     27.538       0.42778       0.42778       0.42778       0.42778       0.42778       0.42778       0.42778       0.42778       0.42778       0.42778       0.42778       0.42778       0.42778       0.42778       0.42778
    iter=96        2.025     7.6564      0.038126      0.038126      0.038126      0.038126      0.038126      0.038126      0.038126      0.038126      0.038126      0.038126      0.038126      0.038126      0.038126      0.038126      0.038126
    iter=97       1.9081     7.0423      0.035438      0.035438      0.035438      0.035438      0.035438      0.035438      0.035438      0.035438      0.035438      0.035438      0.035438      0.035438      0.035438      0.035438      0.035438
    iter=98       1.7975     6.0361      0.031911      0.031911      0.031911      0.031911      0.031911      0.031911      0.031911      0.031911      0.031911      0.031911      0.031911      0.031911      0.031911      0.031911      0.031911
    iter=99       1.6931     5.9927      0.027041      0.027041      0.027041      0.027041      0.027041      0.027041      0.027041      0.027041      0.027041      0.027041      0.027041      0.027041      0.027041      0.027041      0.027041
    iter=100      1.5944     5.6344       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133
    iter=101      1.5013     5.2122      0.020658      0.020658      0.020658      0.020658      0.020658      0.020658      0.020658      0.020658      0.020658      0.020658      0.020658      0.020658      0.020658      0.020658      0.020658
    iter=102      1.4133     4.6803      0.016124      0.016124      0.016124      0.016124      0.016124      0.016124      0.016124      0.016124      0.016124      0.016124      0.016124      0.016124      0.016124      0.016124      0.016124
    iter=103      1.3303     14.989       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133       0.02133
    iter=104       1.252     5.3398      0.013268      0.013268      0.013268      0.013268      0.013268      0.013268      0.013268      0.013268      0.013268      0.013268      0.013268      0.013268      0.013268      0.013268      0.013268
    iter=105      1.1782     3.1465     0.0087336     0.0087336     0.0087336     0.0087336     0.0087336     0.0087336     0.0087336     0.0087336     0.0087336     0.0087336     0.0087336     0.0087336     0.0087336     0.0087336     0.0087336
    iter=106      1.1085     1.7171      0.003695      0.003695      0.003695      0.003695      0.003695      0.003695      0.003695      0.003695      0.003695      0.003695      0.003695      0.003695      0.003695      0.003695      0.003695
    iter=107      1.0429     1.9799     0.0040309     0.0040309     0.0040309     0.0040309     0.0040309     0.0040309     0.0040309     0.0040309     0.0040309     0.0040309     0.0040309     0.0040309     0.0040309     0.0040309     0.0040309
    iter=108     0.98099     2.6591     0.0033591     0.0033591     0.0033591     0.0033591     0.0033591     0.0033591     0.0033591     0.0033591     0.0033591     0.0033591     0.0033591     0.0033591     0.0033591     0.0033591     0.0033591
    iter=109      0.9227     1.4126     0.0025193     0.0025193     0.0025193     0.0025193     0.0025193     0.0025193     0.0025193     0.0025193     0.0025193     0.0025193     0.0025193     0.0025193     0.0025193     0.0025193     0.0025193
    iter=110      0.8678     1.1829     0.0016795     0.0016795     0.0016795     0.0016795     0.0016795     0.0016795     0.0016795     0.0016795     0.0016795     0.0016795     0.0016795     0.0016795     0.0016795     0.0016795     0.0016795
    iter=111     0.81609     1.6275     0.0023514     0.0023514     0.0023514     0.0023514     0.0023514     0.0023514     0.0023514     0.0023514     0.0023514     0.0023514     0.0023514     0.0023514     0.0023514     0.0023514     0.0023514
    iter=112     0.76741     4.6091     0.0018475     0.0018475     0.0018475     0.0018475     0.0018475     0.0018475     0.0018475     0.0018475     0.0018475     0.0018475     0.0018475     0.0018475     0.0018475     0.0018475     0.0018475
    iter=113      0.7216          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=114     0.67848    0.90436    0.00083977    0.00083977    0.00083977    0.00083977    0.00083977    0.00083977    0.00083977    0.00083977    0.00083977    0.00083977    0.00083977    0.00083977    0.00083977    0.00083977    0.00083977
    iter=115      0.6379     0.8397     0.0010077     0.0010077     0.0010077     0.0010077     0.0010077     0.0010077     0.0010077     0.0010077     0.0010077     0.0010077     0.0010077     0.0010077     0.0010077     0.0010077     0.0010077
    iter=116     0.59973      1.188     0.0015116     0.0015116     0.0015116     0.0015116     0.0015116     0.0015116     0.0015116     0.0015116     0.0015116     0.0015116     0.0015116     0.0015116     0.0015116     0.0015116     0.0015116
    iter=117     0.56383    0.49919    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591
    iter=118     0.53006     4.1181    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591    0.00033591
    iter=119      0.4983     0.3778    0.00016795    0.00016795    0.00016795    0.00016795    0.00016795    0.00016795    0.00016795    0.00016795    0.00016795    0.00016795    0.00016795    0.00016795    0.00016795    0.00016795    0.00016795
    iter=120     0.46843          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=121     0.44035          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=122     0.41395          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=123     0.38912          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=124     0.36578          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=125     0.34384          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=126     0.32322          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=127     0.30383          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=128      0.2856          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=129     0.26846          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=130     0.25236          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=131     0.23722          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=132     0.22298          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=133     0.20961          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=134     0.19703          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=135     0.18521          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=136      0.1741          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=137     0.16365          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=138     0.15383          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=139      0.1446          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=140     0.13592          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=141     0.12777          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=142      0.1201          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=143      0.1129          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=144     0.10612          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=145    0.099756          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=0.2,b=0              -41.2         -41.2         -41.2         -41.2        -41.2         -41.2         -41.2         -41.2         -41.2         -41.2         -41.2         -41.2  
    coh2:k=0.210003,b=0       -41.092       -41.092       -41.092       -41.092      -41.092       -41.092       -41.092       -41.092       -41.092       -41.092       -41.092       -41.092  
    coh3:k=0.220006,b=0       -40.988       -40.988       -40.988       -40.988      -40.988       -40.988       -40.988       -40.988       -40.988       -40.988       -40.988       -40.988  
    coh4:k=0.230009,b=0       -40.887       -40.887       -40.887       -40.887      -40.887       -40.887       -40.887       -40.887       -40.887       -40.887       -40.887       -40.887  
    coh5:k=0.240012,b=0       -40.787       -40.787       -40.787       -40.787      -40.787       -40.787       -40.787       -40.787       -40.787       -40.787       -40.787       -40.787  
    coh6:k=0.250015,b=0        -40.69        -40.69        -40.69        -40.69       -40.69        -40.69        -40.69        -40.69        -40.69        -40.69        -40.69        -40.69  
    coh7:k=0.260018,b=0       -40.594       -40.594       -40.594       -40.594      -40.594       -40.594       -40.594       -40.594       -40.594       -40.594       -40.594       -40.594  
    coh8:k=0.270021,b=0         -40.5         -40.5         -40.5         -40.5        -40.5         -40.5         -40.5         -40.5         -40.5         -40.5         -40.5         -40.5  
    coh9:k=0.280024,b=0       -40.406       -40.406       -40.406       -40.406      -40.406       -40.406       -40.406       -40.406       -40.406       -40.406       -40.406       -40.406  
    coh10:k=0.290027,b=0      -40.314       -40.314       -40.314       -40.314      -40.314       -40.314       -40.314       -40.314       -40.314       -40.314       -40.314       -40.314  
    coh11:k=0.300029,b=0      -40.223       -40.223       -40.223       -40.223      -40.223       -40.223       -40.223       -40.223       -40.223       -40.223       -40.223       -40.223  
    coh12:k=0.310032,b=0      -40.134       -40.134       -40.134       -40.134      -40.134       -40.134       -40.134       -40.134       -40.134       -40.134       -40.134       -40.134  
    coh13:k=0.320035,b=0      -40.045       -40.045       -40.045       -40.045      -40.045       -40.045       -40.045       -40.045       -40.045       -40.045       -40.045       -40.045  
    coh14:k=0.330038,b=0      -39.957       -39.957       -39.957       -39.957      -39.957       -39.957       -39.957       -39.957       -39.957       -39.957       -39.957       -39.957  
    coh15:k=0.340041,b=0       -39.87        -39.87        -39.87        -39.87       -39.87        -39.87        -39.87        -39.87        -39.87        -39.87        -39.87        -39.87  
    coh16:k=0.350044,b=0      -39.783       -39.783       -39.783       -39.783      -39.783       -39.783       -39.783       -39.783       -39.783       -39.783       -39.783       -39.783  
    coh17:k=0.360047,b=0      -39.698       -39.698       -39.698       -39.698      -39.698       -39.698       -39.698       -39.698       -39.698       -39.698       -39.698       -39.698  
    coh18:k=0.37005,b=0       -39.613       -39.613       -39.613       -39.613      -39.613       -39.613       -39.613       -39.613       -39.613       -39.613       -39.613       -39.613  
    coh19:k=0.380053,b=0      -39.529       -39.529       -39.529       -39.529      -39.529       -39.529       -39.529       -39.529       -39.529       -39.529       -39.529       -39.529  
    coh20:k=0.390056,b=0      -39.446       -39.446       -39.446       -39.446      -39.446       -39.446       -39.446       -39.446       -39.446       -39.446       -39.446       -39.446  
    coh21:k=0.400059,b=0      -39.364       -39.364       -39.364       -39.364      -39.364       -39.364       -39.364       -39.364       -39.364       -39.364       -39.364       -39.364  
    coh22:k=0.410062,b=0      -39.282       -39.282       -39.282       -39.282      -39.282       -39.282       -39.282       -39.282       -39.282       -39.282       -39.282       -39.282  
    coh23:k=0.420065,b=0      -39.201       -39.201       -39.201       -39.201      -39.201       -39.201       -39.201       -39.201       -39.201       -39.201       -39.201       -39.201  
    coh24:k=0.430068,b=0       -39.12        -39.12        -39.12        -39.12       -39.12        -39.12        -39.12        -39.12        -39.12        -39.12        -39.12        -39.12  
    coh25:k=0.440071,b=0       -39.04        -39.04        -39.04        -39.04       -39.04        -39.04        -39.04        -39.04        -39.04        -39.04        -39.04        -39.04  
    coh26:k=0.450074,b=0      -38.961       -38.961       -38.961       -38.961      -38.961       -38.961       -38.961       -38.961       -38.961       -38.961       -38.961       -38.961  
    coh27:k=0.460077,b=0      -38.882       -38.882       -38.882       -38.882      -38.882       -38.882       -38.882       -38.882       -38.882       -38.882       -38.882       -38.882  
    coh28:k=0.47008,b=0       -38.804       -38.804       -38.804       -38.804      -38.804       -38.804       -38.804       -38.804       -38.804       -38.804       -38.804       -38.804  
    coh29:k=0.480083,b=0      -38.726       -38.726       -38.726       -38.726      -38.726       -38.726       -38.726       -38.726       -38.726       -38.726       -38.726       -38.726  
    coh30:k=0.490086,b=0      -38.649       -38.649       -38.649       -38.649      -38.649       -38.649       -38.649       -38.649       -38.649       -38.649       -38.649       -38.649  
    coh31:k=0.500088,b=0      -38.573       -38.573       -38.573       -38.573      -38.573       -38.573       -38.573       -38.573       -38.573       -38.573       -38.573       -38.573  
    coh32:k=0.510091,b=0      -38.497       -38.497       -38.497       -38.497      -38.497       -38.497       -38.497       -38.497       -38.497       -38.497       -38.497       -38.497  
    coh33:k=0.520094,b=0      -38.421       -38.421       -38.421       -38.421      -38.421       -38.421       -38.421       -38.421       -38.421       -38.421       -38.421       -38.421  
    coh34:k=0.530097,b=0      -38.346       -38.346       -38.346       -38.346      -38.346       -38.346       -38.346       -38.346       -38.346       -38.346       -38.346       -38.346  
    coh35:k=0.5401,b=0        -38.271       -38.271       -38.271       -38.271      -38.271       -38.271       -38.271       -38.271       -38.271       -38.271       -38.271       -38.271  
    coh36:k=0.550103,b=0      -38.197       -38.197       -38.197       -38.197      -38.197       -38.197       -38.197       -38.197       -38.197       -38.197       -38.197       -38.197  
    coh37:k=0.560106,b=0      -38.124       -38.124       -38.124       -38.124      -38.124       -38.124       -38.124       -38.124       -38.124       -38.124       -38.124       -38.124  
    coh38:k=0.570109,b=0       -38.05        -38.05        -38.05        -38.05       -38.05        -38.05        -38.05        -38.05        -38.05        -38.05        -38.05        -38.05  
    coh39:k=0.580112,b=0      -37.977       -37.977       -37.977       -37.977      -37.977       -37.977       -37.977       -37.977       -37.977       -37.977       -37.977       -37.977  
    coh40:k=0.590115,b=0      -37.905       -37.905       -37.905       -37.905      -37.905       -37.905       -37.905       -37.905       -37.905       -37.905       -37.905       -37.905  
    coh41:k=0.600118,b=0      -37.833       -37.833       -37.833       -37.833      -37.833       -37.833       -37.833       -37.833       -37.833       -37.833       -37.833       -37.833  
    coh42:k=0.610121,b=0      -37.762       -37.762       -37.762       -37.762      -37.762       -37.762       -37.762       -37.762       -37.762       -37.762       -37.762       -37.762  
    coh43:k=0.620124,b=0       -37.69        -37.69        -37.69        -37.69       -37.69        -37.69        -37.69        -37.69        -37.69        -37.69        -37.69        -37.69  
    coh44:k=0.630127,b=0       -37.62        -37.62        -37.62        -37.62       -37.62        -37.62        -37.62        -37.62        -37.62        -37.62        -37.62        -37.62  
    coh45:k=0.64013,b=0       -37.549       -37.549       -37.549       -37.549      -37.549       -37.549       -37.549       -37.549       -37.549       -37.549       -37.549       -37.549  
    coh46:k=0.650133,b=0      -37.479       -37.479       -37.479       -37.479      -37.479       -37.479       -37.479       -37.479       -37.479       -37.479       -37.479       -37.479  
    coh47:k=0.660136,b=0       -37.41        -37.41        -37.41        -37.41       -37.41        -37.41        -37.41        -37.41        -37.41        -37.41        -37.41        -37.41  
    coh48:k=0.670139,b=0      -37.341       -37.341       -37.341       -37.341      -37.341       -37.341       -37.341       -37.341       -37.341       -37.341       -37.341       -37.341  
    coh49:k=0.680142,b=0      -37.272       -37.272       -37.272       -37.272      -37.272       -37.272       -37.272       -37.272       -37.272       -37.272       -37.272       -37.272  
    coh50:k=0.690144,b=0      -37.203       -37.203       -37.203       -37.203      -37.203       -37.203       -37.203       -37.203       -37.203       -37.203       -37.203       -37.203  
    coh5905:k=59.2574,b=0      9.7218        9.7218        9.7218        9.7218       9.7218        9.7218        9.7218        9.7218        9.7218        9.7218        9.7218        9.7218  
    coh5906:k=59.2674,b=0      9.7219        9.7219        9.7219        9.7219       9.7219        9.7219        9.7219        9.7219        9.7219        9.7219        9.7219        9.7219  
    coh5907:k=59.2774,b=0      9.7219        9.7219        9.7219        9.7219       9.7219        9.7219        9.7219        9.7219        9.7219        9.7219        9.7219        9.7219  
    coh5908:k=59.2874,b=0      9.7221        9.7221        9.7221        9.7221       9.7221        9.7221        9.7221        9.7221        9.7221        9.7221        9.7221        9.7221  
    coh5909:k=59.2974,b=0      9.7238        9.7238        9.7238        9.7238       9.7238        9.7238        9.7238        9.7238        9.7238        9.7238        9.7238        9.7238  
    coh5910:k=59.3074,b=0      9.7256        9.7256        9.7256        9.7256       9.7256        9.7256        9.7256        9.7256        9.7256        9.7256        9.7256        9.7256  
    coh5911:k=59.3174,b=0      9.7264        9.7264        9.7264        9.7264       9.7264        9.7264        9.7264        9.7264        9.7264        9.7264        9.7264        9.7264  
    coh5912:k=59.3274,b=0      9.7264        9.7264        9.7264        9.7264       9.7264        9.7264        9.7264        9.7264        9.7264        9.7264        9.7264        9.7264  
    coh5913:k=59.3374,b=0      9.7265        9.7265        9.7265        9.7265       9.7265        9.7265        9.7265        9.7265        9.7265        9.7265        9.7265        9.7265  
    coh5914:k=59.3474,b=0      9.7266        9.7266        9.7266        9.7266       9.7266        9.7266        9.7266        9.7266        9.7266        9.7266        9.7266        9.7266  
    coh5915:k=59.3574,b=0      9.7266        9.7266        9.7266        9.7266       9.7266        9.7266        9.7266        9.7266        9.7266        9.7266        9.7266        9.7266  
    coh5916:k=59.3674,b=0      9.7267        9.7267        9.7267        9.7267       9.7267        9.7267        9.7267        9.7267        9.7267        9.7267        9.7267        9.7267  
    coh5917:k=59.3774,b=0      9.7267        9.7267        9.7267        9.7267       9.7267        9.7267        9.7267        9.7267        9.7267        9.7267        9.7267        9.7267  
    coh5918:k=59.3874,b=0      9.7268        9.7268        9.7268        9.7268       9.7268        9.7268        9.7268        9.7268        9.7268        9.7268        9.7268        9.7268  
    coh5919:k=59.3974,b=0      9.7268        9.7268        9.7268        9.7268       9.7268        9.7268        9.7268        9.7268        9.7268        9.7268        9.7268        9.7268  
    coh5920:k=59.4075,b=0      9.7269        9.7269        9.7269        9.7269       9.7269        9.7269        9.7269        9.7269        9.7269        9.7269        9.7269        9.7269  
    coh5921:k=59.4175,b=0      9.7269        9.7269        9.7269        9.7269       9.7269        9.7269        9.7269        9.7269        9.7269        9.7269        9.7269        9.7269  
    coh5922:k=59.4275,b=0       9.727         9.727         9.727         9.727        9.727         9.727         9.727         9.727         9.727         9.727         9.727         9.727  
    coh5923:k=59.4375,b=0      9.7286        9.7286        9.7286        9.7286       9.7286        9.7286        9.7286        9.7286        9.7286        9.7286        9.7286        9.7286  
    coh5924:k=59.4475,b=0      9.7303        9.7303        9.7303        9.7303       9.7303        9.7303        9.7303        9.7303        9.7303        9.7303        9.7303        9.7303  
    coh5925:k=59.4575,b=0      9.7313        9.7313        9.7313        9.7313       9.7313        9.7313        9.7313        9.7313        9.7313        9.7313        9.7313        9.7313  
    coh5926:k=59.4675,b=0      9.7313        9.7313        9.7313        9.7313       9.7313        9.7313        9.7313        9.7313        9.7313        9.7313        9.7313        9.7313  
    coh5927:k=59.4775,b=0      9.7314        9.7314        9.7314        9.7314       9.7314        9.7314        9.7314        9.7314        9.7314        9.7314        9.7314        9.7314  
    coh5928:k=59.4875,b=0      9.7315        9.7315        9.7315        9.7315       9.7315        9.7315        9.7315        9.7315        9.7315        9.7315        9.7315        9.7315  
    coh5929:k=59.4975,b=0      9.7315        9.7315        9.7315        9.7315       9.7315        9.7315        9.7315        9.7315        9.7315        9.7315        9.7315        9.7315  
    coh5930:k=59.5075,b=0      9.7316        9.7316        9.7316        9.7316       9.7316        9.7316        9.7316        9.7316        9.7316        9.7316        9.7316        9.7316  
    coh5931:k=59.5175,b=0      9.7316        9.7316        9.7316        9.7316       9.7316        9.7316        9.7316        9.7316        9.7316        9.7316        9.7316        9.7316  
    coh5932:k=59.5275,b=0      9.7317        9.7317        9.7317        9.7317       9.7317        9.7317        9.7317        9.7317        9.7317        9.7317        9.7317        9.7317  
    coh5933:k=59.5375,b=0      9.7317        9.7317        9.7317        9.7317       9.7317        9.7317        9.7317        9.7317        9.7317        9.7317        9.7317        9.7317  
    coh5934:k=59.5475,b=0      9.7318        9.7318        9.7318        9.7318       9.7318        9.7318        9.7318        9.7318        9.7318        9.7318        9.7318        9.7318  
    coh5935:k=59.5575,b=0      9.7318        9.7318        9.7318        9.7318       9.7318        9.7318        9.7318        9.7318        9.7318        9.7318        9.7318        9.7318  
    coh5936:k=59.5675,b=0      9.7319        9.7319        9.7319        9.7319       9.7319        9.7319        9.7319        9.7319        9.7319        9.7319        9.7319        9.7319  
    coh5937:k=59.5775,b=0      9.7331        9.7331        9.7331        9.7331       9.7331        9.7331        9.7331        9.7331        9.7331        9.7331        9.7331        9.7331  
    coh5938:k=59.5875,b=0      9.7348        9.7348        9.7348        9.7348       9.7348        9.7348        9.7348        9.7348        9.7348        9.7348        9.7348        9.7348  
    coh5939:k=59.5975,b=0      9.7361        9.7361        9.7361        9.7361       9.7361        9.7361        9.7361        9.7361        9.7361        9.7361        9.7361        9.7361  
    coh5940:k=59.6075,b=0      9.7362        9.7362        9.7362        9.7362       9.7362        9.7362        9.7362        9.7362        9.7362        9.7362        9.7362        9.7362  
    coh5941:k=59.6175,b=0      9.7362        9.7362        9.7362        9.7362       9.7362        9.7362        9.7362        9.7362        9.7362        9.7362        9.7362        9.7362  
    coh5942:k=59.6275,b=0      9.7363        9.7363        9.7363        9.7363       9.7363        9.7363        9.7363        9.7363        9.7363        9.7363        9.7363        9.7363  
    coh5943:k=59.6375,b=0      9.7363        9.7363        9.7363        9.7363       9.7363        9.7363        9.7363        9.7363        9.7363        9.7363        9.7363        9.7363  
    coh5944:k=59.6475,b=0      9.7364        9.7364        9.7364        9.7364       9.7364        9.7364        9.7364        9.7364        9.7364        9.7364        9.7364        9.7364  
    coh5945:k=59.6575,b=0      9.7364        9.7364        9.7364        9.7364       9.7364        9.7364        9.7364        9.7364        9.7364        9.7364        9.7364        9.7364  
    coh5946:k=59.6675,b=0      9.7365        9.7365        9.7365        9.7365       9.7365        9.7365        9.7365        9.7365        9.7365        9.7365        9.7365        9.7365  
    coh5947:k=59.6775,b=0      9.7365        9.7365        9.7365        9.7365       9.7365        9.7365        9.7365        9.7365        9.7365        9.7365        9.7365        9.7365  
    coh5948:k=59.6875,b=0      9.7366        9.7366        9.7366        9.7366       9.7366        9.7366        9.7366        9.7366        9.7366        9.7366        9.7366        9.7366  
    coh5949:k=59.6975,b=0      9.7366        9.7366        9.7366        9.7366       9.7366        9.7366        9.7366        9.7366        9.7366        9.7366        9.7366        9.7366  
    coh5950:k=59.7075,b=0      9.7367        9.7367        9.7367        9.7367       9.7367        9.7367        9.7367        9.7367        9.7367        9.7367        9.7367        9.7367  
    coh5951:k=59.7175,b=0      9.7375        9.7375        9.7375        9.7375       9.7375        9.7375        9.7375        9.7375        9.7375        9.7375        9.7375        9.7375  
    coh5952:k=59.7275,b=0      9.7392        9.7392        9.7392        9.7392       9.7392        9.7392        9.7392        9.7392        9.7392        9.7392        9.7392        9.7392  
    coh5953:k=59.7375,b=0      9.7408        9.7408        9.7408        9.7408       9.7408        9.7408        9.7408        9.7408        9.7408        9.7408        9.7408        9.7408  
    coh5954:k=59.7476,b=0      9.7409        9.7409        9.7409        9.7409       9.7409        9.7409        9.7409        9.7409        9.7409        9.7409        9.7409        9.7409  

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=0.2,b=0                2e-07         2e-07         2e-07         2e-07         2e-07         2e-07         2e-07         2e-07         2e-07         2e-07         2e-07         2e-07
    coh2:k=0.210003,b=0         2.1e-07       2.1e-07       2.1e-07       2.1e-07       2.1e-07       2.1e-07       2.1e-07       2.1e-07       2.1e-07       2.1e-07       2.1e-07       2.1e-07
    coh3:k=0.220006,b=0      9.9003e-06    9.9003e-06    9.9003e-06    9.9003e-06    9.9003e-06    9.9003e-06    9.9003e-06    9.9003e-06    9.9003e-06    9.9003e-06    9.9003e-06    9.9003e-06
    coh4:k=0.230009,b=0      1.4951e-05    1.4951e-05    1.4951e-05    1.4951e-05    1.4951e-05    1.4951e-05    1.4951e-05    1.4951e-05    1.4951e-05    1.4951e-05    1.4951e-05    1.4951e-05
    coh5:k=0.240012,b=0      2.0881e-05    2.0881e-05    2.0881e-05    2.0881e-05    2.0881e-05    2.0881e-05    2.0881e-05    2.0881e-05    2.0881e-05    2.0881e-05    2.0881e-05    2.0881e-05
    coh6:k=0.250015,b=0      2.9752e-05    2.9752e-05    2.9752e-05    2.9752e-05    2.9752e-05    2.9752e-05    2.9752e-05    2.9752e-05    2.9752e-05    2.9752e-05    2.9752e-05    2.9752e-05
    coh7:k=0.260018,b=0      3.8223e-05    3.8223e-05    3.8223e-05    3.8223e-05    3.8223e-05    3.8223e-05    3.8223e-05    3.8223e-05    3.8223e-05    3.8223e-05    3.8223e-05    3.8223e-05
    coh8:k=0.270021,b=0      4.5093e-05    4.5093e-05    4.5093e-05    4.5093e-05    4.5093e-05    4.5093e-05    4.5093e-05    4.5093e-05    4.5093e-05    4.5093e-05    4.5093e-05    4.5093e-05
    coh9:k=0.280024,b=0      5.2924e-05    5.2924e-05    5.2924e-05    5.2924e-05    5.2924e-05    5.2924e-05    5.2924e-05    5.2924e-05    5.2924e-05    5.2924e-05    5.2924e-05    5.2924e-05
    coh10:k=0.290027,b=0     6.0616e-05    6.0616e-05    6.0616e-05    6.0616e-05    6.0616e-05    6.0616e-05    6.0616e-05    6.0616e-05    6.0616e-05    6.0616e-05    6.0616e-05    6.0616e-05
    coh11:k=0.300029,b=0     6.9907e-05    6.9907e-05    6.9907e-05    6.9907e-05    6.9907e-05    6.9907e-05    6.9907e-05    6.9907e-05    6.9907e-05    6.9907e-05    6.9907e-05    6.9907e-05
    coh12:k=0.310032,b=0     7.8438e-05    7.8438e-05    7.8438e-05    7.8438e-05    7.8438e-05    7.8438e-05    7.8438e-05    7.8438e-05    7.8438e-05    7.8438e-05    7.8438e-05    7.8438e-05
    coh13:k=0.320035,b=0      8.609e-05     8.609e-05     8.609e-05     8.609e-05     8.609e-05     8.609e-05     8.609e-05     8.609e-05     8.609e-05     8.609e-05     8.609e-05     8.609e-05
    coh14:k=0.330038,b=0     9.4061e-05    9.4061e-05    9.4061e-05    9.4061e-05    9.4061e-05    9.4061e-05    9.4061e-05    9.4061e-05    9.4061e-05    9.4061e-05    9.4061e-05    9.4061e-05
    coh15:k=0.340041,b=0     0.00010235    0.00010235    0.00010235    0.00010235    0.00010235    0.00010235    0.00010235    0.00010235    0.00010235    0.00010235    0.00010235    0.00010235
    coh16:k=0.350044,b=0     0.00011096    0.00011096    0.00011096    0.00011096    0.00011096    0.00011096    0.00011096    0.00011096    0.00011096    0.00011096    0.00011096    0.00011096
    coh17:k=0.360047,b=0      0.0001199     0.0001199     0.0001199     0.0001199     0.0001199     0.0001199     0.0001199     0.0001199     0.0001199     0.0001199     0.0001199     0.0001199
    coh18:k=0.37005,b=0      0.00012915    0.00012915    0.00012915    0.00012915    0.00012915    0.00012915    0.00012915    0.00012915    0.00012915    0.00012915    0.00012915    0.00012915
    coh19:k=0.380053,b=0      0.0001372     0.0001372     0.0001372     0.0001372     0.0001372     0.0001372     0.0001372     0.0001372     0.0001372     0.0001372     0.0001372     0.0001372
    coh20:k=0.390056,b=0     0.00014549    0.00014549    0.00014549    0.00014549    0.00014549    0.00014549    0.00014549    0.00014549    0.00014549    0.00014549    0.00014549    0.00014549
    coh21:k=0.400059,b=0     0.00015402    0.00015402    0.00015402    0.00015402    0.00015402    0.00015402    0.00015402    0.00015402    0.00015402    0.00015402    0.00015402    0.00015402
    coh22:k=0.410062,b=0     0.00016279    0.00016279    0.00016279    0.00016279    0.00016279    0.00016279    0.00016279    0.00016279    0.00016279    0.00016279    0.00016279    0.00016279
    coh23:k=0.420065,b=0     0.00017097    0.00017097    0.00017097    0.00017097    0.00017097    0.00017097    0.00017097    0.00017097    0.00017097    0.00017097    0.00017097    0.00017097
    coh24:k=0.430068,b=0      0.0001802     0.0001802     0.0001802     0.0001802     0.0001802     0.0001802     0.0001802     0.0001802     0.0001802     0.0001802     0.0001802     0.0001802
    coh25:k=0.440071,b=0     0.00018967    0.00018967    0.00018967    0.00018967    0.00018967    0.00018967    0.00018967    0.00018967    0.00018967    0.00018967    0.00018967    0.00018967
    coh26:k=0.450074,b=0     0.00019848    0.00019848    0.00019848    0.00019848    0.00019848    0.00019848    0.00019848    0.00019848    0.00019848    0.00019848    0.00019848    0.00019848
    coh27:k=0.460077,b=0     0.00020749    0.00020749    0.00020749    0.00020749    0.00020749    0.00020749    0.00020749    0.00020749    0.00020749    0.00020749    0.00020749    0.00020749
    coh28:k=0.47008,b=0      0.00021577    0.00021577    0.00021577    0.00021577    0.00021577    0.00021577    0.00021577    0.00021577    0.00021577    0.00021577    0.00021577    0.00021577
    coh29:k=0.480083,b=0      0.0002242     0.0002242     0.0002242     0.0002242     0.0002242     0.0002242     0.0002242     0.0002242     0.0002242     0.0002242     0.0002242     0.0002242
    coh30:k=0.490086,b=0     0.00023279    0.00023279    0.00023279    0.00023279    0.00023279    0.00023279    0.00023279    0.00023279    0.00023279    0.00023279    0.00023279    0.00023279
    coh31:k=0.500088,b=0     0.00024154    0.00024154    0.00024154    0.00024154    0.00024154    0.00024154    0.00024154    0.00024154    0.00024154    0.00024154    0.00024154    0.00024154
    coh32:k=0.510091,b=0     0.00025045    0.00025045    0.00025045    0.00025045    0.00025045    0.00025045    0.00025045    0.00025045    0.00025045    0.00025045    0.00025045    0.00025045
    coh33:k=0.520094,b=0     0.00025953    0.00025953    0.00025953    0.00025953    0.00025953    0.00025953    0.00025953    0.00025953    0.00025953    0.00025953    0.00025953    0.00025953
    coh34:k=0.530097,b=0     0.00026876    0.00026876    0.00026876    0.00026876    0.00026876    0.00026876    0.00026876    0.00026876    0.00026876    0.00026876    0.00026876    0.00026876
    coh35:k=0.5401,b=0       0.00027815    0.00027815    0.00027815    0.00027815    0.00027815    0.00027815    0.00027815    0.00027815    0.00027815    0.00027815    0.00027815    0.00027815
    coh36:k=0.550103,b=0      0.0002866     0.0002866     0.0002866     0.0002866     0.0002866     0.0002866     0.0002866     0.0002866     0.0002866     0.0002866     0.0002866     0.0002866
    coh37:k=0.560106,b=0     0.00029518    0.00029518    0.00029518    0.00029518    0.00029518    0.00029518    0.00029518    0.00029518    0.00029518    0.00029518    0.00029518    0.00029518
    coh38:k=0.570109,b=0     0.00030387    0.00030387    0.00030387    0.00030387    0.00030387    0.00030387    0.00030387    0.00030387    0.00030387    0.00030387    0.00030387    0.00030387
    coh39:k=0.580112,b=0     0.00031268    0.00031268    0.00031268    0.00031268    0.00031268    0.00031268    0.00031268    0.00031268    0.00031268    0.00031268    0.00031268    0.00031268
    coh40:k=0.590115,b=0     0.00032161    0.00032161    0.00032161    0.00032161    0.00032161    0.00032161    0.00032161    0.00032161    0.00032161    0.00032161    0.00032161    0.00032161
    coh41:k=0.600118,b=0     0.00033066    0.00033066    0.00033066    0.00033066    0.00033066    0.00033066    0.00033066    0.00033066    0.00033066    0.00033066    0.00033066    0.00033066
    coh42:k=0.610121,b=0     0.00033984    0.00033984    0.00033984    0.00033984    0.00033984    0.00033984    0.00033984    0.00033984    0.00033984    0.00033984    0.00033984    0.00033984
    coh43:k=0.620124,b=0     0.00034913    0.00034913    0.00034913    0.00034913    0.00034913    0.00034913    0.00034913    0.00034913    0.00034913    0.00034913    0.00034913    0.00034913
    coh44:k=0.630127,b=0     0.00035854    0.00035854    0.00035854    0.00035854    0.00035854    0.00035854    0.00035854    0.00035854    0.00035854    0.00035854    0.00035854    0.00035854
    coh45:k=0.64013,b=0      0.00036807    0.00036807    0.00036807    0.00036807    0.00036807    0.00036807    0.00036807    0.00036807    0.00036807    0.00036807    0.00036807    0.00036807
    coh46:k=0.650133,b=0     0.00037643    0.00037643    0.00037643    0.00037643    0.00037643    0.00037643    0.00037643    0.00037643    0.00037643    0.00037643    0.00037643    0.00037643
    coh47:k=0.660136,b=0     0.00038486    0.00038486    0.00038486    0.00038486    0.00038486    0.00038486    0.00038486    0.00038486    0.00038486    0.00038486    0.00038486    0.00038486
    coh48:k=0.670139,b=0     0.00039471    0.00039471    0.00039471    0.00039471    0.00039471    0.00039471    0.00039471    0.00039471    0.00039471    0.00039471    0.00039471    0.00039471
    coh49:k=0.680142,b=0     0.00040332    0.00040332    0.00040332    0.00040332    0.00040332    0.00040332    0.00040332    0.00040332    0.00040332    0.00040332    0.00040332    0.00040332
    coh50:k=0.690144,b=0     0.00041202    0.00041202    0.00041202    0.00041202    0.00041202    0.00041202    0.00041202    0.00041202    0.00041202    0.00041202    0.00041202    0.00041202
    coh5905:k=59.2574,b=0        35.723        35.723        35.723        35.723        35.723        35.723        35.723        35.723        35.723        35.723        35.723        35.723
    coh5906:k=59.2674,b=0        35.731        35.731        35.731        35.731        35.731        35.731        35.731        35.731        35.731        35.731        35.731        35.731
    coh5907:k=59.2774,b=0        35.739        35.739        35.739        35.739        35.739        35.739        35.739        35.739        35.739        35.739        35.739        35.739
    coh5908:k=59.2874,b=0        35.528        35.528        35.528        35.528        35.528        35.528        35.528        35.528        35.528        35.528        35.528        35.528
    coh5909:k=59.2974,b=0        35.541        35.541        35.541        35.541        35.541        35.541        35.541        35.541        35.541        35.541        35.541        35.541
    coh5910:k=59.3074,b=0        35.609        35.609        35.609        35.609        35.609        35.609        35.609        35.609        35.609        35.609        35.609        35.609
    coh5911:k=59.3174,b=0        35.655        35.655        35.655        35.655        35.655        35.655        35.655        35.655        35.655        35.655        35.655        35.655
    coh5912:k=59.3274,b=0        35.663        35.663        35.663        35.663        35.663        35.663        35.663        35.663        35.663        35.663        35.663        35.663
    coh5913:k=59.3374,b=0        35.671        35.671        35.671        35.671        35.671        35.671        35.671        35.671        35.671        35.671        35.671        35.671
    coh5914:k=59.3474,b=0         35.68         35.68         35.68         35.68         35.68         35.68         35.68         35.68         35.68         35.68         35.68         35.68
    coh5915:k=59.3574,b=0        35.688        35.688        35.688        35.688        35.688        35.688        35.688        35.688        35.688        35.688        35.688        35.688
    coh5916:k=59.3674,b=0        35.697        35.697        35.697        35.697        35.697        35.697        35.697        35.697        35.697        35.697        35.697        35.697
    coh5917:k=59.3774,b=0        35.705        35.705        35.705        35.705        35.705        35.705        35.705        35.705        35.705        35.705        35.705        35.705
    coh5918:k=59.3874,b=0        35.714        35.714        35.714        35.714        35.714        35.714        35.714        35.714        35.714        35.714        35.714        35.714
    coh5919:k=59.3974,b=0        35.722        35.722        35.722        35.722        35.722        35.722        35.722        35.722        35.722        35.722        35.722        35.722
    coh5920:k=59.4075,b=0         35.73         35.73         35.73         35.73         35.73         35.73         35.73         35.73         35.73         35.73         35.73         35.73
    coh5921:k=59.4175,b=0        35.739        35.739        35.739        35.739        35.739        35.739        35.739        35.739        35.739        35.739        35.739        35.739
    coh5922:k=59.4275,b=0        35.747        35.747        35.747        35.747        35.747        35.747        35.747        35.747        35.747        35.747        35.747        35.747
    coh5923:k=59.4375,b=0        35.541        35.541        35.541        35.541        35.541        35.541        35.541        35.541        35.541        35.541        35.541        35.541
    coh5924:k=59.4475,b=0          35.6          35.6          35.6          35.6          35.6          35.6          35.6          35.6          35.6          35.6          35.6          35.6
    coh5925:k=59.4575,b=0        35.654        35.654        35.654        35.654        35.654        35.654        35.654        35.654        35.654        35.654        35.654        35.654
    coh5926:k=59.4675,b=0        35.662        35.662        35.662        35.662        35.662        35.662        35.662        35.662        35.662        35.662        35.662        35.662
    coh5927:k=59.4775,b=0        35.671        35.671        35.671        35.671        35.671        35.671        35.671        35.671        35.671        35.671        35.671        35.671
    coh5928:k=59.4875,b=0        35.679        35.679        35.679        35.679        35.679        35.679        35.679        35.679        35.679        35.679        35.679        35.679
    coh5929:k=59.4975,b=0        35.687        35.687        35.687        35.687        35.687        35.687        35.687        35.687        35.687        35.687        35.687        35.687
    coh5930:k=59.5075,b=0        35.696        35.696        35.696        35.696        35.696        35.696        35.696        35.696        35.696        35.696        35.696        35.696
    coh5931:k=59.5175,b=0        35.704        35.704        35.704        35.704        35.704        35.704        35.704        35.704        35.704        35.704        35.704        35.704
    coh5932:k=59.5275,b=0        35.713        35.713        35.713        35.713        35.713        35.713        35.713        35.713        35.713        35.713        35.713        35.713
    coh5933:k=59.5375,b=0        35.721        35.721        35.721        35.721        35.721        35.721        35.721        35.721        35.721        35.721        35.721        35.721
    coh5934:k=59.5475,b=0        35.729        35.729        35.729        35.729        35.729        35.729        35.729        35.729        35.729        35.729        35.729        35.729
    coh5935:k=59.5575,b=0        35.738        35.738        35.738        35.738        35.738        35.738        35.738        35.738        35.738        35.738        35.738        35.738
    coh5936:k=59.5675,b=0        35.746        35.746        35.746        35.746        35.746        35.746        35.746        35.746        35.746        35.746        35.746        35.746
    coh5937:k=59.5775,b=0         35.54         35.54         35.54         35.54         35.54         35.54         35.54         35.54         35.54         35.54         35.54         35.54
    coh5938:k=59.5875,b=0        35.585        35.585        35.585        35.585        35.585        35.585        35.585        35.585        35.585        35.585        35.585        35.585
    coh5939:k=59.5975,b=0        35.652        35.652        35.652        35.652        35.652        35.652        35.652        35.652        35.652        35.652        35.652        35.652
    coh5940:k=59.6075,b=0        35.661        35.661        35.661        35.661        35.661        35.661        35.661        35.661        35.661        35.661        35.661        35.661
    coh5941:k=59.6175,b=0        35.669        35.669        35.669        35.669        35.669        35.669        35.669        35.669        35.669        35.669        35.669        35.669
    coh5942:k=59.6275,b=0        35.677        35.677        35.677        35.677        35.677        35.677        35.677        35.677        35.677        35.677        35.677        35.677
    coh5943:k=59.6375,b=0        35.686        35.686        35.686        35.686        35.686        35.686        35.686        35.686        35.686        35.686        35.686        35.686
    coh5944:k=59.6475,b=0        35.694        35.694        35.694        35.694        35.694        35.694        35.694        35.694        35.694        35.694        35.694        35.694
    coh5945:k=59.6575,b=0        35.703        35.703        35.703        35.703        35.703        35.703        35.703        35.703        35.703        35.703        35.703        35.703
    coh5946:k=59.6675,b=0        35.711        35.711        35.711        35.711        35.711        35.711        35.711        35.711        35.711        35.711        35.711        35.711
    coh5947:k=59.6775,b=0        35.719        35.719        35.719        35.719        35.719        35.719        35.719        35.719        35.719        35.719        35.719        35.719
    coh5948:k=59.6875,b=0        35.728        35.728        35.728        35.728        35.728        35.728        35.728        35.728        35.728        35.728        35.728        35.728
    coh5949:k=59.6975,b=0        35.736        35.736        35.736        35.736        35.736        35.736        35.736        35.736        35.736        35.736        35.736        35.736
    coh5950:k=59.7075,b=0        35.745        35.745        35.745        35.745        35.745        35.745        35.745        35.745        35.745        35.745        35.745        35.745
    coh5951:k=59.7175,b=0        35.538        35.538        35.538        35.538        35.538        35.538        35.538        35.538        35.538        35.538        35.538        35.538
    coh5952:k=59.7275,b=0        35.564        35.564        35.564        35.564        35.564        35.564        35.564        35.564        35.564        35.564        35.564        35.564
    coh5953:k=59.7375,b=0         35.65         35.65         35.65         35.65         35.65         35.65         35.65         35.65         35.65         35.65         35.65         35.65
    coh5954:k=59.7476,b=0        35.659        35.659        35.659        35.659        35.659        35.659        35.659        35.659        35.659        35.659        35.659        35.659

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=0.2,b=0            0.0001998     0.0001998     0.0001998     0.0001998     0.0001998     0.0001998     0.0001998     0.0001998     0.0001998     0.0001998     0.0001998     0.0001998
    coh2:k=0.210003,b=0      0.00020979    0.00020979    0.00020979    0.00020979    0.00020979    0.00020979    0.00020979    0.00020979    0.00020979    0.00020979    0.00020979    0.00020979
    coh3:k=0.220006,b=0       0.0098904     0.0098904     0.0098904     0.0098904     0.0098904     0.0098904     0.0098904     0.0098904     0.0098904     0.0098904     0.0098904     0.0098904
    coh4:k=0.230009,b=0        0.014936      0.014936      0.014936      0.014936      0.014936      0.014936      0.014936      0.014936      0.014936      0.014936      0.014936      0.014936
    coh5:k=0.240012,b=0         0.02086       0.02086       0.02086       0.02086       0.02086       0.02086       0.02086       0.02086       0.02086       0.02086       0.02086       0.02086
    coh6:k=0.250015,b=0        0.029722      0.029722      0.029722      0.029722      0.029722      0.029722      0.029722      0.029722      0.029722      0.029722      0.029722      0.029722
    coh7:k=0.260018,b=0        0.038184      0.038184      0.038184      0.038184      0.038184      0.038184      0.038184      0.038184      0.038184      0.038184      0.038184      0.038184
    coh8:k=0.270021,b=0        0.045048      0.045048      0.045048      0.045048      0.045048      0.045048      0.045048      0.045048      0.045048      0.045048      0.045048      0.045048
    coh9:k=0.280024,b=0        0.052872      0.052872      0.052872      0.052872      0.052872      0.052872      0.052872      0.052872      0.052872      0.052872      0.052872      0.052872
    coh10:k=0.290027,b=0       0.060555      0.060555      0.060555      0.060555      0.060555      0.060555      0.060555      0.060555      0.060555      0.060555      0.060555      0.060555
    coh11:k=0.300029,b=0       0.069837      0.069837      0.069837      0.069837      0.069837      0.069837      0.069837      0.069837      0.069837      0.069837      0.069837      0.069837
    coh12:k=0.310032,b=0        0.07836       0.07836       0.07836       0.07836       0.07836       0.07836       0.07836       0.07836       0.07836       0.07836       0.07836       0.07836
    coh13:k=0.320035,b=0       0.086003      0.086003      0.086003      0.086003      0.086003      0.086003      0.086003      0.086003      0.086003      0.086003      0.086003      0.086003
    coh14:k=0.330038,b=0       0.093967      0.093967      0.093967      0.093967      0.093967      0.093967      0.093967      0.093967      0.093967      0.093967      0.093967      0.093967
    coh15:k=0.340041,b=0        0.10225       0.10225       0.10225       0.10225       0.10225       0.10225       0.10225       0.10225       0.10225       0.10225       0.10225       0.10225
    coh16:k=0.350044,b=0        0.11085       0.11085       0.11085       0.11085       0.11085       0.11085       0.11085       0.11085       0.11085       0.11085       0.11085       0.11085
    coh17:k=0.360047,b=0        0.11978       0.11978       0.11978       0.11978       0.11978       0.11978       0.11978       0.11978       0.11978       0.11978       0.11978       0.11978
    coh18:k=0.37005,b=0         0.12902       0.12902       0.12902       0.12902       0.12902       0.12902       0.12902       0.12902       0.12902       0.12902       0.12902       0.12902
    coh19:k=0.380053,b=0        0.13706       0.13706       0.13706       0.13706       0.13706       0.13706       0.13706       0.13706       0.13706       0.13706       0.13706       0.13706
    coh20:k=0.390056,b=0        0.14535       0.14535       0.14535       0.14535       0.14535       0.14535       0.14535       0.14535       0.14535       0.14535       0.14535       0.14535
    coh21:k=0.400059,b=0        0.15387       0.15387       0.15387       0.15387       0.15387       0.15387       0.15387       0.15387       0.15387       0.15387       0.15387       0.15387
    coh22:k=0.410062,b=0        0.16263       0.16263       0.16263       0.16263       0.16263       0.16263       0.16263       0.16263       0.16263       0.16263       0.16263       0.16263
    coh23:k=0.420065,b=0         0.1708        0.1708        0.1708        0.1708        0.1708        0.1708        0.1708        0.1708        0.1708        0.1708        0.1708        0.1708
    coh24:k=0.430068,b=0        0.18002       0.18002       0.18002       0.18002       0.18002       0.18002       0.18002       0.18002       0.18002       0.18002       0.18002       0.18002
    coh25:k=0.440071,b=0        0.18948       0.18948       0.18948       0.18948       0.18948       0.18948       0.18948       0.18948       0.18948       0.18948       0.18948       0.18948
    coh26:k=0.450074,b=0        0.19828       0.19828       0.19828       0.19828       0.19828       0.19828       0.19828       0.19828       0.19828       0.19828       0.19828       0.19828
    coh27:k=0.460077,b=0        0.20729       0.20729       0.20729       0.20729       0.20729       0.20729       0.20729       0.20729       0.20729       0.20729       0.20729       0.20729
    coh28:k=0.47008,b=0         0.21555       0.21555       0.21555       0.21555       0.21555       0.21555       0.21555       0.21555       0.21555       0.21555       0.21555       0.21555
    coh29:k=0.480083,b=0        0.22397       0.22397       0.22397       0.22397       0.22397       0.22397       0.22397       0.22397       0.22397       0.22397       0.22397       0.22397
    coh30:k=0.490086,b=0        0.23256       0.23256       0.23256       0.23256       0.23256       0.23256       0.23256       0.23256       0.23256       0.23256       0.23256       0.23256
    coh31:k=0.500088,b=0         0.2413        0.2413        0.2413        0.2413        0.2413        0.2413        0.2413        0.2413        0.2413        0.2413        0.2413        0.2413
    coh32:k=0.510091,b=0         0.2502        0.2502        0.2502        0.2502        0.2502        0.2502        0.2502        0.2502        0.2502        0.2502        0.2502        0.2502
    coh33:k=0.520094,b=0        0.25927       0.25927       0.25927       0.25927       0.25927       0.25927       0.25927       0.25927       0.25927       0.25927       0.25927       0.25927
    coh34:k=0.530097,b=0        0.26849       0.26849       0.26849       0.26849       0.26849       0.26849       0.26849       0.26849       0.26849       0.26849       0.26849       0.26849
    coh35:k=0.5401,b=0          0.27787       0.27787       0.27787       0.27787       0.27787       0.27787       0.27787       0.27787       0.27787       0.27787       0.27787       0.27787
    coh36:k=0.550103,b=0        0.28632       0.28632       0.28632       0.28632       0.28632       0.28632       0.28632       0.28632       0.28632       0.28632       0.28632       0.28632
    coh37:k=0.560106,b=0        0.29488       0.29488       0.29488       0.29488       0.29488       0.29488       0.29488       0.29488       0.29488       0.29488       0.29488       0.29488
    coh38:k=0.570109,b=0        0.30356       0.30356       0.30356       0.30356       0.30356       0.30356       0.30356       0.30356       0.30356       0.30356       0.30356       0.30356
    coh39:k=0.580112,b=0        0.31237       0.31237       0.31237       0.31237       0.31237       0.31237       0.31237       0.31237       0.31237       0.31237       0.31237       0.31237
    coh40:k=0.590115,b=0        0.32129       0.32129       0.32129       0.32129       0.32129       0.32129       0.32129       0.32129       0.32129       0.32129       0.32129       0.32129
    coh41:k=0.600118,b=0        0.33033       0.33033       0.33033       0.33033       0.33033       0.33033       0.33033       0.33033       0.33033       0.33033       0.33033       0.33033
    coh42:k=0.610121,b=0         0.3395        0.3395        0.3395        0.3395        0.3395        0.3395        0.3395        0.3395        0.3395        0.3395        0.3395        0.3395
    coh43:k=0.620124,b=0        0.34878       0.34878       0.34878       0.34878       0.34878       0.34878       0.34878       0.34878       0.34878       0.34878       0.34878       0.34878
    coh44:k=0.630127,b=0        0.35818       0.35818       0.35818       0.35818       0.35818       0.35818       0.35818       0.35818       0.35818       0.35818       0.35818       0.35818
    coh45:k=0.64013,b=0         0.36771       0.36771       0.36771       0.36771       0.36771       0.36771       0.36771       0.36771       0.36771       0.36771       0.36771       0.36771
    coh46:k=0.650133,b=0        0.37605       0.37605       0.37605       0.37605       0.37605       0.37605       0.37605       0.37605       0.37605       0.37605       0.37605       0.37605
    coh47:k=0.660136,b=0        0.38447       0.38447       0.38447       0.38447       0.38447       0.38447       0.38447       0.38447       0.38447       0.38447       0.38447       0.38447
    coh48:k=0.670139,b=0        0.39432       0.39432       0.39432       0.39432       0.39432       0.39432       0.39432       0.39432       0.39432       0.39432       0.39432       0.39432
    coh49:k=0.680142,b=0        0.40292       0.40292       0.40292       0.40292       0.40292       0.40292       0.40292       0.40292       0.40292       0.40292       0.40292       0.40292
    coh50:k=0.690144,b=0         0.4116        0.4116        0.4116        0.4116        0.4116        0.4116        0.4116        0.4116        0.4116        0.4116        0.4116        0.4116
    coh5905:k=59.2574,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5906:k=59.2674,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5907:k=59.2774,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5908:k=59.2874,b=0        14.451        14.451        14.451        14.451        14.451        14.451        14.451        14.451        14.451        14.451        14.451        14.451
    coh5909:k=59.2974,b=0        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446
    coh5910:k=59.3074,b=0        14.387        14.387        14.387        14.387        14.387        14.387        14.387        14.387        14.387        14.387        14.387        14.387
    coh5911:k=59.3174,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5912:k=59.3274,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5913:k=59.3374,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5914:k=59.3474,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5915:k=59.3574,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5916:k=59.3674,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5917:k=59.3774,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5918:k=59.3874,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5919:k=59.3974,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5920:k=59.4075,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5921:k=59.4175,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5922:k=59.4275,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5923:k=59.4375,b=0        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446
    coh5924:k=59.4475,b=0        14.395        14.395        14.395        14.395        14.395        14.395        14.395        14.395        14.395        14.395        14.395        14.395
    coh5925:k=59.4575,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5926:k=59.4675,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5927:k=59.4775,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5928:k=59.4875,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5929:k=59.4975,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5930:k=59.5075,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5931:k=59.5175,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5932:k=59.5275,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5933:k=59.5375,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5934:k=59.5475,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5935:k=59.5575,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5936:k=59.5675,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5937:k=59.5775,b=0        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446        14.446
    coh5938:k=59.5875,b=0        14.409        14.409        14.409        14.409        14.409        14.409        14.409        14.409        14.409        14.409        14.409        14.409
    coh5939:k=59.5975,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5940:k=59.6075,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5941:k=59.6175,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5942:k=59.6275,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5943:k=59.6375,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5944:k=59.6475,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5945:k=59.6575,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5946:k=59.6675,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5947:k=59.6775,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5948:k=59.6875,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5949:k=59.6975,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5950:k=59.7075,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5951:k=59.7175,b=0        14.445        14.445        14.445        14.445        14.445        14.445        14.445        14.445        14.445        14.445        14.445        14.445
    coh5952:k=59.7275,b=0        14.428        14.428        14.428        14.428        14.428        14.428        14.428        14.428        14.428        14.428        14.428        14.428
    coh5953:k=59.7375,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35
    coh5954:k=59.7476,b=0         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35         14.35

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=0.2,b=0              0.0002        0.0002        0.0002        0.0002       0.0002        0.0002        0.0002        0.0002        0.0002        0.0002        0.0002        0.0002 
    coh2:k=0.210003,b=0        0.00021       0.00021       0.00021       0.00021      0.00021       0.00021       0.00021       0.00021       0.00021       0.00021       0.00021       0.00021 
    coh3:k=0.220006,b=0      0.0099003     0.0099003     0.0099003     0.0099003    0.0099003     0.0099003     0.0099003     0.0099003     0.0099003     0.0099003     0.0099003     0.0099003 
    coh4:k=0.230009,b=0       0.014951      0.014951      0.014951      0.014951     0.014951      0.014951      0.014951      0.014951      0.014951      0.014951      0.014951      0.014951 
    coh5:k=0.240012,b=0       0.020881      0.020881      0.020881      0.020881     0.020881      0.020881      0.020881      0.020881      0.020881      0.020881      0.020881      0.020881 
    coh6:k=0.250015,b=0       0.029752      0.029752      0.029752      0.029752     0.029752      0.029752      0.029752      0.029752      0.029752      0.029752      0.029752      0.029752 
    coh7:k=0.260018,b=0       0.038223      0.038223      0.038223      0.038223     0.038223      0.038223      0.038223      0.038223      0.038223      0.038223      0.038223      0.038223 
    coh8:k=0.270021,b=0       0.045093      0.045093      0.045093      0.045093     0.045093      0.045093      0.045093      0.045093      0.045093      0.045093      0.045093      0.045093 
    coh9:k=0.280024,b=0       0.052924      0.052924      0.052924      0.052924     0.052924      0.052924      0.052924      0.052924      0.052924      0.052924      0.052924      0.052924 
    coh10:k=0.290027,b=0      0.060616      0.060616      0.060616      0.060616     0.060616      0.060616      0.060616      0.060616      0.060616      0.060616      0.060616      0.060616 
    coh11:k=0.300029,b=0      0.069907      0.069907      0.069907      0.069907     0.069907      0.069907      0.069907      0.069907      0.069907      0.069907      0.069907      0.069907 
    coh12:k=0.310032,b=0      0.078438      0.078438      0.078438      0.078438     0.078438      0.078438      0.078438      0.078438      0.078438      0.078438      0.078438      0.078438 
    coh13:k=0.320035,b=0       0.08609       0.08609       0.08609       0.08609      0.08609       0.08609       0.08609       0.08609       0.08609       0.08609       0.08609       0.08609 
    coh14:k=0.330038,b=0      0.094061      0.094061      0.094061      0.094061     0.094061      0.094061      0.094061      0.094061      0.094061      0.094061      0.094061      0.094061 
    coh15:k=0.340041,b=0       0.10235       0.10235       0.10235       0.10235      0.10235       0.10235       0.10235       0.10235       0.10235       0.10235       0.10235       0.10235 
    coh16:k=0.350044,b=0       0.11096       0.11096       0.11096       0.11096      0.11096       0.11096       0.11096       0.11096       0.11096       0.11096       0.11096       0.11096 
    coh17:k=0.360047,b=0        0.1199        0.1199        0.1199        0.1199       0.1199        0.1199        0.1199        0.1199        0.1199        0.1199        0.1199        0.1199 
    coh18:k=0.37005,b=0        0.12915       0.12915       0.12915       0.12915      0.12915       0.12915       0.12915       0.12915       0.12915       0.12915       0.12915       0.12915 
    coh19:k=0.380053,b=0        0.1372        0.1372        0.1372        0.1372       0.1372        0.1372        0.1372        0.1372        0.1372        0.1372        0.1372        0.1372 
    coh20:k=0.390056,b=0       0.14549       0.14549       0.14549       0.14549      0.14549       0.14549       0.14549       0.14549       0.14549       0.14549       0.14549       0.14549 
    coh21:k=0.400059,b=0       0.15402       0.15402       0.15402       0.15402      0.15402       0.15402       0.15402       0.15402       0.15402       0.15402       0.15402       0.15402 
    coh22:k=0.410062,b=0       0.16279       0.16279       0.16279       0.16279      0.16279       0.16279       0.16279       0.16279       0.16279       0.16279       0.16279       0.16279 
    coh23:k=0.420065,b=0       0.17097       0.17097       0.17097       0.17097      0.17097       0.17097       0.17097       0.17097       0.17097       0.17097       0.17097       0.17097 
    coh24:k=0.430068,b=0        0.1802        0.1802        0.1802        0.1802       0.1802        0.1802        0.1802        0.1802        0.1802        0.1802        0.1802        0.1802 
    coh25:k=0.440071,b=0       0.18967       0.18967       0.18967       0.18967      0.18967       0.18967       0.18967       0.18967       0.18967       0.18967       0.18967       0.18967 
    coh26:k=0.450074,b=0       0.19848       0.19848       0.19848       0.19848      0.19848       0.19848       0.19848       0.19848       0.19848       0.19848       0.19848       0.19848 
    coh27:k=0.460077,b=0       0.20749       0.20749       0.20749       0.20749      0.20749       0.20749       0.20749       0.20749       0.20749       0.20749       0.20749       0.20749 
    coh28:k=0.47008,b=0        0.21577       0.21577       0.21577       0.21577      0.21577       0.21577       0.21577       0.21577       0.21577       0.21577       0.21577       0.21577 
    coh29:k=0.480083,b=0        0.2242        0.2242        0.2242        0.2242       0.2242        0.2242        0.2242        0.2242        0.2242        0.2242        0.2242        0.2242 
    coh30:k=0.490086,b=0       0.23279       0.23279       0.23279       0.23279      0.23279       0.23279       0.23279       0.23279       0.23279       0.23279       0.23279       0.23279 
    coh31:k=0.500088,b=0       0.24154       0.24154       0.24154       0.24154      0.24154       0.24154       0.24154       0.24154       0.24154       0.24154       0.24154       0.24154 
    coh32:k=0.510091,b=0       0.25045       0.25045       0.25045       0.25045      0.25045       0.25045       0.25045       0.25045       0.25045       0.25045       0.25045       0.25045 
    coh33:k=0.520094,b=0       0.25953       0.25953       0.25953       0.25953      0.25953       0.25953       0.25953       0.25953       0.25953       0.25953       0.25953       0.25953 
    coh34:k=0.530097,b=0       0.26876       0.26876       0.26876       0.26876      0.26876       0.26876       0.26876       0.26876       0.26876       0.26876       0.26876       0.26876 
    coh35:k=0.5401,b=0         0.27815       0.27815       0.27815       0.27815      0.27815       0.27815       0.27815       0.27815       0.27815       0.27815       0.27815       0.27815 
    coh36:k=0.550103,b=0        0.2866        0.2866        0.2866        0.2866       0.2866        0.2866        0.2866        0.2866        0.2866        0.2866        0.2866        0.2866 
    coh37:k=0.560106,b=0       0.29518       0.29518       0.29518       0.29518      0.29518       0.29518       0.29518       0.29518       0.29518       0.29518       0.29518       0.29518 
    coh38:k=0.570109,b=0       0.30387       0.30387       0.30387       0.30387      0.30387       0.30387       0.30387       0.30387       0.30387       0.30387       0.30387       0.30387 
    coh39:k=0.580112,b=0       0.31268       0.31268       0.31268       0.31268      0.31268       0.31268       0.31268       0.31268       0.31268       0.31268       0.31268       0.31268 
    coh40:k=0.590115,b=0       0.32161       0.32161       0.32161       0.32161      0.32161       0.32161       0.32161       0.32161       0.32161       0.32161       0.32161       0.32161 
    coh41:k=0.600118,b=0       0.33066       0.33066       0.33066       0.33066      0.33066       0.33066       0.33066       0.33066       0.33066       0.33066       0.33066       0.33066 
    coh42:k=0.610121,b=0       0.33984       0.33984       0.33984       0.33984      0.33984       0.33984       0.33984       0.33984       0.33984       0.33984       0.33984       0.33984 
    coh43:k=0.620124,b=0       0.34913       0.34913       0.34913       0.34913      0.34913       0.34913       0.34913       0.34913       0.34913       0.34913       0.34913       0.34913 
    coh44:k=0.630127,b=0       0.35854       0.35854       0.35854       0.35854      0.35854       0.35854       0.35854       0.35854       0.35854       0.35854       0.35854       0.35854 
    coh45:k=0.64013,b=0        0.36807       0.36807       0.36807       0.36807      0.36807       0.36807       0.36807       0.36807       0.36807       0.36807       0.36807       0.36807 
    coh46:k=0.650133,b=0       0.37643       0.37643       0.37643       0.37643      0.37643       0.37643       0.37643       0.37643       0.37643       0.37643       0.37643       0.37643 
    coh47:k=0.660136,b=0       0.38486       0.38486       0.38486       0.38486      0.38486       0.38486       0.38486       0.38486       0.38486       0.38486       0.38486       0.38486 
    coh48:k=0.670139,b=0       0.39471       0.39471       0.39471       0.39471      0.39471       0.39471       0.39471       0.39471       0.39471       0.39471       0.39471       0.39471 
    coh49:k=0.680142,b=0       0.40332       0.40332       0.40332       0.40332      0.40332       0.40332       0.40332       0.40332       0.40332       0.40332       0.40332       0.40332 
    coh50:k=0.690144,b=0       0.41202       0.41202       0.41202       0.41202      0.41202       0.41202       0.41202       0.41202       0.41202       0.41202       0.41202       0.41202 
    coh5905:k=59.2574,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 
    coh5906:k=59.2674,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 
    coh5907:k=59.2774,b=0       50.089        50.089        50.089        50.089       50.089        50.089        50.089        50.089        50.089        50.089        50.089        50.089 
    coh5908:k=59.2874,b=0       49.979        49.979        49.979        49.979       49.979        49.979        49.979        49.979        49.979        49.979        49.979        49.979 
    coh5909:k=59.2974,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 
    coh5910:k=59.3074,b=0       49.996        49.996        49.996        49.996       49.996        49.996        49.996        49.996        49.996        49.996        49.996        49.996 
    coh5911:k=59.3174,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 
    coh5912:k=59.3274,b=0       50.013        50.013        50.013        50.013       50.013        50.013        50.013        50.013        50.013        50.013        50.013        50.013 
    coh5913:k=59.3374,b=0       50.021        50.021        50.021        50.021       50.021        50.021        50.021        50.021        50.021        50.021        50.021        50.021 
    coh5914:k=59.3474,b=0        50.03         50.03         50.03         50.03        50.03         50.03         50.03         50.03         50.03         50.03         50.03         50.03 
    coh5915:k=59.3574,b=0       50.038        50.038        50.038        50.038       50.038        50.038        50.038        50.038        50.038        50.038        50.038        50.038 
    coh5916:k=59.3674,b=0       50.047        50.047        50.047        50.047       50.047        50.047        50.047        50.047        50.047        50.047        50.047        50.047 
    coh5917:k=59.3774,b=0       50.055        50.055        50.055        50.055       50.055        50.055        50.055        50.055        50.055        50.055        50.055        50.055 
    coh5918:k=59.3874,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 
    coh5919:k=59.3974,b=0       50.072        50.072        50.072        50.072       50.072        50.072        50.072        50.072        50.072        50.072        50.072        50.072 
    coh5920:k=59.4075,b=0        50.08         50.08         50.08         50.08        50.08         50.08         50.08         50.08         50.08         50.08         50.08         50.08 
    coh5921:k=59.4175,b=0       50.089        50.089        50.089        50.089       50.089        50.089        50.089        50.089        50.089        50.089        50.089        50.089 
    coh5922:k=59.4275,b=0       50.097        50.097        50.097        50.097       50.097        50.097        50.097        50.097        50.097        50.097        50.097        50.097 
    coh5923:k=59.4375,b=0       49.987        49.987        49.987        49.987       49.987        49.987        49.987        49.987        49.987        49.987        49.987        49.987 
    coh5924:k=59.4475,b=0       49.995        49.995        49.995        49.995       49.995        49.995        49.995        49.995        49.995        49.995        49.995        49.995 
    coh5925:k=59.4575,b=0       50.004        50.004        50.004        50.004       50.004        50.004        50.004        50.004        50.004        50.004        50.004        50.004 
    coh5926:k=59.4675,b=0       50.012        50.012        50.012        50.012       50.012        50.012        50.012        50.012        50.012        50.012        50.012        50.012 
    coh5927:k=59.4775,b=0       50.021        50.021        50.021        50.021       50.021        50.021        50.021        50.021        50.021        50.021        50.021        50.021 
    coh5928:k=59.4875,b=0       50.029        50.029        50.029        50.029       50.029        50.029        50.029        50.029        50.029        50.029        50.029        50.029 
    coh5929:k=59.4975,b=0       50.037        50.037        50.037        50.037       50.037        50.037        50.037        50.037        50.037        50.037        50.037        50.037 
    coh5930:k=59.5075,b=0       50.046        50.046        50.046        50.046       50.046        50.046        50.046        50.046        50.046        50.046        50.046        50.046 
    coh5931:k=59.5175,b=0       50.054        50.054        50.054        50.054       50.054        50.054        50.054        50.054        50.054        50.054        50.054        50.054 
    coh5932:k=59.5275,b=0       50.063        50.063        50.063        50.063       50.063        50.063        50.063        50.063        50.063        50.063        50.063        50.063 
    coh5933:k=59.5375,b=0       50.071        50.071        50.071        50.071       50.071        50.071        50.071        50.071        50.071        50.071        50.071        50.071 
    coh5934:k=59.5475,b=0       50.079        50.079        50.079        50.079       50.079        50.079        50.079        50.079        50.079        50.079        50.079        50.079 
    coh5935:k=59.5575,b=0       50.088        50.088        50.088        50.088       50.088        50.088        50.088        50.088        50.088        50.088        50.088        50.088 
    coh5936:k=59.5675,b=0       50.096        50.096        50.096        50.096       50.096        50.096        50.096        50.096        50.096        50.096        50.096        50.096 
    coh5937:k=59.5775,b=0       49.986        49.986        49.986        49.986       49.986        49.986        49.986        49.986        49.986        49.986        49.986        49.986 
    coh5938:k=59.5875,b=0       49.994        49.994        49.994        49.994       49.994        49.994        49.994        49.994        49.994        49.994        49.994        49.994 
    coh5939:k=59.5975,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 
    coh5940:k=59.6075,b=0       50.011        50.011        50.011        50.011       50.011        50.011        50.011        50.011        50.011        50.011        50.011        50.011 
    coh5941:k=59.6175,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 
    coh5942:k=59.6275,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 
    coh5943:k=59.6375,b=0       50.036        50.036        50.036        50.036       50.036        50.036        50.036        50.036        50.036        50.036        50.036        50.036 
    coh5944:k=59.6475,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 
    coh5945:k=59.6575,b=0       50.053        50.053        50.053        50.053       50.053        50.053        50.053        50.053        50.053        50.053        50.053        50.053 
    coh5946:k=59.6675,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 
    coh5947:k=59.6775,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 
    coh5948:k=59.6875,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 
    coh5949:k=59.6975,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 
    coh5950:k=59.7075,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 
    coh5951:k=59.7175,b=0       49.984        49.984        49.984        49.984       49.984        49.984        49.984        49.984        49.984        49.984        49.984        49.984 
    coh5952:k=59.7275,b=0       49.992        49.992        49.992        49.992       49.992        49.992        49.992        49.992        49.992        49.992        49.992        49.992 
    coh5953:k=59.7375,b=0           50            50            50            50           50            50            50            50            50            50            50            50 
    coh5954:k=59.7476,b=0       50.009        50.009        50.009        50.009       50.009        50.009        50.009        50.009        50.009        50.009        50.009        50.009 

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_ipwkz_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_ipwkz_vf_vecsv(param_map, support_map);
Warning: Ignoring extra legend entries. 
Elapsed time is 177.508843 seconds.
Warning: Directory already exists.