Value Function iteration
This function solves the dynamic infinite horizon optimal asset choice problem in case of 3 occupational choices - worker, formal and informal entrepreneur
Back to full code page https://kritikhanna.github.io/BKS-modified/
Function Inputs : asset grid(1 by M), entrepreneurial productivity grid(1 by N), labor opportunity (1X2), and relevant model parameters required Calls function hhsave_util_2OC.m
Function Outputs : Value function (MXNX1), Policy Function (MXNX1), Optimal Location of asset choice (MXNX1)
Contents
function [mt_vf, mt_pf, mt_oploc,mt_coh, mt_o, mt_k, mt_l]=hhsave_VFI_3OC(varargin)
close all; bl_profile = false; % Switch off profile if running a tester/calling from another function bl_saveimg = true; if(bl_profile) profile off; profile on; end addpath(genpath('/Users/sidhantkhanna/Documents/GitHub/BKS modified/')); if ~isempty(varargin) [ar_a,ar_z,ar_n, ... fl_alpha,fl_theta,fl_delta,fl_kappa,fl_r,fl_w,fl_phi,fl_ahi,fl_zhi, ... fl_risk,it_zgridno, it_agridno,it_ngridno, mt_trans_z,mt_trans_n,fl_beta, fl_mu_z, fl_sig_z, ... fl_rho_z,fl_lambda_z, fl_tolvfi ... ] = varargin{:}; bl_print = false; bl_plot = false; else fl_ahi = 50; fl_alo = 0; fl_zhi = 2.2; it_agridno = 100; it_zgridno = 7; it_ngridno = 2; % ar_a = linspace(fl_alo ,fl_ahi,it_agridno); mp_grid_control = containers.Map('KeyType','char', 'ValueType','any'); mp_grid_control('grid_powerspace_power') = 1; % ar_a = ff_saveborr_grid(fl_alo, fl_ahi, it_agridno, mp_grid_control) ar_a = ff_saveborr_grid(fl_alo, fl_ahi, it_agridno, 'grid_powerspace', mp_grid_control); ar_a = ar_a'; fl_phi = 0; fl_risk = 1.5; fl_alpha = 0.4; fl_theta = 0.79-fl_alpha; fl_delta = 0.05; fl_kappa = 0; fl_mu_z = 0; % mean of AR(1) entrepeneurial productivity process fl_rho_z = 0.9; % persistence parameter of the AR(1) entrepreneurial productivity process fl_sig_z = 0.2; fl_lambda_z = 3; fl_beta = 0.92; % [ar_z, mt_trans_z] = mytauchen_z(fl_mu_z,fl_rho_z,fl_sig_z,it_zgridno,fl_lambda_z); [ar_z, mt_trans_z] = ffy_rouwenhorst(fl_rho_z,fl_sig_z,it_zgridno); ar_z = exp(ar_z); P1=mt_trans_z^1000; sd=P1(1,:); el=sd*ar_z; ar_z =ar_z/el; ar_z = ar_z'; ar_n = [0,1]; mt_trans_n = [0.1,0.9;0.9,0.1]; [fl_r,fl_w] = ... deal(0.02,1); fl_tolvfi = 10^-12; bl_print = true; bl_plot = true; end fl_R = fl_r + fl_delta;
Value function Iteration
mt_vf = zeros(it_agridno,it_zgridno,it_ngridno); % Value function mt_Tv = zeros(it_agridno,it_zgridno,it_ngridno); % Next guess of value function mt_oploc = zeros(it_agridno,it_zgridno,it_ngridno); % [mt_util,mt_coh, mt_k, mt_l, mt_o] = hhsave_util_3OC(ar_a,ar_z,ar_n, ... fl_alpha,fl_theta,fl_delta,fl_kappa,fl_r,fl_w,fl_phi,fl_ahi,fl_zhi, ... fl_risk,it_zgridno, it_agridno,it_ngridno); fl_crit = 1; % difference between current and next guess of value function mt_evf = zeros(it_agridno,it_zgridno,it_ngridno); while fl_crit > fl_tolvfi for i = 1:it_ngridno mt_evf(:,:,i) = mt_vf(:,:,i)*mt_trans_z'; end for j = 1:it_zgridno mt_evf1(:,j,:) = squeeze(mt_evf(:,j,:))*mt_trans_n'; end mt_ev1 = repmat(mt_evf1, [1 1 1 it_agridno]); mt_ev1 = permute(mt_ev1,[4 2 3 1]); [mt_Tvf,mt_oploc] = max((mt_util+fl_beta.*(mt_ev1)),[], 4); mt_diff = mt_Tvf-mt_vf; mt_diff = mt_diff(:); fl_crit = norm(mt_diff); mt_vf = mt_Tvf; end mt_pf = ar_a(mt_oploc); mt_con = mt_coh - mt_pf; mt_savefraccoh = mt_pf./mt_coh; if(bl_print) st_ar_a = strsplit(num2str(ar_a)); for i=1:it_agridno st_a(i) =strcat('a=', st_ar_a(i)); end st_ar_z = strsplit(num2str(ar_z)); disp('Table for value function at l=0'); tb_VF = array2table(mt_vf(:,:,1),'VariableNames',{'z1','z2','z3','z4','z5','z6','z7'}); tb_VF.Properties.RowNames = st_a; disp(tb_VF); disp('Table for value function at l=1'); tb_VF = array2table(mt_vf(:,:,2),'VariableNames',{'z1','z2','z3','z4','z5','z6','z7'}); tb_VF.Properties.RowNames = st_a; disp(tb_VF); disp('Table for Policy function at l=0'); tb_pf = array2table(mt_pf(:,:,1),'VariableNames',{'z1','z2','z3','z4','z5','z6','z7'}); tb_pf.Properties.RowNames = st_a; disp(tb_pf); disp('Table for Policy function at l=1'); tb_pf = array2table(mt_pf(:,:,2),'VariableNames',{'z1','z2','z3','z4','z5','z6','z7'}); tb_pf.Properties.RowNames = st_a; disp(tb_pf); disp('Table for Consumption at l=0 '); tb_con = array2table(mt_con(:,:,1),'VariableNames',{'z1','z2','z3','z4','z5','z6','z7'}); tb_con.Properties.RowNames = st_a; disp(tb_con); disp('Table for Consumption at l=1 '); tb_con = array2table(mt_con(:,:,2),'VariableNames',{'z1','z2','z3','z4','z5','z6','z7'}); tb_con.Properties.RowNames = st_a; disp(tb_con); disp('Table for Savings/COH at l=0'); tb_savefraccoh = array2table(mt_savefraccoh(:,:,1),'VariableNames',{'z1','z2','z3','z4','z5','z6','z7'}); tb_savefraccoh.Properties.RowNames = st_a; disp(tb_savefraccoh); disp('Table for Savings/COH at l=1'); tb_savefraccoh = array2table(mt_savefraccoh(:,:,2),'VariableNames',{'z1','z2','z3','z4','z5','z6','z7'}); tb_savefraccoh.Properties.RowNames = st_a; disp(tb_savefraccoh); end
Table for value function at l=0
z1 z2 z3 z4 z5 z6 z7
_________ ________ ________ ________ ________ _______ _______
a=0 -128.79 -124.85 -121.5 -118.99 -116.81 -114.58 -111.53
a=0.505051 -26.515 -19.609 -13.545 -9.7678 -6.5858 -3.3507 2.203
a=1.0101 -11.733 -10.247 -8.6593 -7.009 -4.6245 -1.2383 4.6005
a=1.51515 -8.4539 -7.8008 -6.8293 -5.4902 -3.3533 0.4626 5.9294
a=2.0202 -7.344 -6.8036 -5.8915 -4.5294 -2.3416 1.6251 6.9059
a=2.52525 -6.5588 -6.0275 -5.1182 -3.7194 -1.4539 2.5621 7.6243
a=3.0303 -5.8127 -5.3025 -4.424 -2.9735 -0.66393 3.3718 8.2128
a=3.53535 -5.1331 -4.6255 -3.7925 -2.2938 0.014983 4.0532 8.6784
a=4.0404 -4.4937 -3.9992 -3.1955 -1.6825 0.62389 4.6218 9.1061
a=4.54545 -3.9053 -3.4196 -2.6287 -1.1273 1.222 5.1157 9.4554
a=5.05051 -3.3542 -2.8797 -2.1007 -0.60814 1.7587 5.5505 9.775
a=5.55556 -2.8436 -2.3761 -1.6043 -0.12787 2.2442 5.9449 10.059
a=6.06061 -2.3657 -1.9059 -1.1415 0.31558 2.685 6.3026 10.305
a=6.56566 -1.9204 -1.4659 -0.7082 0.72456 3.0891 6.6365 10.543
a=7.07071 -1.5033 -1.0547 -0.30458 1.1027 3.4589 6.9407 10.755
a=7.57576 -1.1133 -0.6696 0.072452 1.4525 3.7986 7.2188 10.946
a=8.08081 -0.74758 -0.30936 0.42385 1.7766 4.1145 7.476 11.131
a=8.58586 -0.4047 0.028293 0.7521 2.0828 4.4117 7.7167 11.302
a=9.09091 -0.082906 0.34464 1.0597 2.3764 4.691 7.941 11.457
a=9.59596 0.21968 0.6429 1.3531 2.6637 4.9558 8.1515 11.604
a=10.101 0.51673 0.93551 1.6365 2.9363 5.2076 8.3518 11.746
a=10.6061 0.80528 1.2224 1.9051 3.1962 5.4459 8.5405 11.877
a=11.1111 1.0828 1.4953 2.1682 3.4467 5.6724 8.7194 12
a=11.6162 1.3462 1.7553 2.4232 3.6856 5.8888 8.8919 12.117
a=12.1212 1.6035 2.0082 2.6687 3.9135 6.0947 9.0555 12.23
a=12.6263 1.8535 2.2531 2.9037 4.1329 6.291 9.2108 12.337
a=13.1313 2.0937 2.4886 3.1316 4.3429 6.4786 9.359 12.438
a=13.6364 2.3235 2.7147 3.3513 4.5438 6.6579 9.5007 12.533
a=14.1414 2.547 2.9337 3.5629 4.7367 6.8295 9.6363 12.627
a=14.6465 2.7636 3.1453 3.7663 4.9218 6.9939 9.7663 12.717
a=15.1515 2.9721 3.3493 3.9628 5.0994 7.1523 9.8914 12.802
a=15.6566 3.1727 3.5459 4.152 5.2699 7.3052 10.012 12.883
a=16.1616 3.3672 3.7359 4.3344 5.4356 7.4524 10.128 12.96
a=16.6667 3.5555 3.9194 4.5105 5.5966 7.5944 10.24 13.037
a=17.1717 3.7372 4.0967 4.6807 5.7542 7.7319 10.348 13.11
a=17.6768 3.9124 4.2677 4.8449 5.9082 7.8653 10.453 13.18
a=18.1818 4.0822 4.4336 5.0057 6.0572 7.9944 10.555 13.247
a=18.6869 4.2484 4.5981 5.162 6.2018 8.1196 10.653 13.311
a=19.1919 4.4118 4.7581 5.3159 6.3422 8.241 10.749 13.374
a=19.697 4.571 4.9152 5.4665 6.4789 8.3589 10.842 13.435
a=20.202 4.7272 5.0674 5.6128 6.612 8.4733 10.932 13.494
a=20.7071 4.8784 5.2151 5.7545 6.7412 8.5845 11.02 13.55
a=21.2121 5.0256 5.3588 5.8929 6.8669 8.6924 11.105 13.605
a=21.7172 5.1699 5.5 6.0277 6.9893 8.7973 11.188 13.658
a=22.2222 5.3112 5.6378 6.1597 7.1086 8.8993 11.269 13.709
a=22.7273 5.4493 5.7728 6.2886 7.2248 8.9989 11.347 13.759
a=23.2323 5.5842 5.9041 6.4141 7.3377 9.0961 11.424 13.808
a=23.7374 5.7153 6.0317 6.5361 7.4477 9.1907 11.499 13.854
a=24.2424 5.8431 6.1562 6.6551 7.5559 9.2836 11.571 13.9
a=24.7475 5.9683 6.2781 6.7714 7.6623 9.3749 11.642 13.944
a=25.2525 6.0908 6.3972 6.8852 7.7677 9.4644 11.711 13.986
a=25.7576 6.2106 6.5137 6.9964 7.8705 9.5521 11.779 14.027
a=26.2626 6.3275 6.6272 7.1047 7.971 9.638 11.845 14.067
a=26.7677 6.4414 6.7379 7.2106 8.0691 9.7222 11.909 14.106
a=27.2727 6.5527 6.8467 7.315 8.1651 9.8047 11.972 14.144
a=27.7778 6.6617 6.9544 7.4184 8.2592 9.8854 12.034 14.181
a=28.2828 6.77 7.0599 7.52 8.3516 9.9644 12.094 14.216
a=28.7879 6.8762 7.1636 7.6195 8.4425 10.042 12.153 14.251
a=29.2929 6.9802 7.2657 7.7172 8.5315 10.117 12.21 14.285
a=29.798 7.0828 7.3657 7.8128 8.6185 10.191 12.267 14.318
a=30.303 7.1834 7.4635 7.9063 8.7036 10.264 12.322 14.35
a=30.8081 7.2817 7.559 7.9978 8.7872 10.335 12.376 14.381
a=31.3131 7.3781 7.653 8.0877 8.8692 10.405 12.429 14.411
a=31.8182 7.473 7.7458 8.1764 8.9497 10.473 12.48 14.439
a=32.3232 7.5666 7.8368 8.2635 9.0287 10.54 12.531 14.468
a=32.8283 7.6584 7.9262 8.349 9.1064 10.606 12.581 14.495
a=33.3333 7.7486 8.0142 8.4329 9.1837 10.671 12.63 14.522
a=33.8384 7.8374 8.1004 8.5153 9.2602 10.734 12.677 14.547
a=34.3434 7.9244 8.1849 8.5961 9.3351 10.796 12.724 14.572
a=34.8485 8.0096 8.2677 8.6752 9.4085 10.857 12.77 14.597
a=35.3535 8.0933 8.349 8.753 9.4805 10.917 12.815 14.62
a=35.8586 8.1757 8.4293 8.8296 9.5514 10.976 12.86 14.643
a=36.3636 8.2568 8.5081 8.9051 9.6211 11.034 12.903 14.666
a=36.8687 8.3365 8.5861 8.9803 9.6897 11.091 12.945 14.687
a=37.3737 8.4149 8.6633 9.0546 9.757 11.147 12.987 14.708
a=37.8788 8.4927 8.7393 9.1275 9.8236 11.203 13.028 14.729
a=38.3838 8.5694 8.8139 9.1994 9.8892 11.257 13.068 14.748
a=38.8889 8.6449 8.8878 9.2702 9.9536 11.311 13.107 14.767
a=39.3939 8.7189 8.9605 9.3396 10.017 11.364 13.146 14.785
a=39.899 8.7923 9.0318 9.4079 10.079 11.417 13.183 14.803
a=40.404 8.8645 9.1019 9.475 10.14 11.468 13.22 14.82
a=40.9091 8.9357 9.1712 9.5411 10.201 11.519 13.257 14.836
a=41.4141 9.0057 9.2393 9.6062 10.26 11.569 13.292 14.852
a=41.9192 9.0749 9.3066 9.6707 10.32 11.618 13.326 14.867
a=42.4242 9.143 9.3732 9.7344 10.378 11.667 13.36 14.882
a=42.9293 9.2103 9.4387 9.7971 10.436 11.714 13.393 14.896
a=43.4343 9.2766 9.5032 9.8588 10.493 11.761 13.426 14.91
a=43.9394 9.3421 9.567 9.9197 10.549 11.808 13.457 14.923
a=44.4444 9.4064 9.6297 9.9796 10.605 11.853 13.488 14.935
a=44.9495 9.4699 9.6913 10.039 10.659 11.898 13.518 14.948
a=45.4545 9.5324 9.7521 10.097 10.713 11.942 13.547 14.959
a=45.9596 9.5941 9.812 10.155 10.766 11.986 13.575 14.97
a=46.4646 9.6549 9.8715 10.212 10.818 12.029 13.603 14.98
a=46.9697 9.7149 9.9304 10.268 10.87 12.071 13.629 14.99
a=47.4747 9.7742 9.9886 10.324 10.921 12.112 13.655 14.999
a=47.9798 9.8332 10.046 10.379 10.971 12.154 13.68 15.007
a=48.4848 9.8915 10.103 10.434 11.021 12.194 13.703 15.015
a=48.9899 9.949 10.159 10.488 11.07 12.234 13.725 15.023
a=49.4949 10.006 10.214 10.541 11.118 12.274 13.746 15.03
a=50 10.062 10.269 10.593 11.166 12.313 13.764 15.037
Table for value function at l=1
z1 z2 z3 z4 z5 z6 z7
_________ ________ ________ ________ ________ ________ ________
a=0 -22.08 -17.015 -12.671 -9.6475 -7.0415 -4.3433 -0.25402
a=0.505051 -10.9 -9.7191 -8.4489 -7.0461 -5.125 -2.4328 2.1385
a=1.0101 -8.268 -7.6712 -6.8278 -5.6655 -3.8893 -0.96185 4.5978
a=1.51515 -7.2748 -6.7517 -5.9429 -4.7651 -2.9443 0.4311 5.9276
a=2.0202 -6.476 -5.9529 -5.1441 -3.9663 -2.1188 1.6079 6.9053
a=2.52525 -5.7519 -5.234 -4.4212 -3.2122 -1.3367 2.548 7.6245
a=3.0303 -5.0677 -4.567 -3.7664 -2.5223 -0.61235 3.3707 8.2129
a=3.53535 -4.44 -3.9478 -3.169 -1.8979 0.016487 4.0522 8.6785
a=4.0404 -3.8514 -3.3721 -2.6101 -1.3324 0.5945 4.621 9.1062
a=4.54545 -3.3078 -2.8363 -2.0837 -0.80973 1.1953 5.1149 9.4554
a=5.05051 -2.7989 -2.3366 -1.5927 -0.3235 1.7346 5.5498 9.775
a=5.55556 -2.3258 -1.8694 -1.1325 0.12649 2.2226 5.9443 10.059
a=6.06061 -1.8829 -1.4329 -0.70364 0.5422 2.6657 6.3021 10.305
a=6.56566 -1.4692 -1.0242 -0.30276 0.92669 3.0717 6.6361 10.543
a=7.07071 -1.0816 -0.64205 0.070786 1.2827 3.4434 6.9403 10.755
a=7.57576 -0.71851 -0.28397 0.41969 1.6127 3.7847 7.2184 10.946
a=8.08081 -0.37793 0.051167 0.7453 1.9272 4.1015 7.4756 11.131
a=8.58586 -0.058169 0.36572 1.0564 2.231 4.399 7.7163 11.302
a=9.09091 0.24959 0.67348 1.3588 2.5182 4.6787 7.9407 11.457
a=9.59596 0.54745 0.96651 1.6442 2.7921 4.9446 8.1512 11.604
a=10.101 0.82871 1.244 1.917 3.0578 5.1969 8.3515 11.746
a=10.6061 1.1037 1.5145 2.18 3.3111 5.4357 8.5402 11.877
a=11.1111 1.3714 1.7778 2.4308 3.5525 5.663 8.7191 12
a=11.6162 1.6287 2.0303 2.6746 3.7852 5.8799 8.8917 12.117
a=12.1212 1.8739 2.272 2.9103 4.0078 6.0863 9.0553 12.23
a=12.6263 2.1125 2.5063 3.1372 4.2207 6.2831 9.2106 12.337
a=13.1313 2.3443 2.7329 3.355 4.4253 6.4712 9.3588 12.438
a=13.6364 2.5673 2.9513 3.5656 4.6215 6.6509 9.5005 12.533
a=14.1414 2.7812 3.1613 3.7686 4.8097 6.8228 9.6361 12.627
a=14.6465 2.9888 3.3645 3.9641 4.9905 6.9876 9.7661 12.717
a=15.1515 3.1899 3.5607 4.1523 5.167 7.1462 9.8912 12.802
a=15.6566 3.3838 3.7501 4.3339 5.3394 7.2989 10.012 12.883
a=16.1616 3.5706 3.9327 4.509 5.5053 7.4462 10.128 12.96
a=16.6667 3.7515 4.1092 4.6816 5.665 7.5883 10.24 13.037
a=17.1717 3.9267 4.2805 4.8494 5.8198 7.7261 10.348 13.11
a=17.6768 4.0969 4.45 5.0119 5.97 7.8597 10.453 13.18
a=18.1818 4.2649 4.6139 5.1693 6.1166 7.9891 10.555 13.247
a=18.6869 4.4272 4.7723 5.3216 6.2594 8.1144 10.653 13.311
a=19.1919 4.5849 4.9264 5.4703 6.398 8.2359 10.749 13.374
a=19.697 4.7394 5.0781 5.615 6.5326 8.354 10.842 13.435
a=20.202 4.8908 5.2261 5.757 6.6635 8.4687 10.932 13.494
a=20.7071 5.0388 5.3713 5.8959 6.7911 8.5799 11.02 13.55
a=21.2121 5.1836 5.5123 6.0309 6.9153 8.688 11.105 13.605
a=21.7172 5.324 5.6493 6.162 7.036 8.793 11.188 13.658
a=22.2222 5.4608 5.7827 6.2899 7.1534 8.8952 11.269 13.709
a=22.7273 5.5948 5.9134 6.4146 7.2694 8.9949 11.347 13.759
a=23.2323 5.726 6.0411 6.5365 7.3832 9.092 11.424 13.808
a=23.7374 5.8543 6.1661 6.6557 7.4941 9.1865 11.499 13.854
a=24.2424 5.9795 6.2879 6.7717 7.6023 9.2795 11.571 13.9
a=24.7475 6.1014 6.4064 6.8852 7.7077 9.3708 11.642 13.944
a=25.2525 6.2203 6.5219 6.9975 7.8104 9.4604 11.711 13.986
a=25.7576 6.3367 6.6351 7.1072 7.9113 9.5483 11.779 14.027
a=26.2626 6.4506 6.7467 7.2146 8.0103 9.6344 11.845 14.067
a=26.7677 6.562 6.8567 7.3199 8.108 9.7187 11.909 14.106
a=27.2727 6.6723 6.9644 7.4228 8.2035 9.8012 11.972 14.144
a=27.7778 6.7805 7.0695 7.5233 8.2969 9.8819 12.034 14.181
a=28.2828 6.8861 7.1722 7.6215 8.3882 9.961 12.094 14.216
a=28.7879 6.9894 7.273 7.7181 8.4776 10.038 12.153 14.251
a=29.2929 7.091 7.3727 7.8136 8.5654 10.114 12.21 14.285
a=29.798 7.1914 7.4705 7.9073 8.6516 10.188 12.266 14.318
a=30.303 7.2899 7.5666 7.9993 8.7363 10.261 12.322 14.35
a=30.8081 7.3866 7.6611 8.0896 8.8195 10.332 12.376 14.381
a=31.3131 7.4819 7.7538 8.1781 8.9021 10.402 12.428 14.411
a=31.8182 7.5752 7.8444 8.2646 8.983 10.47 12.48 14.44
a=32.3232 7.6665 7.9332 8.3494 9.0621 10.537 12.531 14.468
a=32.8283 7.7561 8.0204 8.4328 9.1398 10.603 12.581 14.495
a=33.3333 7.8444 8.1065 8.5149 9.2161 10.668 12.63 14.522
a=33.8384 7.9313 8.1909 8.5959 9.2911 10.731 12.677 14.547
a=34.3434 8.0166 8.2739 8.6762 9.3647 10.793 12.724 14.572
a=34.8485 8.1005 8.3555 8.755 9.437 10.855 12.77 14.597
a=35.3535 8.183 8.4356 8.8326 9.5087 10.915 12.815 14.62
a=35.8586 8.2639 8.5153 8.9089 9.5794 10.974 12.86 14.643
a=36.3636 8.3432 8.5935 8.9837 9.6488 11.032 12.903 14.666
a=36.8687 8.4221 8.6702 9.0571 9.7168 11.089 12.945 14.687
a=37.3737 8.4997 8.7456 9.1292 9.7837 11.145 12.987 14.708
a=37.8788 8.5761 8.8199 9.2003 9.8495 11.2 13.028 14.729
a=38.3838 8.6513 8.8931 9.2703 9.9142 11.255 13.068 14.748
a=38.8889 8.7254 8.9655 9.3398 9.978 11.309 13.107 14.767
a=39.3939 8.7984 9.037 9.4083 10.041 11.362 13.146 14.785
a=39.899 8.8707 9.1074 9.4757 10.104 11.414 13.183 14.803
a=40.404 8.9418 9.1767 9.5421 10.165 11.466 13.22 14.82
a=40.9091 9.012 9.2452 9.6076 10.226 11.517 13.257 14.836
a=41.4141 9.0809 9.3125 9.6719 10.286 11.567 13.292 14.852
a=41.9192 9.1491 9.3787 9.7351 10.345 11.616 13.326 14.867
a=42.4242 9.2161 9.4439 9.7974 10.403 11.664 13.36 14.882
a=42.9293 9.2823 9.5081 9.859 10.46 11.712 13.393 14.896
a=43.4343 9.3474 9.5714 9.9201 10.516 11.759 13.426 14.91
a=43.9394 9.4116 9.634 9.9803 10.571 11.806 13.457 14.923
a=44.4444 9.475 9.6958 10.04 10.626 11.851 13.488 14.936
a=44.9495 9.5376 9.7566 10.099 10.68 11.896 13.518 14.948
a=45.4545 9.5993 9.8171 10.157 10.734 11.94 13.547 14.959
a=45.9596 9.6602 9.8769 10.214 10.787 11.984 13.575 14.97
a=46.4646 9.7202 9.9359 10.27 10.839 12.027 13.603 14.98
a=46.9697 9.7799 9.9939 10.326 10.89 12.069 13.629 14.99
a=47.4747 9.8388 10.051 10.381 10.941 12.111 13.655 14.999
a=47.9798 9.8968 10.107 10.434 10.991 12.152 13.679 15.007
a=48.4848 9.9541 10.163 10.488 11.04 12.192 13.703 15.015
a=48.9899 10.011 10.218 10.541 11.089 12.233 13.725 15.023
a=49.4949 10.066 10.273 10.593 11.138 12.272 13.746 15.03
a=50 10.122 10.327 10.645 11.186 12.311 13.764 15.037
Table for Policy function at l=0
z1 z2 z3 z4 z5 z6 z7
_______ _______ _______ _______ _______ _______ ______
a=0 0 0 0 0 0 0 0
a=0.505051 0.50505 0.50505 0.50505 0.50505 0.50505 0.50505 1.0101
a=1.0101 0.50505 0.50505 0.50505 0.50505 1.0101 1.5152 2.0202
a=1.51515 1.0101 1.0101 1.0101 1.0101 1.5152 2.0202 3.0303
a=2.0202 1.5152 1.5152 1.5152 1.5152 2.0202 2.5253 3.5354
a=2.52525 2.0202 2.0202 2.0202 2.0202 2.5253 3.0303 4.5455
a=3.0303 2.0202 2.0202 2.5253 2.5253 3.0303 4.0404 5.0505
a=3.53535 2.5253 2.5253 3.0303 3.0303 3.5354 4.5455 6.0606
a=4.0404 3.0303 3.0303 3.0303 3.5354 4.0404 5.0505 6.5657
a=4.54545 3.5354 3.5354 3.5354 4.0404 4.5455 5.5556 7.0707
a=5.05051 4.0404 4.0404 4.0404 4.5455 5.0505 6.0606 8.0808
a=5.55556 4.5455 4.5455 4.5455 5.0505 5.5556 6.5657 8.5859
a=6.06061 5.0505 5.0505 5.0505 5.5556 6.0606 7.5758 9.0909
a=6.56566 5.5556 5.5556 5.5556 6.0606 6.5657 8.0808 10.101
a=7.07071 6.0606 6.0606 6.0606 6.5657 7.0707 8.5859 10.606
a=7.57576 6.5657 6.5657 6.5657 7.0707 7.5758 9.0909 11.111
a=8.08081 7.0707 7.0707 7.0707 7.5758 8.0808 9.596 12.121
a=8.58586 7.5758 7.5758 7.5758 8.0808 8.5859 10.101 12.626
a=9.09091 8.0808 8.0808 8.0808 8.0808 9.0909 10.606 13.131
a=9.59596 8.5859 8.5859 8.5859 8.5859 9.596 11.111 14.141
a=10.101 8.5859 8.5859 9.0909 9.0909 10.101 11.616 14.646
a=10.6061 9.0909 9.0909 9.596 9.596 10.606 12.121 15.152
a=11.1111 9.596 9.596 9.596 10.101 11.111 13.131 15.657
a=11.6162 10.101 10.101 10.101 10.606 11.616 13.636 16.667
a=12.1212 10.606 10.606 10.606 11.111 12.121 14.141 17.172
a=12.6263 11.111 11.111 11.111 11.616 12.626 14.646 17.677
a=13.1313 11.616 11.616 11.616 12.121 13.131 15.152 18.182
a=13.6364 12.121 12.121 12.121 12.626 13.636 15.657 18.687
a=14.1414 12.626 12.626 12.626 13.131 14.141 16.162 19.697
a=14.6465 13.131 13.131 13.131 13.636 14.646 16.667 20.202
a=15.1515 13.636 13.636 13.636 14.141 15.152 17.172 20.707
a=15.6566 14.141 14.141 14.141 14.646 15.657 17.677 21.212
a=16.1616 14.646 14.646 14.646 15.152 16.162 18.182 22.222
a=16.6667 15.152 15.152 15.152 15.657 16.667 18.687 22.727
a=17.1717 15.657 15.657 15.657 15.657 17.172 19.192 23.232
a=17.6768 16.162 16.162 16.162 16.162 17.677 19.697 23.737
a=18.1818 16.667 16.667 16.667 16.667 18.182 20.202 24.242
a=18.6869 16.667 16.667 17.172 17.172 18.687 20.707 24.747
a=19.1919 17.172 17.172 17.172 17.677 19.192 21.717 25.758
a=19.697 17.677 17.677 17.677 18.182 19.697 22.222 26.263
a=20.202 18.182 18.182 18.182 18.687 20.202 22.727 26.768
a=20.7071 18.687 18.687 18.687 19.192 20.707 23.232 27.273
a=21.2121 19.192 19.192 19.192 19.697 21.212 23.737 27.778
a=21.7172 19.697 19.697 19.697 20.202 21.717 24.242 28.283
a=22.2222 20.202 20.202 20.202 20.707 22.222 24.747 28.788
a=22.7273 20.707 20.707 20.707 21.212 22.727 25.253 29.798
a=23.2323 21.212 21.212 21.212 21.717 23.232 25.758 30.303
a=23.7374 21.717 21.717 21.717 22.222 23.737 26.263 30.808
a=24.2424 22.222 22.222 22.222 22.727 23.737 26.768 31.313
a=24.7475 22.727 22.727 22.727 22.727 24.242 27.273 31.818
a=25.2525 23.232 23.232 23.232 23.232 24.747 27.778 32.323
a=25.7576 23.737 23.737 23.737 23.737 25.253 28.283 32.828
a=26.2626 24.242 24.242 24.242 24.242 25.758 28.788 33.333
a=26.7677 24.747 24.747 24.747 24.747 26.263 29.293 33.838
a=27.2727 25.253 24.747 25.253 25.253 26.768 29.798 34.343
a=27.7778 25.253 25.253 25.253 25.758 27.273 30.303 35.354
a=28.2828 25.758 25.758 25.758 26.263 27.778 30.808 35.859
a=28.7879 26.263 26.263 26.263 26.768 28.283 31.313 36.364
a=29.2929 26.768 26.768 26.768 27.273 28.788 31.818 36.869
a=29.798 27.273 27.273 27.273 27.778 29.293 32.323 37.374
a=30.303 27.778 27.778 27.778 28.283 29.798 32.828 37.879
a=30.8081 28.283 28.283 28.283 28.788 30.303 33.333 38.384
a=31.3131 28.788 28.788 28.788 29.293 30.808 33.838 38.889
a=31.8182 29.293 29.293 29.293 29.798 31.313 34.343 39.394
a=32.3232 29.798 29.798 29.798 30.303 31.818 34.848 39.899
a=32.8283 30.303 30.303 30.303 30.808 32.323 35.354 40.404
a=33.3333 30.808 30.808 30.808 30.808 32.828 35.859 40.909
a=33.8384 31.313 31.313 31.313 31.313 33.333 36.364 41.414
a=34.3434 31.818 31.818 31.818 31.818 33.838 36.869 41.919
a=34.8485 32.323 32.323 32.323 32.323 34.343 37.374 42.424
a=35.3535 32.828 32.828 32.828 32.828 34.848 37.879 42.929
a=35.8586 33.333 33.333 33.333 33.333 35.354 38.384 43.434
a=36.3636 33.838 33.838 33.333 33.838 35.859 38.889 43.939
a=36.8687 34.343 33.838 33.838 34.343 36.364 39.394 43.939
a=37.3737 34.848 34.343 34.343 34.848 36.364 39.899 44.444
a=37.8788 34.848 34.848 34.848 35.354 36.869 40.404 44.949
a=38.3838 35.354 35.354 35.354 35.859 37.374 40.909 45.455
a=38.8889 35.859 35.859 35.859 36.364 37.879 41.414 45.96
a=39.3939 36.364 36.364 36.364 36.869 38.384 41.919 46.465
a=39.899 36.869 36.869 36.869 37.374 38.889 42.424 46.465
a=40.404 37.374 37.374 37.374 37.879 39.394 42.929 46.97
a=40.9091 37.879 37.879 37.879 38.384 39.899 43.434 47.475
a=41.4141 38.384 38.384 38.384 38.384 40.404 43.939 47.98
a=41.9192 38.889 38.889 38.889 38.889 40.909 44.444 47.98
a=42.4242 39.394 39.394 39.394 39.394 41.414 44.949 48.485
a=42.9293 39.899 39.899 39.899 39.899 41.919 45.455 48.99
a=43.4343 40.404 40.404 40.404 40.404 42.424 45.96 48.99
a=43.9394 40.909 40.909 40.909 40.909 42.929 46.465 49.495
a=44.4444 41.414 41.414 41.414 41.414 43.434 46.97 50
a=44.9495 41.919 41.919 41.919 41.919 43.939 46.97 50
a=45.4545 42.424 42.424 41.919 42.424 44.444 47.475 50
a=45.9596 42.929 42.929 42.424 42.929 44.949 47.98 50
a=46.4646 43.434 42.929 42.929 43.434 45.455 48.485 50
a=46.9697 43.939 43.434 43.434 43.939 45.96 48.99 50
a=47.4747 43.939 43.939 43.939 44.444 45.96 48.99 50
a=47.9798 44.444 44.444 44.444 44.949 46.465 49.495 50
a=48.4848 44.949 44.949 44.949 45.455 46.97 50 50
a=48.9899 45.455 45.455 45.455 45.96 47.475 50 50
a=49.4949 45.96 45.96 45.96 46.465 47.98 50 50
a=50 46.465 46.465 46.465 46.465 48.485 50 50
Table for Policy function at l=1
z1 z2 z3 z4 z5 z6 z7
_______ _______ _______ _______ _______ _______ _______
a=0 0.50505 0.50505 0.50505 0.50505 0.50505 0.50505 0.50505
a=0.505051 1.0101 1.0101 1.0101 1.0101 1.0101 1.0101 1.0101
a=1.0101 1.5152 1.5152 1.5152 1.5152 1.5152 1.5152 2.0202
a=1.51515 2.0202 2.0202 2.0202 2.0202 2.0202 2.0202 3.0303
a=2.0202 2.0202 2.0202 2.0202 2.0202 2.5253 2.5253 3.5354
a=2.52525 2.5253 2.5253 2.5253 2.5253 2.5253 3.0303 4.5455
a=3.0303 3.0303 3.0303 3.0303 3.0303 3.0303 4.0404 5.0505
a=3.53535 3.5354 3.5354 3.5354 3.5354 3.5354 4.5455 6.0606
a=4.0404 4.0404 4.0404 4.0404 4.0404 4.0404 5.0505 6.5657
a=4.54545 4.5455 4.5455 4.5455 4.5455 4.5455 5.5556 7.0707
a=5.05051 5.0505 5.0505 5.0505 5.0505 5.0505 6.0606 8.0808
a=5.55556 5.5556 5.5556 5.5556 5.5556 5.5556 6.5657 8.5859
a=6.06061 6.0606 6.0606 6.0606 6.0606 6.0606 7.5758 9.0909
a=6.56566 6.5657 6.5657 6.5657 6.5657 6.5657 8.0808 10.101
a=7.07071 7.0707 7.0707 7.0707 7.0707 7.0707 8.5859 10.606
a=7.57576 7.5758 7.5758 7.5758 7.5758 7.5758 9.0909 11.111
a=8.08081 8.0808 8.0808 8.0808 7.5758 8.0808 9.596 12.121
a=8.58586 8.5859 8.5859 8.0808 8.0808 8.5859 10.101 12.626
a=9.09091 8.5859 8.5859 8.5859 8.5859 9.0909 10.606 13.131
a=9.59596 9.0909 9.0909 9.0909 9.0909 9.596 11.111 14.141
a=10.101 9.596 9.596 9.596 9.596 10.101 11.616 14.646
a=10.6061 10.101 10.101 10.101 10.101 10.606 12.121 15.152
a=11.1111 10.606 10.606 10.606 10.606 11.111 13.131 15.657
a=11.6162 11.111 11.111 11.111 11.111 11.616 13.636 16.667
a=12.1212 11.616 11.616 11.616 11.616 12.121 14.141 17.172
a=12.6263 12.121 12.121 12.121 12.121 12.626 14.646 17.677
a=13.1313 12.626 12.626 12.626 12.626 13.131 15.152 18.182
a=13.6364 13.131 13.131 13.131 13.131 13.636 15.657 18.687
a=14.1414 13.636 13.636 13.636 13.636 14.141 16.162 19.697
a=14.6465 14.141 14.141 14.141 14.141 14.646 16.667 20.202
a=15.1515 14.646 14.646 14.646 14.141 15.152 17.172 20.707
a=15.6566 15.152 15.152 15.152 14.646 15.657 17.677 21.212
a=16.1616 15.657 15.657 15.657 15.152 16.162 18.182 22.222
a=16.6667 16.162 16.162 15.657 15.657 16.667 18.687 22.727
a=17.1717 16.667 16.162 16.162 16.162 17.172 19.192 23.232
a=17.6768 16.667 16.667 16.667 16.667 17.677 19.697 23.737
a=18.1818 17.172 17.172 17.172 17.172 18.182 20.202 24.242
a=18.6869 17.677 17.677 17.677 17.677 18.687 20.707 24.747
a=19.1919 18.182 18.182 18.182 18.182 19.192 21.717 25.758
a=19.697 18.687 18.687 18.687 18.687 19.697 22.222 26.263
a=20.202 19.192 19.192 19.192 19.192 20.202 22.727 26.768
a=20.7071 19.697 19.697 19.697 19.697 20.707 23.232 27.273
a=21.2121 20.202 20.202 20.202 20.202 21.212 23.737 27.778
a=21.7172 20.707 20.707 20.707 20.707 21.717 24.242 28.283
a=22.2222 21.212 21.212 21.212 21.212 22.222 24.747 28.788
a=22.7273 21.717 21.717 21.717 21.212 22.727 25.253 29.798
a=23.2323 22.222 22.222 22.222 21.717 23.232 25.758 30.303
a=23.7374 22.727 22.727 22.727 22.222 23.737 26.263 30.808
a=24.2424 23.232 23.232 23.232 22.727 23.737 26.768 31.313
a=24.7475 23.737 23.737 23.232 23.232 24.242 27.273 31.818
a=25.2525 24.242 24.242 23.737 23.737 24.747 27.778 32.323
a=25.7576 24.747 24.747 24.242 24.242 25.253 28.283 32.828
a=26.2626 25.253 24.747 24.747 24.747 25.758 28.788 33.333
a=26.7677 25.758 25.253 25.253 25.253 26.263 29.293 33.838
a=27.2727 25.758 25.758 25.758 25.758 26.768 29.798 34.343
a=27.7778 26.263 26.263 26.263 26.263 27.273 30.303 35.354
a=28.2828 26.768 26.768 26.768 26.768 27.778 30.808 35.859
a=28.7879 27.273 27.273 27.273 27.273 28.283 31.313 36.364
a=29.2929 27.778 27.778 27.778 27.778 28.788 31.818 36.869
a=29.798 28.283 28.283 28.283 28.283 29.293 32.323 37.374
a=30.303 28.788 28.788 28.788 28.788 29.798 32.828 37.879
a=30.8081 29.293 29.293 29.293 28.788 30.303 33.333 38.384
a=31.3131 29.798 29.798 29.798 29.293 30.808 33.838 38.889
a=31.8182 30.303 30.303 30.303 29.798 31.313 34.343 39.394
a=32.3232 30.808 30.808 30.808 30.303 31.818 34.848 39.899
a=32.8283 31.313 31.313 31.313 30.808 32.323 35.354 40.404
a=33.3333 31.818 31.818 31.818 31.313 32.828 35.859 40.909
a=33.8384 32.323 32.323 31.818 31.818 33.333 36.364 41.414
a=34.3434 32.828 32.828 32.323 32.323 33.838 36.869 41.919
a=34.8485 33.333 33.333 32.828 32.828 34.343 37.374 42.424
a=35.3535 33.838 33.838 33.333 33.333 34.848 37.879 42.929
a=35.8586 34.343 33.838 33.838 33.838 35.354 38.384 43.434
a=36.3636 34.848 34.343 34.343 34.343 35.859 38.889 43.939
a=36.8687 34.848 34.848 34.848 34.848 36.364 39.394 43.939
a=37.3737 35.354 35.354 35.354 35.354 36.364 39.899 44.444
a=37.8788 35.859 35.859 35.859 35.859 36.869 40.404 44.949
a=38.3838 36.364 36.364 36.364 36.364 37.374 40.909 45.455
a=38.8889 36.869 36.869 36.869 36.869 37.879 41.414 45.96
a=39.3939 37.374 37.374 37.374 36.869 38.384 41.919 46.465
a=39.899 37.879 37.879 37.879 37.374 38.889 42.424 46.465
a=40.404 38.384 38.384 38.384 37.879 39.394 42.929 46.97
a=40.9091 38.889 38.889 38.889 38.384 39.899 43.434 47.475
a=41.4141 39.394 39.394 39.394 38.889 40.404 43.939 47.98
a=41.9192 39.899 39.899 39.899 39.394 40.909 44.444 47.98
a=42.4242 40.404 40.404 40.404 39.899 41.414 44.949 48.485
a=42.9293 40.909 40.909 40.404 40.404 41.919 45.455 48.99
a=43.4343 41.414 41.414 40.909 40.909 42.424 45.96 48.99
a=43.9394 41.919 41.919 41.414 41.414 42.929 46.465 49.495
a=44.4444 42.424 42.424 41.919 41.919 43.434 46.97 50
a=44.9495 42.929 42.929 42.424 42.424 43.939 46.97 50
a=45.4545 43.434 42.929 42.929 42.929 44.444 47.475 50
a=45.9596 43.939 43.434 43.434 43.434 44.949 47.98 50
a=46.4646 43.939 43.939 43.939 43.939 45.455 48.485 50
a=46.9697 44.444 44.444 44.444 44.444 45.96 48.99 50
a=47.4747 44.949 44.949 44.949 44.949 45.96 48.99 50
a=47.9798 45.455 45.455 45.455 44.949 46.465 49.495 50
a=48.4848 45.96 45.96 45.96 45.455 46.97 50 50
a=48.9899 46.465 46.465 46.465 45.96 47.475 50 50
a=49.4949 46.97 46.97 46.97 46.465 47.98 50 50
a=50 47.475 47.475 47.475 46.97 48.485 50 50
Table for Consumption at l=0
z1 z2 z3 z4 z5 z6 z7
________ ________ ________ _______ _______ _______ _______
a=0 0 0 0 0 0 0 0
a=0.505051 0.013011 0.027424 0.072047 0.15457 0.30707 0.58891 0.60473
a=1.0101 0.52816 0.54258 0.60783 0.73784 0.47304 0.41201 0.72755
a=1.51515 0.53826 0.55268 0.62927 0.79887 0.60725 0.68146 0.74187
a=2.0202 0.54836 0.56278 0.64553 0.85034 0.7238 0.91826 1.2009
a=2.52525 0.55847 0.57288 0.65833 0.89541 0.82851 1.1332 1.1145
a=3.0303 1.0736 1.088 0.66879 0.93577 0.92451 0.82698 1.5034
a=3.53535 1.0837 1.0981 0.67889 0.97245 1.0137 1.0133 1.364
a=4.0404 1.0938 1.1082 1.194 1.0061 1.0974 1.1893 1.7108
a=4.54545 1.1039 1.1183 1.2041 1.0373 1.1765 1.3569 2.0419
a=5.05051 1.114 1.1284 1.2142 1.0664 1.2516 1.5172 1.8546
a=5.55556 1.1241 1.1385 1.2243 1.0937 1.3234 1.6713 2.1607
a=6.06061 1.1342 1.1486 1.2344 1.1193 1.3922 1.3147 2.4565
a=6.56566 1.1443 1.1587 1.2445 1.1435 1.4583 1.4583 2.2382
a=7.07071 1.1544 1.1688 1.2546 1.1664 1.522 1.5974 2.5168
a=7.57576 1.1645 1.1789 1.2647 1.1881 1.5835 1.7326 2.788
a=8.08081 1.1746 1.189 1.2749 1.2087 1.6431 1.8641 2.5474
a=8.58586 1.1847 1.1991 1.285 1.2284 1.7009 1.9923 2.8058
a=9.09091 1.1948 1.2092 1.2951 1.7522 1.757 2.1174 3.0584
a=9.59596 1.2049 1.2193 1.3052 1.7701 1.8115 2.2396 2.8007
a=10.101 1.7201 1.7345 1.3153 1.7873 1.8647 2.3592 3.0431
a=10.6061 1.7302 1.7446 1.3254 1.8037 1.9164 2.4764 3.281
a=11.1111 1.7403 1.7547 1.8405 1.8195 1.967 2.0861 3.5146
a=11.6162 1.7504 1.7648 1.8506 1.8346 2.0163 2.1988 3.2391
a=12.1212 1.7605 1.7749 1.8607 1.8491 2.0646 2.3093 3.4649
a=12.6263 1.7706 1.785 1.8708 1.8631 2.1118 2.418 3.6871
a=13.1313 1.7807 1.7951 1.8809 1.8765 2.158 2.5248 3.9059
a=13.6364 1.7908 1.8052 1.891 1.8894 2.2033 2.6299 4.1216
a=14.1414 1.8009 1.8153 1.9011 1.9018 2.2476 2.7333 3.8291
a=14.6465 1.811 1.8254 1.9112 1.9138 2.2912 2.8352 4.0387
a=15.1515 1.8211 1.8355 1.9213 1.9253 2.3339 2.9355 4.2456
a=15.6566 1.8312 1.8456 1.9314 1.9364 2.3758 3.0345 4.4498
a=16.1616 1.8413 1.8557 1.9415 1.9471 2.417 3.132 4.1465
a=16.6667 1.8514 1.8658 1.9516 1.9574 2.4575 3.2282 4.3458
a=17.1717 1.8615 1.8759 1.9617 2.4726 2.4973 3.3232 4.5427
a=17.6768 1.8716 1.886 1.9718 2.4827 2.5365 3.417 4.7374
a=18.1818 1.8817 1.8961 1.9819 2.4928 2.575 3.5096 4.93
a=18.6869 2.3968 2.4113 1.992 2.5029 2.6129 3.6011 5.1204
a=19.1919 2.407 2.4214 2.5072 2.513 2.6502 3.1864 4.8039
a=19.697 2.4171 2.4315 2.5173 2.5231 2.687 3.2758 4.9904
a=20.202 2.4272 2.4416 2.5274 2.5332 2.7232 3.3641 5.1751
a=20.7071 2.4373 2.4517 2.5375 2.5433 2.7589 3.4515 5.358
a=21.2121 2.4474 2.4618 2.5476 2.5534 2.7941 3.5379 5.5391
a=21.7172 2.4575 2.4719 2.5577 2.5635 2.8287 3.6234 5.7186
a=22.2222 2.4676 2.482 2.5678 2.5736 2.8629 3.7081 5.8964
a=22.7273 2.4777 2.4921 2.5779 2.5837 2.8967 3.7918 5.5676
a=23.2323 2.4878 2.5022 2.588 2.5938 2.93 3.8748 5.7423
a=23.7374 2.4979 2.5123 2.5981 2.6039 2.9628 3.9569 5.9155
a=24.2424 2.508 2.5224 2.6082 2.614 3.5003 4.0383 6.0873
a=24.7475 2.5181 2.5325 2.6183 3.1292 3.5323 4.1189 6.2577
a=25.2525 2.5282 2.5426 2.6284 3.1393 3.564 4.1987 6.4267
a=25.7576 2.5383 2.5527 2.6385 3.1494 3.5952 4.2779 6.5943
a=26.2626 2.5484 2.5628 2.6486 3.1595 3.626 4.3563 6.7607
a=26.7677 2.5585 2.5729 2.6587 3.1696 3.6565 4.434 6.9258
a=27.2727 2.5686 3.088 2.6688 3.1797 3.6866 4.5111 7.0897
a=27.7778 3.0837 3.0981 3.1839 3.1898 3.7164 4.5876 6.7473
a=28.2828 3.0938 3.1082 3.194 3.1999 3.7458 4.6633 6.9088
a=28.7879 3.1039 3.1183 3.2041 3.21 3.7749 4.7385 7.0691
a=29.2929 3.114 3.1284 3.2142 3.2201 3.8037 4.8131 7.2283
a=29.798 3.1241 3.1385 3.2243 3.2302 3.8321 4.8871 7.3865
a=30.303 3.1342 3.1486 3.2344 3.2403 3.8603 4.9605 7.5436
a=30.8081 3.1443 3.1587 3.2445 3.2504 3.8881 5.0333 7.6996
a=31.3131 3.1544 3.1688 3.2546 3.2605 3.9156 5.1056 7.8546
a=31.8182 3.1645 3.1789 3.2647 3.2706 3.9428 5.1773 8.0086
a=32.3232 3.1746 3.189 3.2749 3.2807 3.9698 5.2486 8.1616
a=32.8283 3.1847 3.1991 3.285 3.2908 3.9965 5.3193 8.3137
a=33.3333 3.1948 3.2092 3.2951 3.8059 4.0228 5.3894 8.4648
a=33.8384 3.2049 3.2193 3.3052 3.816 4.049 5.4591 8.6151
a=34.3434 3.215 3.2294 3.3153 3.8261 4.0748 5.5284 8.7644
a=34.8485 3.2251 3.2395 3.3254 3.8362 4.1004 5.5971 8.9128
a=35.3535 3.2352 3.2496 3.3355 3.8463 4.1258 5.6654 9.0604
a=35.8586 3.2453 3.2597 3.3456 3.8564 4.1509 5.7332 9.2071
a=36.3636 3.2554 3.2698 3.8607 3.8665 4.1757 5.8005 9.353
a=36.8687 3.2655 3.785 3.8708 3.8766 4.2004 5.8674 10.003
a=37.3737 3.2756 3.7951 3.8809 3.8867 4.7298 5.9339 10.147
a=37.8788 3.7908 3.8052 3.891 3.8968 4.7539 6 10.291
a=38.3838 3.8009 3.8153 3.9011 3.9069 4.7779 6.0656 10.434
a=38.8889 3.811 3.8254 3.9112 3.917 4.8016 6.1308 10.576
a=39.3939 3.8211 3.8355 3.9213 3.9271 4.8251 6.1956 10.717
a=39.899 3.8312 3.8456 3.9314 3.9372 4.8483 6.2601 11.362
a=40.404 3.8413 3.8557 3.9415 3.9473 4.8714 6.3241 11.502
a=40.9091 3.8514 3.8658 3.9516 3.9574 4.8942 6.3878 11.641
a=41.4141 3.8615 3.8759 3.9617 4.4726 4.9169 6.451 11.78
a=41.9192 3.8716 3.886 3.9718 4.4827 4.9393 6.5139 12.422
a=42.4242 3.8817 3.8961 3.9819 4.4928 4.9616 6.5765 12.559
a=42.9293 3.8918 3.9062 3.992 4.5029 4.9837 6.6386 12.696
a=43.4343 3.9019 3.9163 4.0021 4.513 5.0055 6.7005 13.336
a=43.9394 3.912 3.9264 4.0122 4.5231 5.0272 6.762 13.471
a=44.4444 3.9221 3.9365 4.0223 4.5332 5.0487 6.8231 13.606
a=44.9495 3.9322 3.9466 4.0324 4.5433 5.07 7.3889 14.245
a=45.4545 3.9423 3.9567 4.5476 4.5534 5.0911 7.4494 14.883
a=45.9596 3.9524 3.9668 4.5577 4.5635 5.1121 7.5095 15.52
a=46.4646 3.9625 4.482 4.5678 4.5736 5.1329 7.5694 16.157
a=46.9697 3.9726 4.4921 4.5779 4.5837 5.1535 7.6289 16.794
a=47.4747 4.4878 4.5022 4.588 4.5938 5.679 8.1931 17.43
a=47.9798 4.4979 4.5123 4.5981 4.6039 5.6992 8.252 18.065
a=48.4848 4.508 4.5224 4.6082 4.614 5.7193 8.3105 18.7
a=48.9899 4.5181 4.5325 4.6183 4.6241 5.7393 8.8739 19.334
a=49.4949 4.5282 4.5426 4.6284 4.6342 5.7591 9.4369 19.967
a=50 4.5383 4.5527 4.6385 5.1494 5.7787 9.9996 20.6
Table for Consumption at l=1
z1 z2 z3 z4 z5 z6 z7
_______ _______ _______ _______ _______ _______ _______
a=0 0.49495 0.49495 0.49495 0.49495 0.49495 0.49495 0.49495
a=0.505051 0.50505 0.50505 0.50505 0.50505 0.50505 0.50505 0.60473
a=1.0101 0.51515 0.51515 0.51515 0.51515 0.51515 0.51515 0.72755
a=1.51515 0.52525 0.52525 0.52525 0.52525 0.52525 0.68146 0.74187
a=2.0202 1.0404 1.0404 1.0404 1.0404 0.53535 0.91826 1.2009
a=2.52525 1.0505 1.0505 1.0505 1.0505 1.0505 1.1332 1.1145
a=3.0303 1.0606 1.0606 1.0606 1.0606 1.0606 0.82698 1.5034
a=3.53535 1.0707 1.0707 1.0707 1.0707 1.0707 1.0133 1.364
a=4.0404 1.0808 1.0808 1.0808 1.0808 1.0974 1.1893 1.7108
a=4.54545 1.0909 1.0909 1.0909 1.0909 1.1765 1.3569 2.0419
a=5.05051 1.101 1.101 1.101 1.101 1.2516 1.5172 1.8546
a=5.55556 1.1111 1.1111 1.1111 1.1111 1.3234 1.6713 2.1607
a=6.06061 1.1212 1.1212 1.1212 1.1212 1.3922 1.3147 2.4565
a=6.56566 1.1313 1.1313 1.1313 1.1313 1.4583 1.4583 2.2382
a=7.07071 1.1414 1.1414 1.1414 1.1414 1.522 1.5974 2.5168
a=7.57576 1.1515 1.1515 1.1515 1.1515 1.5835 1.7326 2.788
a=8.08081 1.1616 1.1616 1.1616 1.6667 1.6431 1.8641 2.5474
a=8.58586 1.1717 1.1717 1.6768 1.6768 1.7009 1.9923 2.8058
a=9.09091 1.6869 1.6869 1.6869 1.6869 1.757 2.1174 3.0584
a=9.59596 1.697 1.697 1.697 1.697 1.8115 2.2396 2.8007
a=10.101 1.7071 1.7071 1.7071 1.7071 1.8647 2.3592 3.0431
a=10.6061 1.7172 1.7172 1.7172 1.7172 1.9164 2.4764 3.281
a=11.1111 1.7273 1.7273 1.7273 1.7273 1.967 2.0861 3.5146
a=11.6162 1.7374 1.7374 1.7374 1.7374 2.0163 2.1988 3.2391
a=12.1212 1.7475 1.7475 1.7475 1.7475 2.0646 2.3093 3.4649
a=12.6263 1.7576 1.7576 1.7576 1.7576 2.1118 2.418 3.6871
a=13.1313 1.7677 1.7677 1.7677 1.7677 2.158 2.5248 3.9059
a=13.6364 1.7778 1.7778 1.7778 1.7778 2.2033 2.6299 4.1216
a=14.1414 1.7879 1.7879 1.7879 1.7879 2.2476 2.7333 3.8291
a=14.6465 1.798 1.798 1.798 1.798 2.2912 2.8352 4.0387
a=15.1515 1.8081 1.8081 1.8081 2.3131 2.3339 2.9355 4.2456
a=15.6566 1.8182 1.8182 1.8182 2.3232 2.3758 3.0345 4.4498
a=16.1616 1.8283 1.8283 1.8283 2.3333 2.417 3.132 4.1465
a=16.6667 1.8384 1.8384 2.3434 2.3434 2.4575 3.2282 4.3458
a=17.1717 1.8485 2.3535 2.3535 2.3535 2.4973 3.3232 4.5427
a=17.6768 2.3636 2.3636 2.3636 2.3636 2.5365 3.417 4.7374
a=18.1818 2.3737 2.3737 2.3737 2.3737 2.575 3.5096 4.93
a=18.6869 2.3838 2.3838 2.3838 2.3838 2.6129 3.6011 5.1204
a=19.1919 2.3939 2.3939 2.3939 2.3939 2.6502 3.1864 4.8039
a=19.697 2.404 2.404 2.404 2.404 2.687 3.2758 4.9904
a=20.202 2.4141 2.4141 2.4141 2.4141 2.7232 3.3641 5.1751
a=20.7071 2.4242 2.4242 2.4242 2.4242 2.7589 3.4515 5.358
a=21.2121 2.4343 2.4343 2.4343 2.4343 2.7941 3.5379 5.5391
a=21.7172 2.4444 2.4444 2.4444 2.4444 2.8287 3.6234 5.7186
a=22.2222 2.4545 2.4545 2.4545 2.4545 2.8629 3.7081 5.8964
a=22.7273 2.4646 2.4646 2.4646 2.9697 2.8967 3.7918 5.5676
a=23.2323 2.4747 2.4747 2.4747 2.9798 2.93 3.8748 5.7423
a=23.7374 2.4848 2.4848 2.4848 2.9899 2.9628 3.9569 5.9155
a=24.2424 2.4949 2.4949 2.4949 3 3.5003 4.0383 6.0873
a=24.7475 2.5051 2.5051 3.0101 3.0101 3.5323 4.1189 6.2577
a=25.2525 2.5152 2.5152 3.0202 3.0202 3.564 4.1987 6.4267
a=25.7576 2.5253 2.5253 3.0303 3.0303 3.5952 4.2779 6.5943
a=26.2626 2.5354 3.0404 3.0404 3.0404 3.626 4.3563 6.7607
a=26.7677 2.5455 3.0505 3.0505 3.0505 3.6565 4.434 6.9258
a=27.2727 3.0606 3.0606 3.0606 3.0606 3.6866 4.5111 7.0897
a=27.7778 3.0707 3.0707 3.0707 3.0707 3.7164 4.5876 6.7473
a=28.2828 3.0808 3.0808 3.0808 3.0808 3.7458 4.6633 6.9088
a=28.7879 3.0909 3.0909 3.0909 3.0909 3.7749 4.7385 7.0691
a=29.2929 3.101 3.101 3.101 3.101 3.8037 4.8131 7.2283
a=29.798 3.1111 3.1111 3.1111 3.1111 3.8321 4.8871 7.3865
a=30.303 3.1212 3.1212 3.1212 3.1212 3.8603 4.9605 7.5436
a=30.8081 3.1313 3.1313 3.1313 3.6364 3.8881 5.0333 7.6996
a=31.3131 3.1414 3.1414 3.1414 3.6465 3.9156 5.1056 7.8546
a=31.8182 3.1515 3.1515 3.1515 3.6566 3.9428 5.1773 8.0086
a=32.3232 3.1616 3.1616 3.1616 3.6667 3.9698 5.2486 8.1616
a=32.8283 3.1717 3.1717 3.1717 3.6768 3.9965 5.3193 8.3137
a=33.3333 3.1818 3.1818 3.1818 3.6869 4.0228 5.3894 8.4648
a=33.8384 3.1919 3.1919 3.697 3.697 4.049 5.4591 8.6151
a=34.3434 3.202 3.202 3.7071 3.7071 4.0748 5.5284 8.7644
a=34.8485 3.2121 3.2121 3.7172 3.7172 4.1004 5.5971 8.9128
a=35.3535 3.2222 3.2222 3.7273 3.7273 4.1258 5.6654 9.0604
a=35.8586 3.2323 3.7374 3.7374 3.7374 4.1509 5.7332 9.2071
a=36.3636 3.2424 3.7475 3.7475 3.7475 4.1757 5.8005 9.353
a=36.8687 3.7576 3.7576 3.7576 3.7576 4.2004 5.8674 10.003
a=37.3737 3.7677 3.7677 3.7677 3.7677 4.7298 5.9339 10.147
a=37.8788 3.7778 3.7778 3.7778 3.7778 4.7539 6 10.291
a=38.3838 3.7879 3.7879 3.7879 3.7879 4.7779 6.0656 10.434
a=38.8889 3.798 3.798 3.798 3.798 4.8016 6.1308 10.576
a=39.3939 3.8081 3.8081 3.8081 4.3131 4.8251 6.1956 10.717
a=39.899 3.8182 3.8182 3.8182 4.3232 4.8483 6.2601 11.362
a=40.404 3.8283 3.8283 3.8283 4.3333 4.8714 6.3241 11.502
a=40.9091 3.8384 3.8384 3.8384 4.3434 4.8942 6.3878 11.641
a=41.4141 3.8485 3.8485 3.8485 4.3535 4.9169 6.451 11.78
a=41.9192 3.8586 3.8586 3.8586 4.3636 4.9393 6.5139 12.422
a=42.4242 3.8687 3.8687 3.8687 4.3737 4.9616 6.5765 12.559
a=42.9293 3.8788 3.8788 4.3838 4.3838 4.9837 6.6386 12.696
a=43.4343 3.8889 3.8889 4.3939 4.3939 5.0055 6.7005 13.336
a=43.9394 3.899 3.899 4.404 4.404 5.0272 6.762 13.471
a=44.4444 3.9091 3.9091 4.4141 4.4141 5.0487 6.8231 13.606
a=44.9495 3.9192 3.9192 4.4242 4.4242 5.07 7.3889 14.245
a=45.4545 3.9293 4.4343 4.4343 4.4343 5.0911 7.4494 14.883
a=45.9596 3.9394 4.4444 4.4444 4.4444 5.1121 7.5095 15.52
a=46.4646 4.4545 4.4545 4.4545 4.4545 5.1329 7.5694 16.157
a=46.9697 4.4646 4.4646 4.4646 4.4646 5.1535 7.6289 16.794
a=47.4747 4.4747 4.4747 4.4747 4.4747 5.679 8.1931 17.43
a=47.9798 4.4848 4.4848 4.4848 4.9899 5.6992 8.252 18.065
a=48.4848 4.4949 4.4949 4.4949 5 5.7193 8.3105 18.7
a=48.9899 4.5051 4.5051 4.5051 5.0101 5.7393 8.8739 19.334
a=49.4949 4.5152 4.5152 4.5152 5.0202 5.7591 9.4369 19.967
a=50 4.5253 4.5253 4.5253 5.0303 5.7787 9.9996 20.6
Table for Savings/COH at l=0
z1 z2 z3 z4 z5 z6 z7
_______ _______ _______ _______ _______ _______ _______
a=0 NaN NaN NaN NaN NaN NaN NaN
a=0.505051 0.97489 0.9485 0.87516 0.76567 0.62189 0.46167 0.62552
a=1.0101 0.48882 0.48209 0.45382 0.40635 0.68105 0.78621 0.73522
a=1.51515 0.65237 0.64635 0.61615 0.55838 0.71388 0.74776 0.80333
a=2.0202 0.73426 0.72916 0.70124 0.64052 0.73622 0.73334 0.74644
a=2.52525 0.78343 0.77907 0.75422 0.69289 0.75296 0.72783 0.80309
a=3.0303 0.65298 0.64995 0.79061 0.72963 0.76623 0.8301 0.77061
a=3.53535 0.69972 0.69693 0.81697 0.75705 0.77716 0.81772 0.81629
a=4.0404 0.73478 0.73222 0.71734 0.77846 0.78641 0.8094 0.79329
a=4.54545 0.76205 0.75969 0.74593 0.79571 0.79439 0.8037 0.77592
a=5.05051 0.78387 0.78169 0.76892 0.80997 0.80139 0.79978 0.81333
a=5.55556 0.80173 0.79969 0.7878 0.822 0.80762 0.7971 0.79894
a=6.06061 0.81661 0.81471 0.80359 0.83231 0.8132 0.85212 0.78727
a=6.56566 0.8292 0.82742 0.81698 0.84128 0.81826 0.84713 0.81861
a=7.07071 0.84 0.83832 0.82849 0.84915 0.82288 0.84313 0.80821
a=7.57576 0.84935 0.84777 0.83848 0.85615 0.82711 0.83992 0.79941
a=8.08081 0.85754 0.85604 0.84724 0.8624 0.83103 0.83734 0.82633
a=8.58586 0.86476 0.86334 0.85498 0.86805 0.83465 0.83526 0.81819
a=9.09091 0.87119 0.86983 0.86187 0.8218 0.83804 0.83358 0.81109
a=9.59596 0.87693 0.87564 0.86805 0.82907 0.8412 0.83225 0.83469
a=10.101 0.8331 0.83193 0.87361 0.8357 0.84417 0.83119 0.82797
a=10.6061 0.84011 0.83899 0.87865 0.84177 0.84696 0.83036 0.822
a=11.1111 0.84648 0.84541 0.83907 0.84737 0.8496 0.86291 0.81667
a=11.6162 0.85231 0.85127 0.84516 0.85253 0.85209 0.86115 0.83728
a=12.1212 0.85764 0.85664 0.85075 0.85732 0.85446 0.85962 0.8321
a=12.6263 0.86255 0.86159 0.85589 0.86178 0.85671 0.8583 0.82741
a=13.1313 0.86708 0.86615 0.86064 0.86594 0.85886 0.85716 0.82316
a=13.6364 0.87128 0.87038 0.86505 0.86984 0.8609 0.85618 0.8193
a=14.1414 0.87517 0.8743 0.86914 0.87349 0.86286 0.85534 0.83724
a=14.6465 0.8788 0.87795 0.87295 0.87693 0.86473 0.85462 0.83339
a=15.1515 0.88219 0.88137 0.8765 0.88017 0.86652 0.85401 0.82985
a=15.6566 0.88535 0.88456 0.87983 0.88323 0.86825 0.85349 0.8266
a=16.1616 0.88832 0.88755 0.88296 0.88612 0.8699 0.85305 0.84275
a=16.6667 0.89111 0.89036 0.88589 0.88887 0.8715 0.85269 0.83948
a=17.1717 0.89374 0.893 0.88865 0.86361 0.87303 0.8524 0.83645
a=17.6768 0.89621 0.8955 0.89126 0.86684 0.87451 0.85217 0.83363
a=18.1818 0.89855 0.89785 0.89372 0.86989 0.87595 0.85199 0.83101
a=18.6869 0.87427 0.87361 0.89605 0.87279 0.87733 0.85186 0.82856
a=19.1919 0.87706 0.87642 0.8726 0.87553 0.87867 0.87205 0.84281
a=19.697 0.87971 0.87908 0.87535 0.87814 0.87996 0.87153 0.84032
a=20.202 0.88223 0.88161 0.87796 0.88062 0.88121 0.87106 0.83799
a=20.7071 0.88462 0.88402 0.88044 0.88299 0.88243 0.87065 0.8358
a=21.2121 0.8869 0.88631 0.88281 0.88524 0.88361 0.87029 0.83374
a=21.7172 0.88908 0.8885 0.88507 0.8874 0.88476 0.86997 0.83181
a=22.2222 0.89115 0.89058 0.88723 0.88945 0.88587 0.86969 0.83
a=22.7273 0.89313 0.89258 0.88929 0.89142 0.88695 0.86945 0.84257
a=23.2323 0.89503 0.89449 0.89126 0.89331 0.88801 0.86924 0.84069
a=23.7374 0.89685 0.89631 0.89315 0.89511 0.88903 0.86906 0.83892
a=24.2424 0.89859 0.89806 0.89496 0.89685 0.87149 0.86891 0.83724
a=24.7475 0.90026 0.89974 0.8967 0.87898 0.87282 0.86879 0.83565
a=25.2525 0.90186 0.90135 0.89836 0.88096 0.87412 0.86869 0.83415
a=25.7576 0.9034 0.9029 0.89997 0.88287 0.87537 0.86862 0.83273
a=26.2626 0.90488 0.90439 0.90151 0.8847 0.8766 0.86857 0.83138
a=26.7677 0.9063 0.90583 0.90299 0.88647 0.87779 0.86853 0.8301
a=27.2727 0.90768 0.88906 0.90442 0.88817 0.87895 0.86852 0.82889
a=27.7778 0.89117 0.89072 0.88803 0.88981 0.88007 0.86852 0.83974
a=28.2828 0.89277 0.89232 0.88968 0.89139 0.88117 0.86853 0.83846
a=28.7879 0.8943 0.89387 0.89126 0.89292 0.88225 0.86856 0.83724
a=29.2929 0.89579 0.89536 0.89279 0.8944 0.88329 0.86861 0.83608
a=29.798 0.89722 0.8968 0.89427 0.89583 0.88431 0.86866 0.83498
a=30.303 0.89861 0.89819 0.8957 0.89721 0.88531 0.86873 0.83392
a=30.8081 0.89995 0.89954 0.89709 0.89855 0.88628 0.86881 0.83292
a=31.3131 0.90125 0.90084 0.89843 0.89984 0.88724 0.8689 0.83196
a=31.8182 0.9025 0.9021 0.89972 0.9011 0.88817 0.869 0.83105
a=32.3232 0.90372 0.90332 0.90098 0.90231 0.88907 0.8691 0.83018
a=32.8283 0.9049 0.90451 0.9022 0.90349 0.88996 0.86922 0.82935
a=33.3333 0.90604 0.90566 0.90338 0.89005 0.89084 0.86934 0.82856
a=33.8384 0.90715 0.90677 0.90453 0.89137 0.89169 0.86947 0.8278
a=34.3434 0.90823 0.90786 0.90564 0.89266 0.89252 0.86961 0.82708
a=34.8485 0.90927 0.90891 0.90672 0.89391 0.89334 0.86975 0.82639
a=35.3535 0.91029 0.90993 0.90777 0.89512 0.89414 0.86989 0.82573
a=35.8586 0.91128 0.91092 0.90879 0.8963 0.89493 0.87005 0.8251
a=36.3636 0.91224 0.91188 0.8962 0.89745 0.8957 0.8702 0.8245
a=36.8687 0.91317 0.8994 0.89735 0.89857 0.89645 0.87037 0.81456
a=37.3737 0.91408 0.90049 0.89847 0.89966 0.8849 0.87053 0.81412
a=37.8788 0.90189 0.90156 0.89956 0.90072 0.88578 0.8707 0.81371
a=38.3838 0.90293 0.90259 0.90062 0.90175 0.88665 0.87088 0.81331
a=38.8889 0.90393 0.9036 0.90165 0.90276 0.8875 0.87105 0.81294
a=39.3939 0.90491 0.90459 0.90266 0.90374 0.88833 0.87123 0.81258
a=39.899 0.90587 0.90555 0.90364 0.90469 0.88915 0.87142 0.80351
a=40.404 0.9068 0.90648 0.9046 0.90563 0.88995 0.8716 0.80329
a=40.9091 0.90771 0.90739 0.90553 0.90653 0.89074 0.87179 0.80308
a=41.4141 0.90859 0.90828 0.90644 0.89564 0.89151 0.87198 0.80288
a=41.9192 0.90946 0.90915 0.90733 0.89664 0.89227 0.87217 0.79434
a=42.4242 0.9103 0.91 0.9082 0.89763 0.89301 0.87237 0.79426
a=42.9293 0.91113 0.91083 0.90905 0.89859 0.89375 0.87256 0.79419
a=43.4343 0.91193 0.91164 0.90987 0.89953 0.89446 0.87276 0.78602
a=43.9394 0.91272 0.91243 0.91068 0.90044 0.89517 0.87296 0.78606
a=44.4444 0.91349 0.9132 0.91147 0.90134 0.89587 0.87316 0.78609
a=44.9495 0.91424 0.91395 0.91225 0.90222 0.89655 0.86407 0.77828
a=45.4545 0.91498 0.91469 0.90213 0.90307 0.89722 0.86437 0.77062
a=45.9596 0.91569 0.91541 0.90299 0.90391 0.89788 0.86467 0.76312
a=46.4646 0.9164 0.90547 0.90383 0.90473 0.89853 0.86496 0.75577
a=46.9697 0.91709 0.90627 0.90465 0.90554 0.89918 0.86526 0.74857
a=47.4747 0.90733 0.90706 0.90546 0.90632 0.89002 0.85672 0.74151
a=47.9798 0.9081 0.90783 0.90624 0.90709 0.89074 0.8571 0.73459
a=48.4848 0.90885 0.90859 0.90701 0.90785 0.89145 0.85748 0.72781
a=48.9899 0.90959 0.90933 0.90777 0.90859 0.89215 0.84927 0.72115
a=49.4949 0.91031 0.91005 0.90851 0.90931 0.89283 0.84123 0.71462
a=50 0.91102 0.91076 0.90923 0.90023 0.89351 0.83334 0.70821
Table for Savings/COH at l=1
z1 z2 z3 z4 z5 z6 z7
_______ _______ _______ _______ _______ _______ _______
a=0 0.50505 0.50505 0.50505 0.50505 0.50505 0.50505 0.50505
a=0.505051 0.66667 0.66667 0.66667 0.66667 0.66667 0.66667 0.62552
a=1.0101 0.74627 0.74627 0.74627 0.74627 0.74627 0.74627 0.73522
a=1.51515 0.79365 0.79365 0.79365 0.79365 0.79365 0.74776 0.80333
a=2.0202 0.66007 0.66007 0.66007 0.66007 0.82508 0.73334 0.74644
a=2.52525 0.70621 0.70621 0.70621 0.70621 0.70621 0.72783 0.80309
a=3.0303 0.74074 0.74074 0.74074 0.74074 0.74074 0.8301 0.77061
a=3.53535 0.76754 0.76754 0.76754 0.76754 0.76754 0.81772 0.81629
a=4.0404 0.78895 0.78895 0.78895 0.78895 0.78641 0.8094 0.79329
a=4.54545 0.80645 0.80645 0.80645 0.80645 0.79439 0.8037 0.77592
a=5.05051 0.82102 0.82102 0.82102 0.82102 0.80139 0.79978 0.81333
a=5.55556 0.83333 0.83333 0.83333 0.83333 0.80762 0.7971 0.79894
a=6.06061 0.84388 0.84388 0.84388 0.84388 0.8132 0.85212 0.78727
a=6.56566 0.85302 0.85302 0.85302 0.85302 0.81826 0.84713 0.81861
a=7.07071 0.86101 0.86101 0.86101 0.86101 0.82288 0.84313 0.80821
a=7.57576 0.86806 0.86806 0.86806 0.86806 0.82711 0.83992 0.79941
a=8.08081 0.87432 0.87432 0.87432 0.81967 0.83103 0.83734 0.82633
a=8.58586 0.87992 0.87992 0.82816 0.82816 0.83465 0.83526 0.81819
a=9.09091 0.83579 0.83579 0.83579 0.83579 0.83804 0.83358 0.81109
a=9.59596 0.8427 0.8427 0.8427 0.8427 0.8412 0.83225 0.83469
a=10.101 0.84897 0.84897 0.84897 0.84897 0.84417 0.83119 0.82797
a=10.6061 0.8547 0.8547 0.8547 0.8547 0.84696 0.83036 0.822
a=11.1111 0.85995 0.85995 0.85995 0.85995 0.8496 0.86291 0.81667
a=11.6162 0.86478 0.86478 0.86478 0.86478 0.85209 0.86115 0.83728
a=12.1212 0.86924 0.86924 0.86924 0.86924 0.85446 0.85962 0.8321
a=12.6263 0.87336 0.87336 0.87336 0.87336 0.85671 0.8583 0.82741
a=13.1313 0.87719 0.87719 0.87719 0.87719 0.85886 0.85716 0.82316
a=13.6364 0.88076 0.88076 0.88076 0.88076 0.8609 0.85618 0.8193
a=14.1414 0.88409 0.88409 0.88409 0.88409 0.86286 0.85534 0.83724
a=14.6465 0.8872 0.8872 0.8872 0.8872 0.86473 0.85462 0.83339
a=15.1515 0.89012 0.89012 0.89012 0.85942 0.86652 0.85401 0.82985
a=15.6566 0.89286 0.89286 0.89286 0.8631 0.86825 0.85349 0.8266
a=16.1616 0.89544 0.89544 0.89544 0.86655 0.8699 0.85305 0.84275
a=16.6667 0.89787 0.89787 0.86981 0.86981 0.8715 0.85269 0.83948
a=17.1717 0.90016 0.87289 0.87289 0.87289 0.87303 0.8524 0.83645
a=17.6768 0.8758 0.8758 0.8758 0.8758 0.87451 0.85217 0.83363
a=18.1818 0.87855 0.87855 0.87855 0.87855 0.87595 0.85199 0.83101
a=18.6869 0.88117 0.88117 0.88117 0.88117 0.87733 0.85186 0.82856
a=19.1919 0.88365 0.88365 0.88365 0.88365 0.87867 0.87205 0.84281
a=19.697 0.88602 0.88602 0.88602 0.88602 0.87996 0.87153 0.84032
a=20.202 0.88827 0.88827 0.88827 0.88827 0.88121 0.87106 0.83799
a=20.7071 0.89041 0.89041 0.89041 0.89041 0.88243 0.87065 0.8358
a=21.2121 0.89246 0.89246 0.89246 0.89246 0.88361 0.87029 0.83374
a=21.7172 0.89442 0.89442 0.89442 0.89442 0.88476 0.86997 0.83181
a=22.2222 0.89629 0.89629 0.89629 0.89629 0.88587 0.86969 0.83
a=22.7273 0.89808 0.89808 0.89808 0.87719 0.88695 0.86945 0.84257
a=23.2323 0.8998 0.8998 0.8998 0.87935 0.88801 0.86924 0.84069
a=23.7374 0.90144 0.90144 0.90144 0.88141 0.88903 0.86906 0.83892
a=24.2424 0.90302 0.90302 0.90302 0.88339 0.87149 0.86891 0.83724
a=24.7475 0.90454 0.90454 0.8853 0.8853 0.87282 0.86879 0.83565
a=25.2525 0.906 0.906 0.88713 0.88713 0.87412 0.86869 0.83415
a=25.7576 0.90741 0.90741 0.88889 0.88889 0.87537 0.86862 0.83273
a=26.2626 0.90876 0.89059 0.89059 0.89059 0.8766 0.86857 0.83138
a=26.7677 0.91006 0.89222 0.89222 0.89222 0.87779 0.86853 0.8301
a=27.2727 0.8938 0.8938 0.8938 0.8938 0.87895 0.86852 0.82889
a=27.7778 0.89532 0.89532 0.89532 0.89532 0.88007 0.86852 0.83974
a=28.2828 0.89679 0.89679 0.89679 0.89679 0.88117 0.86853 0.83846
a=28.7879 0.8982 0.8982 0.8982 0.8982 0.88225 0.86856 0.83724
a=29.2929 0.89957 0.89957 0.89957 0.89957 0.88329 0.86861 0.83608
a=29.798 0.9009 0.9009 0.9009 0.9009 0.88431 0.86866 0.83498
a=30.303 0.90218 0.90218 0.90218 0.90218 0.88531 0.86873 0.83392
a=30.8081 0.90343 0.90343 0.90343 0.88785 0.88628 0.86881 0.83292
a=31.3131 0.90463 0.90463 0.90463 0.8893 0.88724 0.8689 0.83196
a=31.8182 0.9058 0.9058 0.9058 0.8907 0.88817 0.869 0.83105
a=32.3232 0.90693 0.90693 0.90693 0.89206 0.88907 0.8691 0.83018
a=32.8283 0.90803 0.90803 0.90803 0.89338 0.88996 0.86922 0.82935
a=33.3333 0.90909 0.90909 0.90909 0.89466 0.89084 0.86934 0.82856
a=33.8384 0.91013 0.91013 0.8959 0.8959 0.89169 0.86947 0.8278
a=34.3434 0.91113 0.91113 0.89711 0.89711 0.89252 0.86961 0.82708
a=34.8485 0.91211 0.91211 0.89829 0.89829 0.89334 0.86975 0.82639
a=35.3535 0.91306 0.91306 0.89943 0.89943 0.89414 0.86989 0.82573
a=35.8586 0.91398 0.90054 0.90054 0.90054 0.89493 0.87005 0.8251
a=36.3636 0.91488 0.90162 0.90162 0.90162 0.8957 0.8702 0.8245
a=36.8687 0.90267 0.90267 0.90267 0.90267 0.89645 0.87037 0.81456
a=37.3737 0.90369 0.90369 0.90369 0.90369 0.8849 0.87053 0.81412
a=37.8788 0.90469 0.90469 0.90469 0.90469 0.88578 0.8707 0.81371
a=38.3838 0.90566 0.90566 0.90566 0.90566 0.88665 0.87088 0.81331
a=38.8889 0.90661 0.90661 0.90661 0.90661 0.8875 0.87105 0.81294
a=39.3939 0.90753 0.90753 0.90753 0.89527 0.88833 0.87123 0.81258
a=39.899 0.90843 0.90843 0.90843 0.89632 0.88915 0.87142 0.80351
a=40.404 0.90931 0.90931 0.90931 0.89734 0.88995 0.8716 0.80329
a=40.9091 0.91017 0.91017 0.91017 0.89835 0.89074 0.87179 0.80308
a=41.4141 0.911 0.911 0.911 0.89932 0.89151 0.87198 0.80288
a=41.9192 0.91182 0.91182 0.91182 0.90028 0.89227 0.87217 0.79434
a=42.4242 0.91262 0.91262 0.91262 0.90121 0.89301 0.87237 0.79426
a=42.9293 0.9134 0.9134 0.90212 0.90212 0.89375 0.87256 0.79419
a=43.4343 0.91416 0.91416 0.90301 0.90301 0.89446 0.87276 0.78602
a=43.9394 0.9149 0.9149 0.90388 0.90388 0.89517 0.87296 0.78606
a=44.4444 0.91563 0.91563 0.90473 0.90473 0.89587 0.87316 0.78609
a=44.9495 0.91634 0.91634 0.90556 0.90556 0.89655 0.86407 0.77828
a=45.4545 0.91704 0.90638 0.90638 0.90638 0.89722 0.86437 0.77062
a=45.9596 0.91772 0.90717 0.90717 0.90717 0.89788 0.86467 0.76312
a=46.4646 0.90795 0.90795 0.90795 0.90795 0.89853 0.86496 0.75577
a=46.9697 0.90872 0.90872 0.90872 0.90872 0.89918 0.86526 0.74857
a=47.4747 0.90946 0.90946 0.90946 0.90946 0.89002 0.85672 0.74151
a=47.9798 0.91019 0.91019 0.91019 0.90008 0.89074 0.8571 0.73459
a=48.4848 0.91091 0.91091 0.91091 0.9009 0.89145 0.85748 0.72781
a=48.9899 0.91161 0.91161 0.91161 0.9017 0.89215 0.84927 0.72115
a=49.4949 0.9123 0.9123 0.9123 0.90249 0.89283 0.84123 0.71462
a=50 0.91298 0.91298 0.91298 0.90326 0.89351 0.83334 0.70821
Plot Value function and Policy Function
% 3D Value function - l=0 if(bl_plot) figure(1) surf(ar_z,ar_a,mt_vf(:,:,1)); title('Value Function, l=0'); xlabel('Productivity') ; ylabel('Assets'); zlabel('Value'); xt = get(gca, 'XTick'); set(gca, 'XTick', xt, 'XTickLabel', xt*fl_ahi/it_agridno); if(bl_saveimg) saveas(gcf, '/Users/sidhantkhanna/Documents/GitHub/BKS modified/code/Households/figures/hhsave_VFI_3OC/Valuefunction_3D_0.png') end figure(2) surf(ar_z,ar_a,mt_vf(:,:,2)); title('Value Function, l=1'); xlabel('Productivity') ; ylabel('Assets'); zlabel('Value'); xt = get(gca, 'XTick'); set(gca, 'XTick', xt, 'XTickLabel', xt*fl_ahi/it_agridno); if(bl_saveimg) saveas(gcf, '/Users/sidhantkhanna/Documents/GitHub/BKS modified/code/Households/figures/hhsave_VFI_3OC/Valuefunction_3D_1.png') end % 2D Value function figure(3) %chart = plot(mt_vf); for i = 1: it_zgridno a(i) = plot((mt_vf(:,i,1)), '.'); b{i} = sprintf(['Z shock =', num2str(ar_z(i))]); hold on end clr = jet(numel(a)); for m = 1:numel(a) set(a(m),'Color',clr(m,:)) end legend(b,'Location','South') %clr = jet(numel(chart)); %for m = 1:numel(chart) % set(chart(m),'Color',clr(m,:)) %end title('Value Function, l=0'); xlabel('Asset') ; ylabel('Value'); xt = get(gca, 'XTick'); set(gca, 'XTick', xt, 'XTickLabel', xt*fl_ahi/it_agridno); if(bl_saveimg) saveas(gcf, '/Users/sidhantkhanna/Documents/GitHub/BKS modified/code/Households/figures/hhsave_VFI_3OC/Valuefunction_2D_0.png') end figure(4) %chart = plot(mt_vf); for i = 1: it_zgridno a(i) = plot((mt_vf(:,i,2)),'.'); b{i} = sprintf(['Z shock =', num2str(ar_z(i))]); hold on end clr = jet(numel(a)); for m = 1:numel(a) set(a(m),'Color',clr(m,:)) end legend(b,'Location','South') %clr = jet(numel(chart)); %for m = 1:numel(chart) % set(chart(m),'Color',clr(m,:)) %end title('Value Function, l=1'); xlabel('Asset') ; ylabel('Value'); xt = get(gca, 'XTick'); set(gca, 'XTick', xt, 'XTickLabel', xt*fl_ahi/it_agridno); if(bl_saveimg) saveas(gcf, '/Users/sidhantkhanna/Documents/GitHub/BKS modified/code/Households/figures/hhsave_VFI_3OC/Valuefunction_2D_1.png') end figure(5) for i = 1: it_zgridno a(i) = plot((mt_pf(:,i,1)),'.'); b{i} = sprintf(['Z shock =', num2str(ar_z(i))]); hold on end clr = jet(numel(a)); for m = 1:numel(a) set(a(m),'Color',clr(m,:)) end legend(b,'Location','north') title('Policy Function,l=0'); xlabel('Asset') ; ylabel('Asset Choice'); xt = get(gca, 'XTick'); set(gca, 'XTick', xt, 'XTickLabel', xt*fl_ahi/it_agridno); if(bl_saveimg) saveas(gcf, '/Users/sidhantkhanna/Documents/GitHub/BKS modified/code/Households/figures/hhsave_VFI_3OC/Policyfunction_0.png') end figure(6) for i = 1: it_zgridno a(i) = plot((mt_pf(:,i,2)),'.'); b{i} = sprintf(['Z shock =', num2str(ar_z(i))]); hold on end clr = jet(numel(a)); for m = 1:numel(a) set(a(m),'Color',clr(m,:)) end legend(b,'Location','north') title('Policy Function, l=1'); xlabel('Asset') ; ylabel('Asset choice'); xt = get(gca, 'XTick'); set(gca, 'XTick', xt, 'XTickLabel', xt*fl_ahi/it_agridno); if(bl_saveimg) saveas(gcf, '/Users/sidhantkhanna/Documents/GitHub/BKS modified/code/Households/figures/hhsave_VFI_3OC/Policyfunction_1.png') end figure(7) for i = 1: it_zgridno a(i) = plot((mt_con(:,i,1)),'-o'); b{i} = sprintf(['Z shock =', num2str(ar_z(i))]); hold on end clr = jet(numel(a)); for m = 1:numel(a) set(a(m),'Color',clr(m,:)) end legend(b,'Location','north') title('Consumption,l=0'); xlabel('Asset') ; ylabel('Consumption'); xt = get(gca, 'XTick'); set(gca, 'XTick', xt, 'XTickLabel', xt*fl_ahi/it_agridno); if(bl_saveimg) saveas(gcf, '/Users/sidhantkhanna/Documents/GitHub/BKS modified/code/Households/figures/hhsave_VFI_3OC/Consumption_0.png') end figure(8) for i = 1: it_zgridno a(i) = plot((mt_con(:,i,2)),'-o'); b{i} = sprintf(['Z shock =', num2str(ar_z(i))]); hold on end clr = jet(numel(a)); for m = 1:numel(a) set(a(m),'Color',clr(m,:)) end legend(b,'Location','north') title('Consumption, l=1'); xlabel('Asset') ; ylabel('Consumption'); xt = get(gca, 'XTick'); set(gca, 'XTick', xt, 'XTickLabel', xt*fl_ahi/it_agridno); if(bl_saveimg) saveas(gcf, '/Users/sidhantkhanna/Documents/GitHub/BKS modified/code/Households/figures/hhsave_VFI_3OC/Consumption_1.png') end figure(9) for i = 1: it_zgridno a(i) = plot(mt_savefraccoh(:,i,1),'-o'); b{i} = sprintf(['Z shock =', num2str(ar_z(i))]); hold on end clr = jet(numel(a)); for m = 1:numel(a) set(a(m),'Color',clr(m,:)) end legend(b,'Location','South') title('SavingsfracCOH,l=0'); xlabel('Asset') ; ylabel('SavingsfracCOH'); ylim([0 1]) xt = get(gca, 'XTick'); set(gca, 'XTick', xt, 'XTickLabel', xt*fl_ahi/it_agridno); if(bl_saveimg) saveas(gcf, '/Users/sidhantkhanna/Documents/GitHub/BKS modified/code/Households/figures/hhsave_VFI_3OC/Savings_COH_0.png') end figure(10) for i = 1: it_zgridno a(i) = plot(mt_savefraccoh(:,i,2),'-o'); b{i} = sprintf(['Z shock =', num2str(ar_z(i))]); hold on end clr = jet(numel(a)); for m = 1:numel(a) set(a(m),'Color',clr(m,:)) end legend(b,'Location','South') title('SavingsfracCOH, l=1'); xlabel('Asset') ; ylabel('SavingsfracCOH'); ylim([0 1]) xt = get(gca, 'XTick'); set(gca, 'XTick', xt, 'XTickLabel', xt*fl_ahi/it_agridno); if(bl_saveimg) saveas(gcf, '/Users/sidhantkhanna/Documents/GitHub/BKS modified/code/Households/figures/hhsave_VFI_3OC/Savings_COH_1.png') end end
Plot Value function and Policy Function
if(bl_profile) profile off; profile viewer; st_file_name = '/Users/sidhantkhanna/Documents/GitHub/BKS modified/code/Profile/Households/hhsave_VFI_3OC'; profsave(profile('info'), st_file_name); end