Contents

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_akz_set_default_param(it_param_set);

% Simulation Accuracy
param_map('it_w_n') = 1250;
param_map('it_ak_n') = param_map('it_w_n');
param_map('it_z_n') = 15;
param_map('fl_coh_interp_grid_gap') = 0.0125;
param_map('it_c_interp_grid_gap') = 10^-4;

% 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_iwkz_vf_vecsv(param_map, support_map);
Elapsed time is 98.754855 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       416.42     534.64             1             1             1             1             1             1             1             1             1             1             1             1             1             1             1
    iter=2        317.8     4481.7       0.99706       0.99706       0.99706       0.99706       0.99706       0.99706       0.99706       0.99706       0.99706       0.99706       0.99706       0.99706       0.99706       0.99706       0.99706
    iter=3       261.55     1493.2       0.99538       0.99538       0.99538       0.99538       0.99538       0.99538       0.99538       0.99538       0.99538       0.99538       0.99538       0.99538       0.99538       0.99538       0.99538
    iter=4       222.05     748.76        0.9916        0.9916        0.9916        0.9916        0.9916        0.9916        0.9916        0.9916        0.9916        0.9916        0.9916        0.9916        0.9916        0.9916        0.9916
    iter=5       192.67     454.88       0.98573       0.98573       0.98573       0.98573       0.98573       0.98573       0.98573       0.98573       0.98573       0.98573       0.98573       0.98573       0.98573       0.98573       0.98573
    iter=6        170.2     312.35       0.97649       0.97649       0.97649       0.97649       0.97649       0.97649       0.97649       0.97649       0.97649       0.97649       0.97649       0.97649       0.97649       0.97649       0.97649
    iter=7       152.75     237.11       0.96526       0.96526       0.96526       0.96526       0.96526       0.96526       0.96526       0.96526       0.96526       0.96526       0.96526       0.96526       0.96526       0.96526       0.96526
    iter=8       139.04     195.38       0.95854       0.95854       0.95854       0.95854       0.95854       0.95854       0.95854       0.95854       0.95854       0.95854       0.95854       0.95854       0.95854       0.95854       0.95854
    iter=9       128.13     167.84       0.94553       0.94553       0.94553       0.94553       0.94553       0.94553       0.94553       0.94553       0.94553       0.94553       0.94553       0.94553       0.94553       0.94553       0.94553
    iter=10      119.34     158.01       0.94301       0.94301       0.94301       0.94301       0.94301       0.94301       0.94301       0.94301       0.94301       0.94301       0.94301       0.94301       0.94301       0.94301       0.94301
    iter=11      112.17     147.45       0.94752       0.94752       0.94752       0.94752       0.94752       0.94752       0.94752       0.94752       0.94752       0.94752       0.94752       0.94752       0.94752       0.94752       0.94752
    iter=12      106.21     142.43       0.94165       0.94165       0.94165       0.94165       0.94165       0.94165       0.94165       0.94165       0.94165       0.94165       0.94165       0.94165       0.94165       0.94165       0.94165
    iter=13      101.16     156.36       0.94836       0.94836       0.94836       0.94836       0.94836       0.94836       0.94836       0.94836       0.94836       0.94836       0.94836       0.94836       0.94836       0.94836       0.94836
    iter=14      96.781     153.15       0.92643       0.92643       0.92643       0.92643       0.92643       0.92643       0.92643       0.92643       0.92643       0.92643       0.92643       0.92643       0.92643       0.92643       0.92643
    iter=15      92.907     162.04       0.92706       0.92706       0.92706       0.92706       0.92706       0.92706       0.92706       0.92706       0.92706       0.92706       0.92706       0.92706       0.92706       0.92706       0.92706
    iter=16      89.402     145.42       0.92664       0.92664       0.92664       0.92664       0.92664       0.92664       0.92664       0.92664       0.92664       0.92664       0.92664       0.92664       0.92664       0.92664       0.92664
    iter=17      86.172     144.29       0.84939       0.84939       0.84939       0.84939       0.84939       0.84939       0.84939       0.84939       0.84939       0.84939       0.84939       0.84939       0.84939       0.84939       0.84939
    iter=18      83.144     156.35        0.9152        0.9152        0.9152        0.9152        0.9152        0.9152        0.9152        0.9152        0.9152        0.9152        0.9152        0.9152        0.9152        0.9152        0.9152
    iter=19       80.27     149.81       0.93272       0.93272       0.93272       0.93272       0.93272       0.93272       0.93272       0.93272       0.93272       0.93272       0.93272       0.93272       0.93272       0.93272       0.93272
    iter=20      77.514     157.09       0.92433       0.92433       0.92433       0.92433       0.92433       0.92433       0.92433       0.92433       0.92433       0.92433       0.92433       0.92433       0.92433       0.92433       0.92433
    iter=21      74.851     149.88       0.92181       0.92181       0.92181       0.92181       0.92181       0.92181       0.92181       0.92181       0.92181       0.92181       0.92181       0.92181       0.92181       0.92181       0.92181
    iter=22      72.264      153.6       0.92045       0.92045       0.92045       0.92045       0.92045       0.92045       0.92045       0.92045       0.92045       0.92045       0.92045       0.92045       0.92045       0.92045       0.92045
    iter=23      69.743     157.75       0.92118       0.92118       0.92118       0.92118       0.92118       0.92118       0.92118       0.92118       0.92118       0.92118       0.92118       0.92118       0.92118       0.92118       0.92118
    iter=24      67.279     166.67       0.92097       0.92097       0.92097       0.92097       0.92097       0.92097       0.92097       0.92097       0.92097       0.92097       0.92097       0.92097       0.92097       0.92097       0.92097
    iter=25       64.87     169.09        0.9194        0.9194        0.9194        0.9194        0.9194        0.9194        0.9194        0.9194        0.9194        0.9194        0.9194        0.9194        0.9194        0.9194        0.9194
    iter=26      62.512      167.5       0.90554       0.90554       0.90554       0.90554       0.90554       0.90554       0.90554       0.90554       0.90554       0.90554       0.90554       0.90554       0.90554       0.90554       0.90554
    iter=27      60.205     155.52        0.9067        0.9067        0.9067        0.9067        0.9067        0.9067        0.9067        0.9067        0.9067        0.9067        0.9067        0.9067        0.9067        0.9067        0.9067
    iter=28      57.949      170.5       0.89085       0.89085       0.89085       0.89085       0.89085       0.89085       0.89085       0.89085       0.89085       0.89085       0.89085       0.89085       0.89085       0.89085       0.89085
    iter=29      55.744     150.97       0.78243       0.78243       0.78243       0.78243       0.78243       0.78243       0.78243       0.78243       0.78243       0.78243       0.78243       0.78243       0.78243       0.78243       0.78243
    iter=30      53.591     156.61       0.87542       0.87542       0.87542       0.87542       0.87542       0.87542       0.87542       0.87542       0.87542       0.87542       0.87542       0.87542       0.87542       0.87542       0.87542
    iter=31      51.492      152.2       0.87783       0.87783       0.87783       0.87783       0.87783       0.87783       0.87783       0.87783       0.87783       0.87783       0.87783       0.87783       0.87783       0.87783       0.87783
    iter=32      49.446     203.74       0.87615       0.87615       0.87615       0.87615       0.87615       0.87615       0.87615       0.87615       0.87615       0.87615       0.87615       0.87615       0.87615       0.87615       0.87615
    iter=33      47.456     166.26       0.87794       0.87794       0.87794       0.87794       0.87794       0.87794       0.87794       0.87794       0.87794       0.87794       0.87794       0.87794       0.87794       0.87794       0.87794
    iter=34      45.521     167.36       0.88046       0.88046       0.88046       0.88046       0.88046       0.88046       0.88046       0.88046       0.88046       0.88046       0.88046       0.88046       0.88046       0.88046       0.88046
    iter=35      43.642      177.7       0.87049       0.87049       0.87049       0.87049       0.87049       0.87049       0.87049       0.87049       0.87049       0.87049       0.87049       0.87049       0.87049       0.87049       0.87049
    iter=36      41.819     168.08       0.86839       0.86839       0.86839       0.86839       0.86839       0.86839       0.86839       0.86839       0.86839       0.86839       0.86839       0.86839       0.86839       0.86839       0.86839
    iter=37      40.053     164.74       0.86629       0.86629       0.86629       0.86629       0.86629       0.86629       0.86629       0.86629       0.86629       0.86629       0.86629       0.86629       0.86629       0.86629       0.86629
    iter=38      38.343     190.25       0.84887       0.84887       0.84887       0.84887       0.84887       0.84887       0.84887       0.84887       0.84887       0.84887       0.84887       0.84887       0.84887       0.84887       0.84887
    iter=39      36.689      167.3       0.84824       0.84824       0.84824       0.84824       0.84824       0.84824       0.84824       0.84824       0.84824       0.84824       0.84824       0.84824       0.84824       0.84824       0.84824
    iter=40      35.091     169.98       0.84037       0.84037       0.84037       0.84037       0.84037       0.84037       0.84037       0.84037       0.84037       0.84037       0.84037       0.84037       0.84037       0.84037       0.84037
    iter=41      33.549     170.11       0.83428       0.83428       0.83428       0.83428       0.83428       0.83428       0.83428       0.83428       0.83428       0.83428       0.83428       0.83428       0.83428       0.83428       0.83428
    iter=42      32.061     155.44       0.81791       0.81791       0.81791       0.81791       0.81791       0.81791       0.81791       0.81791       0.81791       0.81791       0.81791       0.81791       0.81791       0.81791       0.81791
    iter=43      30.626     156.55       0.78967       0.78967       0.78967       0.78967       0.78967       0.78967       0.78967       0.78967       0.78967       0.78967       0.78967       0.78967       0.78967       0.78967       0.78967
    iter=44      29.245     165.31       0.78023       0.78023       0.78023       0.78023       0.78023       0.78023       0.78023       0.78023       0.78023       0.78023       0.78023       0.78023       0.78023       0.78023       0.78023
    iter=45      27.915     155.23        0.7755        0.7755        0.7755        0.7755        0.7755        0.7755        0.7755        0.7755        0.7755        0.7755        0.7755        0.7755        0.7755        0.7755        0.7755
    iter=46      26.636     146.39       0.63927       0.63927       0.63927       0.63927       0.63927       0.63927       0.63927       0.63927       0.63927       0.63927       0.63927       0.63927       0.63927       0.63927       0.63927
    iter=47      25.407     153.25       0.73363       0.73363       0.73363       0.73363       0.73363       0.73363       0.73363       0.73363       0.73363       0.73363       0.73363       0.73363       0.73363       0.73363       0.73363
    iter=48      24.226     139.82       0.71411       0.71411       0.71411       0.71411       0.71411       0.71411       0.71411       0.71411       0.71411       0.71411       0.71411       0.71411       0.71411       0.71411       0.71411
    iter=49      23.092     142.45       0.69962       0.69962       0.69962       0.69962       0.69962       0.69962       0.69962       0.69962       0.69962       0.69962       0.69962       0.69962       0.69962       0.69962       0.69962
    iter=50      22.004     129.68       0.68419       0.68419       0.68419       0.68419       0.68419       0.68419       0.68419       0.68419       0.68419       0.68419       0.68419       0.68419       0.68419       0.68419       0.68419
    iter=92      2.2948      33.58       0.03327       0.03327       0.03327       0.03327       0.03327       0.03327       0.03327       0.03327       0.03327       0.03327       0.03327       0.03327       0.03327       0.03327       0.03327
    iter=93      2.1639     25.724      0.024349      0.024349      0.024349      0.024349      0.024349      0.024349      0.024349      0.024349      0.024349      0.024349      0.024349      0.024349      0.024349      0.024349      0.024349
    iter=94      2.0401      23.08      0.020256      0.020256      0.020256      0.020256      0.020256      0.020256      0.020256      0.020256      0.020256      0.020256      0.020256      0.020256      0.020256      0.020256      0.020256
    iter=95       1.923     22.252      0.024454      0.024454      0.024454      0.024454      0.024454      0.024454      0.024454      0.024454      0.024454      0.024454      0.024454      0.024454      0.024454      0.024454      0.024454
    iter=96      1.8123      31.47      0.022145      0.022145      0.022145      0.022145      0.022145      0.022145      0.022145      0.022145      0.022145      0.022145      0.022145      0.022145      0.022145      0.022145      0.022145
    iter=97      1.7076      16.94      0.013329      0.013329      0.013329      0.013329      0.013329      0.013329      0.013329      0.013329      0.013329      0.013329      0.013329      0.013329      0.013329      0.013329      0.013329
    iter=98      1.6087     18.358      0.014589      0.014589      0.014589      0.014589      0.014589      0.014589      0.014589      0.014589      0.014589      0.014589      0.014589      0.014589      0.014589      0.014589      0.014589
    iter=99      1.5152     27.625      0.012385      0.012385      0.012385      0.012385      0.012385      0.012385      0.012385      0.012385      0.012385      0.012385      0.012385      0.012385      0.012385      0.012385      0.012385
    iter=100     1.4269     16.822      0.013959      0.013959      0.013959      0.013959      0.013959      0.013959      0.013959      0.013959      0.013959      0.013959      0.013959      0.013959      0.013959      0.013959      0.013959
    iter=101     1.3435     17.409     0.0051427     0.0051427     0.0051427     0.0051427     0.0051427     0.0051427     0.0051427     0.0051427     0.0051427     0.0051427     0.0051427     0.0051427     0.0051427     0.0051427     0.0051427
    iter=102     1.2648     19.305      0.010915      0.010915      0.010915      0.010915      0.010915      0.010915      0.010915      0.010915      0.010915      0.010915      0.010915      0.010915      0.010915      0.010915      0.010915
    iter=103     1.1905     8.4666     0.0075567     0.0075567     0.0075567     0.0075567     0.0075567     0.0075567     0.0075567     0.0075567     0.0075567     0.0075567     0.0075567     0.0075567     0.0075567     0.0075567     0.0075567
    iter=104     1.1204     13.866     0.0058774     0.0058774     0.0058774     0.0058774     0.0058774     0.0058774     0.0058774     0.0058774     0.0058774     0.0058774     0.0058774     0.0058774     0.0058774     0.0058774     0.0058774
    iter=105     1.0542     6.5779     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387
    iter=106    0.99189     7.8047     0.0038833     0.0038833     0.0038833     0.0038833     0.0038833     0.0038833     0.0038833     0.0038833     0.0038833     0.0038833     0.0038833     0.0038833     0.0038833     0.0038833     0.0038833
    iter=107    0.93312      11.09     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387     0.0029387
    iter=108    0.87774     16.166     0.0025189     0.0025189     0.0025189     0.0025189     0.0025189     0.0025189     0.0025189     0.0025189     0.0025189     0.0025189     0.0025189     0.0025189     0.0025189     0.0025189     0.0025189
    iter=109    0.82556     8.0801     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694
    iter=110    0.77641     6.2249     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694     0.0014694
    iter=111    0.73013     6.2056    0.00094458    0.00094458    0.00094458    0.00094458    0.00094458    0.00094458    0.00094458    0.00094458    0.00094458    0.00094458    0.00094458    0.00094458    0.00094458    0.00094458    0.00094458
    iter=112    0.68656     6.2017    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963
    iter=113    0.64556    0.57863    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963
    iter=114    0.60697          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=115    0.57066     1.2403    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963    0.00083963
    iter=116     0.5365          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=117    0.50438          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=118    0.47416          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=119    0.44574          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=120    0.41902          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=121     0.3939          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=122    0.37028          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=123    0.34807          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=124    0.32719          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=125    0.30757          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=126    0.28912          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=127    0.27177          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=128    0.25547          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=129    0.24014          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=130    0.22573          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=131    0.21219          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=132    0.19946          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=133    0.18749          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=134    0.17624          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=135    0.16567          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=136    0.15573          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=137    0.14638          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=138     0.1376          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=139    0.12934          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=140    0.12158          0             0             0             0             0             0             0             0             0             0             0             0             0             0             0             0
    iter=141    0.11429          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.196       -41.196       -41.196       -41.196      -41.196       -41.196       -41.196       -41.196       -41.196       -41.196       -41.196       -41.196  
    coh2:k=0.212504,b=0       -41.062       -41.062       -41.062       -41.062      -41.062       -41.062       -41.062       -41.062       -41.062       -41.062       -41.062       -41.062  
    coh3:k=0.225008,b=0        -40.94        -40.94        -40.94        -40.94       -40.94        -40.94        -40.94        -40.94        -40.94        -40.94        -40.94        -40.94  
    coh4:k=0.237512,b=0       -40.828       -40.828       -40.828       -40.828      -40.828       -40.828       -40.828       -40.828       -40.828       -40.828       -40.828       -40.828  
    coh5:k=0.250015,b=0       -40.711       -40.711       -40.711       -40.711      -40.711       -40.711       -40.711       -40.711       -40.711       -40.711       -40.711       -40.711  
    coh6:k=0.262519,b=0       -40.586       -40.586       -40.586       -40.586      -40.586       -40.586       -40.586       -40.586       -40.586       -40.586       -40.586       -40.586  
    coh7:k=0.275023,b=0       -40.472       -40.472       -40.472       -40.472      -40.472       -40.472       -40.472       -40.472       -40.472       -40.472       -40.472       -40.472  
    coh8:k=0.287527,b=0       -40.366       -40.366       -40.366       -40.366      -40.366       -40.366       -40.366       -40.366       -40.366       -40.366       -40.366       -40.366  
    coh9:k=0.300031,b=0       -40.251       -40.251       -40.251       -40.251      -40.251       -40.251       -40.251       -40.251       -40.251       -40.251       -40.251       -40.251  
    coh10:k=0.312535,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  
    coh11:k=0.325039,b=0      -40.027       -40.027       -40.027       -40.027      -40.027       -40.027       -40.027       -40.027       -40.027       -40.027       -40.027       -40.027  
    coh12:k=0.337542,b=0      -39.928       -39.928       -39.928       -39.928      -39.928       -39.928       -39.928       -39.928       -39.928       -39.928       -39.928       -39.928  
    coh13:k=0.350046,b=0      -39.811       -39.811       -39.811       -39.811      -39.811       -39.811       -39.811       -39.811       -39.811       -39.811       -39.811       -39.811  
    coh14:k=0.36255,b=0       -39.702       -39.702       -39.702       -39.702      -39.702       -39.702       -39.702       -39.702       -39.702       -39.702       -39.702       -39.702  
    coh15:k=0.375054,b=0      -39.601       -39.601       -39.601       -39.601      -39.601       -39.601       -39.601       -39.601       -39.601       -39.601       -39.601       -39.601  
    coh16:k=0.387558,b=0      -39.501       -39.501       -39.501       -39.501      -39.501       -39.501       -39.501       -39.501       -39.501       -39.501       -39.501       -39.501  
    coh17:k=0.400062,b=0       -39.39        -39.39        -39.39        -39.39       -39.39        -39.39        -39.39        -39.39        -39.39        -39.39        -39.39        -39.39  
    coh18:k=0.412566,b=0      -39.288       -39.288       -39.288       -39.288      -39.288       -39.288       -39.288       -39.288       -39.288       -39.288       -39.288       -39.288  
    coh19:k=0.425069,b=0      -39.193       -39.193       -39.193       -39.193      -39.193       -39.193       -39.193       -39.193       -39.193       -39.193       -39.193       -39.193  
    coh20:k=0.437573,b=0      -39.091       -39.091       -39.091       -39.091      -39.091       -39.091       -39.091       -39.091       -39.091       -39.091       -39.091       -39.091  
    coh21:k=0.450077,b=0      -38.987       -38.987       -38.987       -38.987      -38.987       -38.987       -38.987       -38.987       -38.987       -38.987       -38.987       -38.987  
    coh22:k=0.462581,b=0      -38.891       -38.891       -38.891       -38.891      -38.891       -38.891       -38.891       -38.891       -38.891       -38.891       -38.891       -38.891  
    coh23:k=0.475085,b=0      -38.801       -38.801       -38.801       -38.801      -38.801       -38.801       -38.801       -38.801       -38.801       -38.801       -38.801       -38.801  
    coh24:k=0.487589,b=0      -38.698       -38.698       -38.698       -38.698      -38.698       -38.698       -38.698       -38.698       -38.698       -38.698       -38.698       -38.698  
    coh25:k=0.500092,b=0        -38.6         -38.6         -38.6         -38.6        -38.6         -38.6         -38.6         -38.6         -38.6         -38.6         -38.6         -38.6  
    coh26:k=0.512596,b=0      -38.509       -38.509       -38.509       -38.509      -38.509       -38.509       -38.509       -38.509       -38.509       -38.509       -38.509       -38.509  
    coh27:k=0.5251,b=0        -38.418       -38.418       -38.418       -38.418      -38.418       -38.418       -38.418       -38.418       -38.418       -38.418       -38.418       -38.418  
    coh28:k=0.537604,b=0      -38.319       -38.319       -38.319       -38.319      -38.319       -38.319       -38.319       -38.319       -38.319       -38.319       -38.319       -38.319  
    coh29:k=0.550108,b=0      -38.227       -38.227       -38.227       -38.227      -38.227       -38.227       -38.227       -38.227       -38.227       -38.227       -38.227       -38.227  
    coh30:k=0.562612,b=0       -38.14        -38.14        -38.14        -38.14       -38.14        -38.14        -38.14        -38.14        -38.14        -38.14        -38.14        -38.14  
    coh31:k=0.575116,b=0      -38.046       -38.046       -38.046       -38.046      -38.046       -38.046       -38.046       -38.046       -38.046       -38.046       -38.046       -38.046  
    coh32:k=0.587619,b=0      -37.952       -37.952       -37.952       -37.952      -37.952       -37.952       -37.952       -37.952       -37.952       -37.952       -37.952       -37.952  
    coh33:k=0.600123,b=0      -37.865       -37.865       -37.865       -37.865      -37.865       -37.865       -37.865       -37.865       -37.865       -37.865       -37.865       -37.865  
    coh34:k=0.612627,b=0      -37.781       -37.781       -37.781       -37.781      -37.781       -37.781       -37.781       -37.781       -37.781       -37.781       -37.781       -37.781  
    coh35:k=0.625131,b=0      -37.686       -37.686       -37.686       -37.686      -37.686       -37.686       -37.686       -37.686       -37.686       -37.686       -37.686       -37.686  
    coh36:k=0.637635,b=0      -37.597       -37.597       -37.597       -37.597      -37.597       -37.597       -37.597       -37.597       -37.597       -37.597       -37.597       -37.597  
    coh37:k=0.650139,b=0      -37.515       -37.515       -37.515       -37.515      -37.515       -37.515       -37.515       -37.515       -37.515       -37.515       -37.515       -37.515  
    coh38:k=0.662643,b=0      -37.427       -37.427       -37.427       -37.427      -37.427       -37.427       -37.427       -37.427       -37.427       -37.427       -37.427       -37.427  
    coh39:k=0.675146,b=0      -37.337       -37.337       -37.337       -37.337      -37.337       -37.337       -37.337       -37.337       -37.337       -37.337       -37.337       -37.337  
    coh40:k=0.68765,b=0       -37.253       -37.253       -37.253       -37.253      -37.253       -37.253       -37.253       -37.253       -37.253       -37.253       -37.253       -37.253  
    coh41:k=0.700154,b=0      -37.174       -37.174       -37.174       -37.174      -37.174       -37.174       -37.174       -37.174       -37.174       -37.174       -37.174       -37.174  
    coh42:k=0.712658,b=0      -37.083       -37.083       -37.083       -37.083      -37.083       -37.083       -37.083       -37.083       -37.083       -37.083       -37.083       -37.083  
    coh43:k=0.725162,b=0      -36.998       -36.998       -36.998       -36.998      -36.998       -36.998       -36.998       -36.998       -36.998       -36.998       -36.998       -36.998  
    coh44:k=0.737666,b=0      -36.918       -36.918       -36.918       -36.918      -36.918       -36.918       -36.918       -36.918       -36.918       -36.918       -36.918       -36.918  
    coh45:k=0.75017,b=0       -36.835       -36.835       -36.835       -36.835      -36.835       -36.835       -36.835       -36.835       -36.835       -36.835       -36.835       -36.835  
    coh46:k=0.762673,b=0      -36.749       -36.749       -36.749       -36.749      -36.749       -36.749       -36.749       -36.749       -36.749       -36.749       -36.749       -36.749  
    coh47:k=0.775177,b=0      -36.668       -36.668       -36.668       -36.668      -36.668       -36.668       -36.668       -36.668       -36.668       -36.668       -36.668       -36.668  
    coh48:k=0.787681,b=0      -36.592       -36.592       -36.592       -36.592      -36.592       -36.592       -36.592       -36.592       -36.592       -36.592       -36.592       -36.592  
    coh49:k=0.800185,b=0      -36.505       -36.505       -36.505       -36.505      -36.505       -36.505       -36.505       -36.505       -36.505       -36.505       -36.505       -36.505  
    coh50:k=0.812689,b=0      -36.423       -36.423       -36.423       -36.423      -36.423       -36.423       -36.423       -36.423       -36.423       -36.423       -36.423       -36.423  
    coh4715:k=59.1432,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  
    coh4716:k=59.1557,b=0      9.7224        9.7224        9.7224        9.7224       9.7224        9.7224        9.7224        9.7224        9.7224        9.7224        9.7224        9.7224  
    coh4717:k=59.1682,b=0      9.7228        9.7228        9.7228        9.7228       9.7228        9.7228        9.7228        9.7228        9.7228        9.7228        9.7228        9.7228  
    coh4718:k=59.1807,b=0      9.7233        9.7233        9.7233        9.7233       9.7233        9.7233        9.7233        9.7233        9.7233        9.7233        9.7233        9.7233  
    coh4719:k=59.1932,b=0      9.7237        9.7237        9.7237        9.7237       9.7237        9.7237        9.7237        9.7237        9.7237        9.7237        9.7237        9.7237  
    coh4720:k=59.2057,b=0      9.7242        9.7242        9.7242        9.7242       9.7242        9.7242        9.7242        9.7242        9.7242        9.7242        9.7242        9.7242  
    coh4721:k=59.2182,b=0      9.7246        9.7246        9.7246        9.7246       9.7246        9.7246        9.7246        9.7246        9.7246        9.7246        9.7246        9.7246  
    coh4722:k=59.2307,b=0      9.7251        9.7251        9.7251        9.7251       9.7251        9.7251        9.7251        9.7251        9.7251        9.7251        9.7251        9.7251  
    coh4723:k=59.2432,b=0      9.7255        9.7255        9.7255        9.7255       9.7255        9.7255        9.7255        9.7255        9.7255        9.7255        9.7255        9.7255  
    coh4724:k=59.2557,b=0       9.726         9.726         9.726         9.726        9.726         9.726         9.726         9.726         9.726         9.726         9.726         9.726  
    coh4725:k=59.2682,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  
    coh4726:k=59.2807,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  
    coh4727:k=59.2932,b=0      9.7273        9.7273        9.7273        9.7273       9.7273        9.7273        9.7273        9.7273        9.7273        9.7273        9.7273        9.7273  
    coh4728:k=59.3057,b=0      9.7277        9.7277        9.7277        9.7277       9.7277        9.7277        9.7277        9.7277        9.7277        9.7277        9.7277        9.7277  
    coh4729:k=59.3182,b=0      9.7282        9.7282        9.7282        9.7282       9.7282        9.7282        9.7282        9.7282        9.7282        9.7282        9.7282        9.7282  
    coh4730:k=59.3307,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  
    coh4731:k=59.3432,b=0      9.7291        9.7291        9.7291        9.7291       9.7291        9.7291        9.7291        9.7291        9.7291        9.7291        9.7291        9.7291  
    coh4732:k=59.3557,b=0      9.7295        9.7295        9.7295        9.7295       9.7295        9.7295        9.7295        9.7295        9.7295        9.7295        9.7295        9.7295  
    coh4733:k=59.3682,b=0      9.7299        9.7299        9.7299        9.7299       9.7299        9.7299        9.7299        9.7299        9.7299        9.7299        9.7299        9.7299  
    coh4734:k=59.3807,b=0      9.7304        9.7304        9.7304        9.7304       9.7304        9.7304        9.7304        9.7304        9.7304        9.7304        9.7304        9.7304  
    coh4735:k=59.3932,b=0      9.7308        9.7308        9.7308        9.7308       9.7308        9.7308        9.7308        9.7308        9.7308        9.7308        9.7308        9.7308  
    coh4736:k=59.4057,b=0      9.7312        9.7312        9.7312        9.7312       9.7312        9.7312        9.7312        9.7312        9.7312        9.7312        9.7312        9.7312  
    coh4737:k=59.4183,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  
    coh4738:k=59.4308,b=0      9.7321        9.7321        9.7321        9.7321       9.7321        9.7321        9.7321        9.7321        9.7321        9.7321        9.7321        9.7321  
    coh4739:k=59.4433,b=0      9.7325        9.7325        9.7325        9.7325       9.7325        9.7325        9.7325        9.7325        9.7325        9.7325        9.7325        9.7325  
    coh4740:k=59.4558,b=0       9.733         9.733         9.733         9.733        9.733         9.733         9.733         9.733         9.733         9.733         9.733         9.733  
    coh4741:k=59.4683,b=0      9.7334        9.7334        9.7334        9.7334       9.7334        9.7334        9.7334        9.7334        9.7334        9.7334        9.7334        9.7334  
    coh4742:k=59.4808,b=0      9.7338        9.7338        9.7338        9.7338       9.7338        9.7338        9.7338        9.7338        9.7338        9.7338        9.7338        9.7338  
    coh4743:k=59.4933,b=0      9.7342        9.7342        9.7342        9.7342       9.7342        9.7342        9.7342        9.7342        9.7342        9.7342        9.7342        9.7342  
    coh4744:k=59.5058,b=0      9.7347        9.7347        9.7347        9.7347       9.7347        9.7347        9.7347        9.7347        9.7347        9.7347        9.7347        9.7347  
    coh4745:k=59.5183,b=0      9.7351        9.7351        9.7351        9.7351       9.7351        9.7351        9.7351        9.7351        9.7351        9.7351        9.7351        9.7351  
    coh4746:k=59.5308,b=0      9.7355        9.7355        9.7355        9.7355       9.7355        9.7355        9.7355        9.7355        9.7355        9.7355        9.7355        9.7355  
    coh4747:k=59.5433,b=0       9.736         9.736         9.736         9.736        9.736         9.736         9.736         9.736         9.736         9.736         9.736         9.736  
    coh4748:k=59.5558,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  
    coh4749:k=59.5683,b=0      9.7368        9.7368        9.7368        9.7368       9.7368        9.7368        9.7368        9.7368        9.7368        9.7368        9.7368        9.7368  
    coh4750:k=59.5808,b=0      9.7372        9.7372        9.7372        9.7372       9.7372        9.7372        9.7372        9.7372        9.7372        9.7372        9.7372        9.7372  
    coh4751:k=59.5933,b=0      9.7376        9.7376        9.7376        9.7376       9.7376        9.7376        9.7376        9.7376        9.7376        9.7376        9.7376        9.7376  
    coh4752:k=59.6058,b=0      9.7381        9.7381        9.7381        9.7381       9.7381        9.7381        9.7381        9.7381        9.7381        9.7381        9.7381        9.7381  
    coh4753:k=59.6183,b=0      9.7385        9.7385        9.7385        9.7385       9.7385        9.7385        9.7385        9.7385        9.7385        9.7385        9.7385        9.7385  
    coh4754:k=59.6308,b=0      9.7389        9.7389        9.7389        9.7389       9.7389        9.7389        9.7389        9.7389        9.7389        9.7389        9.7389        9.7389  
    coh4755:k=59.6433,b=0      9.7393        9.7393        9.7393        9.7393       9.7393        9.7393        9.7393        9.7393        9.7393        9.7393        9.7393        9.7393  
    coh4756:k=59.6558,b=0      9.7397        9.7397        9.7397        9.7397       9.7397        9.7397        9.7397        9.7397        9.7397        9.7397        9.7397        9.7397  
    coh4757:k=59.6683,b=0      9.7402        9.7402        9.7402        9.7402       9.7402        9.7402        9.7402        9.7402        9.7402        9.7402        9.7402        9.7402  
    coh4758:k=59.6808,b=0      9.7406        9.7406        9.7406        9.7406       9.7406        9.7406        9.7406        9.7406        9.7406        9.7406        9.7406        9.7406  
    coh4759:k=59.6933,b=0       9.741         9.741         9.741         9.741        9.741         9.741         9.741         9.741         9.741         9.741         9.741         9.741  
    coh4760:k=59.7058,b=0      9.7414        9.7414        9.7414        9.7414       9.7414        9.7414        9.7414        9.7414        9.7414        9.7414        9.7414        9.7414  
    coh4761:k=59.7183,b=0      9.7418        9.7418        9.7418        9.7418       9.7418        9.7418        9.7418        9.7418        9.7418        9.7418        9.7418        9.7418  
    coh4762:k=59.7308,b=0      9.7422        9.7422        9.7422        9.7422       9.7422        9.7422        9.7422        9.7422        9.7422        9.7422        9.7422        9.7422  
    coh4763:k=59.7434,b=0      9.7426        9.7426        9.7426        9.7426       9.7426        9.7426        9.7426        9.7426        9.7426        9.7426        9.7426        9.7426  
    coh4764:k=59.7559,b=0       9.743         9.743         9.743         9.743        9.743         9.743         9.743         9.743         9.743         9.743         9.743         9.743  

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                   0             0             0            0             0             0             0             0             0             0             0             0 
    coh2:k=0.212504,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh3:k=0.225008,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh4:k=0.237512,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh5:k=0.250015,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh6:k=0.262519,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh7:k=0.275023,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh8:k=0.287527,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh9:k=0.300031,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh10:k=0.312535,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh11:k=0.325039,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh12:k=0.337542,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh13:k=0.350046,b=0      0.040032      0.040032      0.040032     0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032 
    coh14:k=0.36255,b=0       0.040032      0.040032      0.040032     0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032 
    coh15:k=0.375054,b=0      0.040032      0.040032      0.040032     0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032 
    coh16:k=0.387558,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh17:k=0.400062,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh18:k=0.412566,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh19:k=0.425069,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh20:k=0.437573,b=0        0.1201        0.1201        0.1201       0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201 
    coh21:k=0.450077,b=0        0.1201        0.1201        0.1201       0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201 
    coh22:k=0.462581,b=0        0.1201        0.1201        0.1201       0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201 
    coh23:k=0.475085,b=0        0.1201        0.1201        0.1201       0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201 
    coh24:k=0.487589,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh25:k=0.500092,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh26:k=0.512596,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh27:k=0.5251,b=0               0             0             0            0             0             0             0             0             0             0             0             0 
    coh28:k=0.537604,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh29:k=0.550108,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh30:k=0.562612,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh31:k=0.575116,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh32:k=0.587619,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh33:k=0.600123,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh34:k=0.612627,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh35:k=0.625131,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh36:k=0.637635,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh37:k=0.650139,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh38:k=0.662643,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh39:k=0.675146,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh40:k=0.68765,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh41:k=0.700154,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh42:k=0.712658,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh43:k=0.725162,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh44:k=0.737666,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh45:k=0.75017,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh46:k=0.762673,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh47:k=0.775177,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh48:k=0.787681,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh49:k=0.800185,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh50:k=0.812689,b=0             0             0             0            0             0             0             0             0             0             0             0             0 
    coh4715:k=59.1432,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4716:k=59.1557,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4717:k=59.1682,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4718:k=59.1807,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4719:k=59.1932,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4720:k=59.2057,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4721:k=59.2182,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4722:k=59.2307,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4723:k=59.2432,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4724:k=59.2557,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4725:k=59.2682,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4726:k=59.2807,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4727:k=59.2932,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4728:k=59.3057,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4729:k=59.3182,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4730:k=59.3307,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4731:k=59.3432,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4732:k=59.3557,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4733:k=59.3682,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4734:k=59.3807,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4735:k=59.3932,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4736:k=59.4057,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4737:k=59.4183,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4738:k=59.4308,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4739:k=59.4433,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4740:k=59.4558,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4741:k=59.4683,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4742:k=59.4808,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4743:k=59.4933,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4744:k=59.5058,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4745:k=59.5183,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4746:k=59.5308,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4747:k=59.5433,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4748:k=59.5558,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4749:k=59.5683,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4750:k=59.5808,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4751:k=59.5933,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4752:k=59.6058,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4753:k=59.6183,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4754:k=59.6308,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4755:k=59.6433,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4756:k=59.6558,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4757:k=59.6683,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4758:k=59.6808,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4759:k=59.6933,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4760:k=59.7058,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4761:k=59.7183,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4762:k=59.7308,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4763:k=59.7434,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 
    coh4764:k=59.7559,b=0       35.588        35.588        35.588       35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588        35.588 

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             0             0            0             0             0             0             0             0             0             0             0 
    coh2:k=0.212504,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh3:k=0.225008,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh4:k=0.237512,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh5:k=0.250015,b=0       0.040032      0.040032      0.040032     0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032 
    coh6:k=0.262519,b=0       0.040032      0.040032      0.040032     0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032 
    coh7:k=0.275023,b=0       0.040032      0.040032      0.040032     0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032 
    coh8:k=0.287527,b=0       0.040032      0.040032      0.040032     0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032 
    coh9:k=0.300031,b=0       0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh10:k=0.312535,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh11:k=0.325039,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh12:k=0.337542,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh13:k=0.350046,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh14:k=0.36255,b=0       0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh15:k=0.375054,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh16:k=0.387558,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh17:k=0.400062,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh18:k=0.412566,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh19:k=0.425069,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh20:k=0.437573,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh21:k=0.450077,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh22:k=0.462581,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh23:k=0.475085,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh24:k=0.487589,b=0       0.24019       0.24019       0.24019      0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019 
    coh25:k=0.500092,b=0       0.24019       0.24019       0.24019      0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019 
    coh26:k=0.512596,b=0       0.24019       0.24019       0.24019      0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019 
    coh27:k=0.5251,b=0         0.28022       0.28022       0.28022      0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022 
    coh28:k=0.537604,b=0       0.28022       0.28022       0.28022      0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022 
    coh29:k=0.550108,b=0       0.28022       0.28022       0.28022      0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022 
    coh30:k=0.562612,b=0       0.28022       0.28022       0.28022      0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022 
    coh31:k=0.575116,b=0       0.32026       0.32026       0.32026      0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026 
    coh32:k=0.587619,b=0       0.32026       0.32026       0.32026      0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026 
    coh33:k=0.600123,b=0       0.32026       0.32026       0.32026      0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026 
    coh34:k=0.612627,b=0       0.36029       0.36029       0.36029      0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029 
    coh35:k=0.625131,b=0       0.36029       0.36029       0.36029      0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029 
    coh36:k=0.637635,b=0       0.36029       0.36029       0.36029      0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029 
    coh37:k=0.650139,b=0       0.36029       0.36029       0.36029      0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029 
    coh38:k=0.662643,b=0       0.40032       0.40032       0.40032      0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032 
    coh39:k=0.675146,b=0       0.40032       0.40032       0.40032      0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032 
    coh40:k=0.68765,b=0        0.40032       0.40032       0.40032      0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032 
    coh41:k=0.700154,b=0       0.44035       0.44035       0.44035      0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035 
    coh42:k=0.712658,b=0       0.44035       0.44035       0.44035      0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035 
    coh43:k=0.725162,b=0       0.44035       0.44035       0.44035      0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035 
    coh44:k=0.737666,b=0       0.44035       0.44035       0.44035      0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035 
    coh45:k=0.75017,b=0        0.48038       0.48038       0.48038      0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038 
    coh46:k=0.762673,b=0       0.48038       0.48038       0.48038      0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038 
    coh47:k=0.775177,b=0       0.48038       0.48038       0.48038      0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038 
    coh48:k=0.787681,b=0       0.48038       0.48038       0.48038      0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038 
    coh49:k=0.800185,b=0       0.52042       0.52042       0.52042      0.52042       0.52042       0.52042       0.52042       0.52042       0.52042       0.52042       0.52042       0.52042 
    coh50:k=0.812689,b=0       0.52042       0.52042       0.52042      0.52042       0.52042       0.52042       0.52042       0.52042       0.52042       0.52042       0.52042       0.52042 
    coh4715:k=59.1432,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4716:k=59.1557,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4717:k=59.1682,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4718:k=59.1807,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4719:k=59.1932,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4720:k=59.2057,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4721:k=59.2182,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4722:k=59.2307,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4723:k=59.2432,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4724:k=59.2557,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4725:k=59.2682,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4726:k=59.2807,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4727:k=59.2932,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4728:k=59.3057,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4729:k=59.3182,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4730:k=59.3307,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4731:k=59.3432,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4732:k=59.3557,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4733:k=59.3682,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4734:k=59.3807,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4735:k=59.3932,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4736:k=59.4057,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4737:k=59.4183,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4738:k=59.4308,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4739:k=59.4433,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4740:k=59.4558,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4741:k=59.4683,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4742:k=59.4808,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4743:k=59.4933,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4744:k=59.5058,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4745:k=59.5183,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4746:k=59.5308,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4747:k=59.5433,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4748:k=59.5558,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4749:k=59.5683,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4750:k=59.5808,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4751:k=59.5933,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4752:k=59.6058,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4753:k=59.6183,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4754:k=59.6308,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4755:k=59.6433,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4756:k=59.6558,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4757:k=59.6683,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4758:k=59.6808,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4759:k=59.6933,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4760:k=59.7058,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4761:k=59.7183,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4762:k=59.7308,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4763:k=59.7434,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 
    coh4764:k=59.7559,b=0       14.412        14.412        14.412       14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412        14.412 

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             0             0            0             0             0             0             0             0             0             0             0 
    coh2:k=0.212504,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh3:k=0.225008,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh4:k=0.237512,b=0              0             0             0            0             0             0             0             0             0             0             0             0 
    coh5:k=0.250015,b=0       0.040032      0.040032      0.040032     0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032 
    coh6:k=0.262519,b=0       0.040032      0.040032      0.040032     0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032 
    coh7:k=0.275023,b=0       0.040032      0.040032      0.040032     0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032 
    coh8:k=0.287527,b=0       0.040032      0.040032      0.040032     0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032      0.040032 
    coh9:k=0.300031,b=0       0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh10:k=0.312535,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh11:k=0.325039,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh12:k=0.337542,b=0      0.080064      0.080064      0.080064     0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064      0.080064 
    coh13:k=0.350046,b=0        0.1201        0.1201        0.1201       0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201 
    coh14:k=0.36255,b=0         0.1201        0.1201        0.1201       0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201 
    coh15:k=0.375054,b=0        0.1201        0.1201        0.1201       0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201        0.1201 
    coh16:k=0.387558,b=0       0.16013       0.16013       0.16013      0.16013       0.16013       0.16013       0.16013       0.16013       0.16013       0.16013       0.16013       0.16013 
    coh17:k=0.400062,b=0       0.16013       0.16013       0.16013      0.16013       0.16013       0.16013       0.16013       0.16013       0.16013       0.16013       0.16013       0.16013 
    coh18:k=0.412566,b=0       0.16013       0.16013       0.16013      0.16013       0.16013       0.16013       0.16013       0.16013       0.16013       0.16013       0.16013       0.16013 
    coh19:k=0.425069,b=0       0.16013       0.16013       0.16013      0.16013       0.16013       0.16013       0.16013       0.16013       0.16013       0.16013       0.16013       0.16013 
    coh20:k=0.437573,b=0       0.20016       0.20016       0.20016      0.20016       0.20016       0.20016       0.20016       0.20016       0.20016       0.20016       0.20016       0.20016 
    coh21:k=0.450077,b=0       0.20016       0.20016       0.20016      0.20016       0.20016       0.20016       0.20016       0.20016       0.20016       0.20016       0.20016       0.20016 
    coh22:k=0.462581,b=0       0.20016       0.20016       0.20016      0.20016       0.20016       0.20016       0.20016       0.20016       0.20016       0.20016       0.20016       0.20016 
    coh23:k=0.475085,b=0       0.20016       0.20016       0.20016      0.20016       0.20016       0.20016       0.20016       0.20016       0.20016       0.20016       0.20016       0.20016 
    coh24:k=0.487589,b=0       0.24019       0.24019       0.24019      0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019 
    coh25:k=0.500092,b=0       0.24019       0.24019       0.24019      0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019 
    coh26:k=0.512596,b=0       0.24019       0.24019       0.24019      0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019       0.24019 
    coh27:k=0.5251,b=0         0.28022       0.28022       0.28022      0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022 
    coh28:k=0.537604,b=0       0.28022       0.28022       0.28022      0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022 
    coh29:k=0.550108,b=0       0.28022       0.28022       0.28022      0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022 
    coh30:k=0.562612,b=0       0.28022       0.28022       0.28022      0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022       0.28022 
    coh31:k=0.575116,b=0       0.32026       0.32026       0.32026      0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026 
    coh32:k=0.587619,b=0       0.32026       0.32026       0.32026      0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026 
    coh33:k=0.600123,b=0       0.32026       0.32026       0.32026      0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026       0.32026 
    coh34:k=0.612627,b=0       0.36029       0.36029       0.36029      0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029 
    coh35:k=0.625131,b=0       0.36029       0.36029       0.36029      0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029 
    coh36:k=0.637635,b=0       0.36029       0.36029       0.36029      0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029 
    coh37:k=0.650139,b=0       0.36029       0.36029       0.36029      0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029       0.36029 
    coh38:k=0.662643,b=0       0.40032       0.40032       0.40032      0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032 
    coh39:k=0.675146,b=0       0.40032       0.40032       0.40032      0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032 
    coh40:k=0.68765,b=0        0.40032       0.40032       0.40032      0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032       0.40032 
    coh41:k=0.700154,b=0       0.44035       0.44035       0.44035      0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035 
    coh42:k=0.712658,b=0       0.44035       0.44035       0.44035      0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035 
    coh43:k=0.725162,b=0       0.44035       0.44035       0.44035      0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035 
    coh44:k=0.737666,b=0       0.44035       0.44035       0.44035      0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035       0.44035 
    coh45:k=0.75017,b=0        0.48038       0.48038       0.48038      0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038 
    coh46:k=0.762673,b=0       0.48038       0.48038       0.48038      0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038 
    coh47:k=0.775177,b=0       0.48038       0.48038       0.48038      0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038 
    coh48:k=0.787681,b=0       0.48038       0.48038       0.48038      0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038       0.48038 
    coh49:k=0.800185,b=0       0.52042       0.52042       0.52042      0.52042       0.52042       0.52042       0.52042       0.52042       0.52042       0.52042       0.52042       0.52042 
    coh50:k=0.812689,b=0       0.52042       0.52042       0.52042      0.52042       0.52042       0.52042       0.52042       0.52042       0.52042       0.52042       0.52042       0.52042 
    coh4715:k=59.1432,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4716:k=59.1557,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4717:k=59.1682,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4718:k=59.1807,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4719:k=59.1932,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4720:k=59.2057,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4721:k=59.2182,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4722:k=59.2307,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4723:k=59.2432,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4724:k=59.2557,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4725:k=59.2682,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4726:k=59.2807,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4727:k=59.2932,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4728:k=59.3057,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4729:k=59.3182,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4730:k=59.3307,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4731:k=59.3432,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4732:k=59.3557,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4733:k=59.3682,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4734:k=59.3807,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4735:k=59.3932,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4736:k=59.4057,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4737:k=59.4183,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4738:k=59.4308,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4739:k=59.4433,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4740:k=59.4558,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4741:k=59.4683,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4742:k=59.4808,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4743:k=59.4933,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4744:k=59.5058,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4745:k=59.5183,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4746:k=59.5308,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4747:k=59.5433,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4748:k=59.5558,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4749:k=59.5683,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4750:k=59.5808,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4751:k=59.5933,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4752:k=59.6058,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4753:k=59.6183,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4754:k=59.6308,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4755:k=59.6433,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4756:k=59.6558,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4757:k=59.6683,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4758:k=59.6808,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4759:k=59.6933,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4760:k=59.7058,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4761:k=59.7183,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4762:k=59.7308,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4763:k=59.7434,b=0           50            50            50           50            50            50            50            50            50            50            50            50 
    coh4764:k=59.7559,b=0           50            50            50           50            50            50            50            50            50            50            50            50 

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_akz_set_default_param(it_param_set);

% Simulation Accuracy
param_map('it_w_n') = 1250;
param_map('it_ak_n') = param_map('it_w_n');
param_map('it_z_n') = 15;
param_map('fl_coh_interp_grid_gap') = 0.0125;
param_map('it_c_interp_grid_gap') = 10^-4;

% 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_iwkz_vf_vecsv(param_map, support_map);
Elapsed time is 63.845735 seconds.
Warning: Directory already exists.