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_o, mt_k, mt_l]=hhsave_VFI_powerspace(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_mu_n, fl_sig_n,fl_rho_n,fl_lambda_n, 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(0,fl_ahi,it_agridno); ar_a = ff_saveborr_grid(fl_alo, fl_ahi, it_agridno, 'grid_powerspace'); ar_a = ar_a'; fl_phi = 0.2; 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.02; 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); ar_z = ar_z'; % P1=mt_trans_z^1000; % sd=P1(1,:); % el=sum(sd.*ar_z ); % ar_z =ar_z ./el; ar_n = [0,1]; mt_trans_n = [0.7,0.3;0.7,0.3]; [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_evf(:,j,:) = squeeze(mt_evf(:,j,:))*mt_trans_n'; end mt_ev1 = repmat(mt_evf, [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 -112.62 -103.34 -95.188 -87.935 -81.403 -75.225 -68.518
a=5.15305e-05 -103.69 -96.642 -90.135 -84.101 -78.16 -71.688 -66.048
a=0.000412244 -83.263 -78.453 -73.951 -69.719 -65.72 -61.938 -58.357
a=0.00139132 -68.111 -64.554 -61.2 -58.029 -55.017 -52.154 -49.431
a=0.00329795 -57.592 -54.801 -52.132 -49.589 -47.162 -44.847 -42.64
a=0.00644131 -49.868 -47.663 -45.477 -43.36 -41.321 -39.367 -37.499
a=0.0111306 -43.748 -42.136 -40.359 -38.573 -36.824 -35.132 -33.506
a=0.017675 -38.873 -37.556 -36.235 -34.746 -33.24 -31.756 -30.315
a=0.0263836 -35.035 -33.865 -32.715 -31.565 -30.289 -28.985 -27.698
a=0.0375657 -31.954 -30.889 -29.844 -28.814 -27.782 -26.654 -25.501
a=0.0515305 -29.423 -28.44 -27.475 -26.525 -25.581 -24.645 -23.618
a=0.0685871 -27.304 -26.386 -25.486 -24.599 -23.717 -22.846 -21.978
a=0.0890447 -25.495 -24.634 -23.788 -22.953 -22.122 -21.302 -20.497
a=0.113213 -23.898 -23.084 -22.287 -21.501 -20.719 -19.946 -19.188
a=0.1414 -22.498 -21.723 -20.963 -20.214 -19.468 -18.731 -18.008
a=0.173915 -21.214 -20.516 -19.791 -19.074 -18.359 -17.652 -16.959
a=0.211069 -20.038 -19.406 -18.74 -18.054 -17.367 -16.687 -16.02
a=0.253169 -18.963 -18.378 -17.746 -17.088 -16.428 -15.775 -15.134
a=0.300526 -17.965 -17.412 -16.837 -16.21 -15.573 -14.941 -14.32
a=0.353448 -17.044 -16.518 -15.98 -15.399 -14.788 -14.18 -13.58
a=0.412244 -16.187 -15.686 -15.176 -14.617 -14.023 -13.432 -12.85
a=0.477224 -15.383 -14.902 -14.412 -13.897 -13.331 -12.757 -12.19
a=0.548697 -14.634 -14.17 -13.701 -13.205 -12.652 -12.096 -11.547
a=0.626972 -13.919 -13.469 -13.015 -12.545 -12.023 -11.484 -10.949
a=0.712358 -13.245 -12.809 -12.37 -11.918 -11.391 -10.86 -10.336
a=0.805164 -12.597 -12.172 -11.746 -11.307 -10.8 -10.282 -9.7711
a=0.9057 -11.976 -11.562 -11.147 -10.72 -10.225 -9.7169 -9.2149
a=1.01427 -11.381 -10.977 -10.571 -10.154 -9.6685 -9.1732 -8.684
a=1.1312 -10.802 -10.407 -10.01 -9.6034 -9.146 -8.6645 -8.1858
a=1.25678 -10.247 -9.861 -9.4738 -9.0775 -8.6195 -8.1438 -7.6722
a=1.39132 -9.7088 -9.3303 -8.9507 -8.5618 -8.1041 -7.6356 -7.1726
a=1.53515 -9.181 -8.8099 -8.4377 -8.0557 -7.6144 -7.1569 -6.7041
a=1.68855 -8.6685 -8.3045 -7.9397 -7.5656 -7.1399 -6.6898 -6.2431
a=1.85185 -8.175 -7.8179 -7.4595 -7.0921 -6.6635 -6.2198 -5.7804
a=2.02536 -7.6879 -7.3365 -6.9842 -6.6231 -6.2043 -5.7693 -5.3382
a=2.20937 -7.2106 -6.865 -6.519 -6.1645 -5.7583 -5.3304 -4.9056
a=2.40421 -6.7463 -6.4065 -6.0664 -5.7185 -5.3249 -4.9049 -4.4877
a=2.61017 -6.2946 -5.9603 -5.626 -5.2845 -4.8963 -4.4805 -4.0672
a=2.82758 -5.8553 -5.5254 -5.1957 -4.859 -4.4714 -4.062 -3.6558
a=3.05674 -5.4189 -5.0937 -4.7689 -4.4372 -4.06 -3.6592 -3.2606
a=3.29795 -4.9893 -4.6688 -4.3491 -4.0233 -3.6623 -3.2707 -2.8794
a=3.55153 -4.5689 -4.2533 -3.9388 -3.6198 -3.2769 -2.894 -2.5093
a=3.81779 -4.1595 -3.8487 -3.5396 -3.2278 -2.9013 -2.5254 -2.146
a=4.09704 -3.7619 -3.4557 -3.1519 -2.847 -2.5312 -2.1598 -1.784
a=4.38957 -3.3756 -3.0736 -2.7745 -2.4752 -2.1596 -1.7901 -1.4167
a=4.69572 -2.9968 -2.6981 -2.4026 -2.1069 -1.7927 -1.4262 -1.0562
a=5.01577 -2.622 -2.3266 -2.0345 -1.7421 -1.4342 -1.0723 -0.7067
a=5.35005 -2.2524 -1.9603 -1.6715 -1.3827 -1.0838 -0.72817 -0.36756
a=5.69886 -1.8888 -1.5998 -1.3145 -1.0297 -0.73759 -0.3913 -0.036734
a=6.06251 -1.532 -1.2461 -0.96409 -0.68333 -0.39698 -0.062282 0.28559
a=6.44131 -1.1825 -0.89961 -0.62087 -0.34401 -0.062911 0.26118 0.61521
a=6.83557 -0.84072 -0.56069 -0.28513 -0.011983 0.26433 0.57631 0.93526
a=7.2456 -0.50687 -0.22962 0.042887 0.31251 0.58423 0.88301 1.2458
a=7.67171 -0.18117 0.093432 0.36304 0.62929 0.89652 1.183 1.5487
a=8.1142 0.1363 0.4084 0.67524 0.93829 1.2013 1.487 1.8538
a=8.57339 0.44552 0.71526 0.97951 1.2397 1.5001 1.7851 2.1517
a=9.04958 0.74798 1.016 1.2783 1.5365 1.7951 2.0771 2.4424
a=9.54309 1.0457 1.312 1.5726 1.8289 2.0857 2.3642 2.7268
a=10.0542 1.3385 1.6032 1.862 2.1165 2.3717 2.649 3.0067
a=10.5833 1.6261 1.8893 2.1464 2.3992 2.653 2.9303 3.2815
a=11.1306 1.9084 2.1701 2.4257 2.6771 2.9305 3.2078 3.5512
a=11.6964 2.1876 2.4456 2.6998 2.9501 3.2043 3.4812 3.8156
a=12.2812 2.4629 2.7158 2.9687 3.2183 3.4738 3.7502 4.0746
a=12.885 2.734 2.9806 3.2324 3.4816 3.7387 4.0145 4.3279
a=13.5084 3.001 3.2402 3.4909 3.7398 3.9987 4.2739 4.5757
a=14.1516 3.2635 3.4944 3.7441 3.993 4.2536 4.5282 4.8199
a=14.8148 3.5217 3.7439 3.9921 4.2411 4.5033 4.7774 5.0623
a=15.4985 3.7755 3.9919 4.2349 4.484 4.7478 5.0214 5.3017
a=16.2028 4.025 4.2351 4.4728 4.7218 4.987 5.2603 5.5374
a=16.9282 4.2705 4.4753 4.7057 4.9545 5.221 5.4941 5.7688
a=17.675 4.5119 4.7113 4.9338 5.1821 5.4497 5.7226 5.9957
a=18.4433 4.7496 4.9437 5.1571 5.4048 5.6731 5.946 6.2178
a=19.2337 4.9833 5.1721 5.3799 5.6226 5.8914 6.1642 6.435
a=20.0462 5.2133 5.3966 5.5982 5.8361 6.1047 6.3773 6.6473
a=20.8814 5.4396 5.6177 5.8142 6.0451 6.3131 6.5854 6.8547
a=21.7394 5.6623 5.8352 6.0264 6.2498 6.5167 6.7884 7.0571
a=22.6207 5.8814 6.0494 6.2356 6.4503 6.7156 6.9866 7.2547
a=23.5254 6.097 6.2602 6.4412 6.65 6.9098 7.18 7.4474
a=24.4539 6.3092 6.4677 6.6436 6.8461 7.1 7.3688 7.6354
a=25.4065 6.518 6.6719 6.8431 7.0407 7.2862 7.5532 7.8188
a=26.3836 6.7236 6.8731 7.0396 7.2318 7.4687 7.7332 7.9978
a=27.3854 6.926 7.0713 7.2337 7.4209 7.6481 7.9092 8.1724
a=28.4123 7.1253 7.2669 7.425 7.6068 7.8268 8.0813 8.3428
a=29.4645 7.3219 7.4598 7.6138 7.7905 8.0026 8.2501 8.5092
a=30.5423 7.5158 7.6501 7.7998 7.9714 8.1771 8.4157 8.672
a=31.6462 7.707 7.8378 7.9834 8.1501 8.3489 8.5785 8.8313
a=32.7763 7.8957 8.0228 8.1644 8.3261 8.519 8.7392 8.9873
a=33.933 8.0816 8.2053 8.343 8.5 8.6866 8.8991 9.1403
a=35.1166 8.265 8.3853 8.519 8.6714 8.8523 9.0569 9.2909
a=36.3274 8.4458 8.5627 8.6927 8.8405 9.0157 9.2136 9.439
a=37.5657 8.6241 8.7377 8.8639 9.0074 9.1772 9.3682 9.5851
a=38.8319 8.7998 8.9102 9.0328 9.1721 9.3365 9.5215 9.7291
a=40.1262 8.9731 9.0803 9.1993 9.3345 9.494 9.6728 9.8727
a=41.4489 9.1439 9.2481 9.3636 9.4948 9.6494 9.8226 10.015
a=42.8004 9.3123 9.4134 9.5256 9.6529 9.8028 9.9705 10.156
a=44.181 9.4782 9.5765 9.6854 9.8089 9.9543 10.117 10.296
a=45.5909 9.6418 9.7372 9.843 9.9629 10.104 10.261 10.435
a=47.0305 9.8031 9.8958 9.9984 10.115 10.252 10.404 10.572
a=48.5001 9.9621 10.052 10.152 10.265 10.398 10.546 10.708
a=50 10.119 10.206 10.303 10.413 10.542 10.685 10.842
Table for value function at l=1
z1 z2 z3 z4 z5 z6 z7
________ ________ _________ ________ ________ ________ _________
a=0 -12.808 -12.482 -12.152 -11.807 -11.424 -11.029 -10.636
a=5.15305e-05 -12.808 -12.482 -12.151 -11.807 -11.423 -11.029 -10.636
a=0.000412244 -12.806 -12.479 -12.149 -11.805 -11.421 -11.027 -10.633
a=0.00139132 -12.799 -12.473 -12.143 -11.798 -11.415 -11.02 -10.627
a=0.00329795 -12.787 -12.46 -12.13 -11.786 -11.402 -11.008 -10.614
a=0.00644131 -12.766 -12.44 -12.11 -11.766 -11.382 -10.987 -10.594
a=0.0111306 -12.737 -12.411 -12.08 -11.736 -11.352 -10.958 -10.565
a=0.017675 -12.697 -12.37 -12.04 -11.696 -11.312 -10.918 -10.524
a=0.0263836 -12.645 -12.319 -11.989 -11.644 -11.261 -10.866 -10.473
a=0.0375657 -12.582 -12.256 -11.926 -11.581 -11.198 -10.803 -10.41
a=0.0515305 -12.508 -12.182 -11.852 -11.507 -11.124 -10.729 -10.336
a=0.0685871 -12.424 -12.097 -11.767 -11.423 -11.039 -10.645 -10.251
a=0.0890447 -12.32 -12.002 -11.674 -11.329 -10.946 -10.551 -10.158
a=0.113213 -12.169 -11.85 -11.529 -11.194 -10.821 -10.436 -10.052
a=0.1414 -12.013 -11.695 -11.373 -11.038 -10.666 -10.281 -9.8964
a=0.173915 -11.857 -11.538 -11.216 -10.881 -10.509 -10.124 -9.7398
a=0.211069 -11.69 -11.378 -11.061 -10.726 -10.354 -9.9689 -9.5847
a=0.253169 -11.464 -11.152 -10.837 -10.509 -10.143 -9.7637 -9.3848
a=0.300526 -11.254 -10.942 -10.627 -10.299 -9.9335 -9.5542 -9.1754
a=0.353448 -11.018 -10.713 -10.405 -10.084 -9.7267 -9.3567 -8.982
a=0.412244 -10.757 -10.452 -10.144 -9.8235 -9.4658 -9.0958 -8.7262
a=0.477224 -10.515 -10.217 -9.9152 -9.5951 -9.2374 -8.8674 -8.4978
a=0.548697 -10.216 -9.918 -9.6166 -9.3043 -8.9627 -8.6033 -8.2412
a=0.626972 -9.9455 -9.6524 -9.3565 -9.0495 -8.7087 -8.3521 -7.991
a=0.712358 -9.6306 -9.3374 -9.0415 -8.7345 -8.3937 -8.0371 -7.6791
a=0.805164 -9.3192 -9.0314 -8.7409 -8.4387 -8.0995 -7.7474 -7.3953
a=0.9057 -9.0098 -8.722 -8.4315 -8.1293 -7.7901 -7.438 -7.0859
a=1.01427 -8.6653 -8.3829 -8.0979 -7.8017 -7.4727 -7.1292 -6.7848
a=1.1312 -8.3478 -8.0705 -7.7901 -7.4991 -7.1782 -6.8396 -6.4975
a=1.25678 -8.0121 -7.7348 -7.4544 -7.1634 -6.8425 -6.5039 -6.1634
a=1.39132 -7.6535 -7.3803 -7.1045 -6.818 -6.4983 -6.1639 -5.8285
a=1.53515 -7.3067 -7.038 -6.7669 -6.4856 -6.173 -5.8454 -5.5163
a=1.68855 -6.9764 -6.7119 -6.4421 -6.1608 -5.8482 -5.5206 -5.1915
a=1.85185 -6.6101 -6.3456 -6.0791 -5.803 -5.4976 -5.1753 -4.8507
a=2.02536 -6.2493 -5.9891 -5.727 -5.456 -5.1581 -4.8419 -4.5231
a=2.20937 -5.8999 -5.6429 -5.3842 -5.1168 -4.8221 -4.5085 -4.192
a=2.40421 -5.5529 -5.2995 -5.0445 -4.7809 -4.4889 -4.1803 -3.8693
a=2.61017 -5.2115 -4.958 -4.7031 -4.4394 -4.1475 -3.8388 -3.5279
a=2.82758 -4.8413 -4.5915 -4.3406 -4.0817 -3.7978 -3.4963 -3.1915
a=3.05674 -4.478 -4.2321 -3.9854 -3.7324 -3.4593 -3.1658 -2.867
a=3.29795 -4.1237 -3.8817 -3.6395 -3.3927 -3.131 -2.8449 -2.5515
a=3.55153 -3.7793 -3.541 -3.3032 -3.0624 -2.8104 -2.5301 -2.2408
a=3.81779 -3.4443 -3.2094 -2.9755 -2.7394 -2.4933 -2.2164 -1.9298
a=4.09704 -3.1151 -2.8828 -2.6517 -2.4184 -2.1729 -1.8971 -1.6117
a=4.38957 -2.7817 -2.5494 -2.3183 -2.085 -1.8395 -1.5636 -1.2782
a=4.69572 -2.4386 -2.2088 -1.9804 -1.7497 -1.5056 -1.2318 -0.949
a=5.01577 -2.0987 -1.8715 -1.6458 -1.4179 -1.1782 -0.90822 -0.62881
a=5.35005 -1.7631 -1.5385 -1.3154 -1.091 -0.85735 -0.59252 -0.31731
a=5.69886 -1.4326 -1.2105 -0.99021 -0.76917 -0.54063 -0.28241 -0.012475
a=6.06251 -1.1079 -0.88821 -0.67063 -0.45291 -0.22901 0.021461 0.28559
a=6.44131 -0.78928 -0.57204 -0.35711 -0.14255 0.077129 0.3213 0.61521
a=6.83557 -0.47726 -0.26233 -0.049961 0.16159 0.37723 0.61425 0.93526
a=7.2456 -0.17201 0.040697 0.25061 0.45929 0.67094 0.90016 1.2458
a=7.67171 0.12629 0.33689 0.54447 0.75039 0.95826 1.183 1.5487
a=8.1142 0.41759 0.62618 0.83157 1.0351 1.2417 1.487 1.8538
a=8.57339 0.70651 0.9151 1.1205 1.324 1.5297 1.7851 2.1517
a=9.04958 0.99752 1.2047 1.4086 1.6104 1.8148 2.0771 2.4424
a=9.54309 1.2848 1.4906 1.6929 1.8933 2.0963 2.3642 2.7268
a=10.0542 1.5679 1.7724 1.9733 2.1722 2.374 2.649 3.0067
a=10.5833 1.8468 2.0499 2.2495 2.447 2.653 2.9303 3.2815
a=11.1306 2.1211 2.323 2.5213 2.7177 2.9305 3.2078 3.5512
a=11.6964 2.3907 2.5915 2.7887 2.9843 3.2043 3.4812 3.8156
a=12.2812 2.6567 2.8554 3.0514 3.2466 3.4738 3.7502 4.0746
a=12.885 2.9188 3.1145 3.3096 3.5046 3.7387 4.0145 4.3279
a=13.5084 3.1768 3.369 3.5632 3.7582 3.9987 4.2739 4.5757
a=14.1516 3.4308 3.6188 3.8121 4.0072 4.2536 4.5282 4.8199
a=14.8148 3.6805 3.8638 4.0563 4.2517 4.5033 4.7774 5.0623
a=15.4985 3.9261 4.1044 4.2959 4.4915 4.7478 5.0214 5.3017
a=16.2028 4.1681 4.3421 4.5309 4.7266 4.987 5.2603 5.5374
a=16.9282 4.4061 4.5758 4.7615 4.957 5.221 5.4941 5.7688
a=17.675 4.6406 4.8061 4.9877 5.1828 5.4497 5.7226 5.9957
a=18.4433 4.8714 5.0324 5.2094 5.4048 5.6731 5.946 6.2178
a=19.2337 5.0986 5.2551 5.4269 5.6226 5.8914 6.1642 6.435
a=20.0462 5.3223 5.4745 5.6421 5.8361 6.1047 6.3773 6.6473
a=20.8814 5.5425 5.6905 5.8537 6.0451 6.3131 6.5854 6.8547
a=21.7394 5.7593 5.9032 6.0624 6.2498 6.5167 6.7884 7.0571
a=22.6207 5.9728 6.1127 6.2676 6.4503 6.7156 6.9866 7.2547
a=23.5254 6.1831 6.319 6.4696 6.65 6.9098 7.18 7.4474
a=24.4539 6.3901 6.5222 6.6701 6.8461 7.1 7.3688 7.6354
a=25.4065 6.594 6.7224 6.8687 7.0407 7.2862 7.5532 7.8188
a=26.3836 6.7948 6.9218 7.0645 7.2318 7.4687 7.7332 7.9978
a=27.3854 6.9935 7.1185 7.2579 7.4209 7.6481 7.9092 8.1724
a=28.4123 7.1908 7.3127 7.4484 7.6068 7.8268 8.0813 8.3428
a=29.4645 7.3853 7.5041 7.6365 7.7905 8.0026 8.2501 8.5092
a=30.5423 7.5773 7.6931 7.8218 7.9714 8.1771 8.4157 8.672
a=31.6462 7.7667 7.8795 8.0048 8.1501 8.3489 8.5785 8.8313
a=32.7763 7.9535 8.0633 8.1851 8.3261 8.519 8.7392 8.9873
a=33.933 8.1377 8.2446 8.3631 8.5 8.6866 8.8991 9.1403
a=35.1166 8.3194 8.4233 8.5385 8.6714 8.8523 9.0569 9.2909
a=36.3274 8.4986 8.5997 8.7117 8.8405 9.0157 9.2136 9.439
a=37.5657 8.6753 8.7736 8.8823 9.0074 9.1772 9.3682 9.5851
a=38.8319 8.8495 8.9451 9.0507 9.1721 9.3365 9.5215 9.7291
a=40.1262 9.0213 9.1141 9.2167 9.3345 9.494 9.6728 9.8727
a=41.4489 9.1907 9.2809 9.3806 9.4948 9.6494 9.8226 10.015
a=42.8004 9.3577 9.4453 9.5421 9.6529 9.8028 9.9705 10.156
a=44.181 9.5224 9.6075 9.7014 9.8089 9.9543 10.117 10.296
a=45.5909 9.6847 9.7673 9.8585 9.9629 10.104 10.261 10.435
a=47.0305 9.8447 9.925 10.014 10.115 10.252 10.404 10.572
a=48.5001 10.003 10.08 10.166 10.265 10.398 10.546 10.708
a=50 10.158 10.234 10.317 10.413 10.542 10.685 10.842
Table for Policy function at l=0
z1 z2 z3 z4 z5 z6 z7
__________ __________ __________ __________ __________ _________ _________
a=0 0.00041224 0.00041224 0.00041224 0.00041224 0.00041224 0.0013913 0.0013913
a=5.15305e-05 0.00041224 0.00041224 0.00041224 0.00041224 0.0013913 0.0013913 0.0013913
a=0.000412244 0.0013913 0.0013913 0.0013913 0.0013913 0.0013913 0.0013913 0.0013913
a=0.00139132 0.003298 0.003298 0.003298 0.003298 0.003298 0.003298 0.003298
a=0.00329795 0.0064413 0.0064413 0.0064413 0.0064413 0.0064413 0.0064413 0.0064413
a=0.00644131 0.011131 0.011131 0.011131 0.011131 0.011131 0.011131 0.011131
a=0.0111306 0.011131 0.017675 0.017675 0.017675 0.017675 0.017675 0.017675
a=0.017675 0.017675 0.017675 0.017675 0.026384 0.026384 0.026384 0.026384
a=0.0263836 0.026384 0.026384 0.026384 0.026384 0.037566 0.037566 0.037566
a=0.0375657 0.037566 0.037566 0.037566 0.037566 0.037566 0.051531 0.051531
a=0.0515305 0.051531 0.051531 0.051531 0.051531 0.051531 0.051531 0.068587
a=0.0685871 0.068587 0.068587 0.068587 0.068587 0.068587 0.068587 0.089045
a=0.0890447 0.089045 0.089045 0.089045 0.089045 0.089045 0.089045 0.089045
a=0.113213 0.11321 0.11321 0.11321 0.11321 0.11321 0.11321 0.11321
a=0.1414 0.1414 0.1414 0.1414 0.1414 0.1414 0.1414 0.1414
a=0.173915 0.1414 0.17392 0.17392 0.17392 0.17392 0.17392 0.17392
a=0.211069 0.17392 0.17392 0.21107 0.21107 0.21107 0.21107 0.21107
a=0.253169 0.21107 0.21107 0.25317 0.25317 0.25317 0.25317 0.25317
a=0.300526 0.25317 0.25317 0.25317 0.30053 0.30053 0.30053 0.30053
a=0.353448 0.30053 0.30053 0.30053 0.35345 0.35345 0.35345 0.35345
a=0.412244 0.35345 0.35345 0.35345 0.41224 0.41224 0.41224 0.41224
a=0.477224 0.41224 0.41224 0.41224 0.41224 0.47722 0.47722 0.47722
a=0.548697 0.47722 0.47722 0.47722 0.5487 0.5487 0.5487 0.5487
a=0.626972 0.5487 0.5487 0.5487 0.5487 0.62697 0.62697 0.62697
a=0.712358 0.62697 0.62697 0.62697 0.71236 0.71236 0.71236 0.71236
a=0.805164 0.71236 0.71236 0.71236 0.71236 0.80516 0.80516 0.80516
a=0.9057 0.80516 0.80516 0.80516 0.80516 0.9057 0.9057 0.9057
a=1.01427 0.9057 0.9057 0.9057 0.9057 1.0143 1.0143 1.0143
a=1.1312 1.0143 1.0143 1.0143 1.0143 1.1312 1.1312 1.1312
a=1.25678 1.1312 1.1312 1.1312 1.1312 1.2568 1.2568 1.2568
a=1.39132 1.2568 1.2568 1.2568 1.2568 1.3913 1.3913 1.3913
a=1.53515 1.3913 1.3913 1.3913 1.3913 1.5351 1.5351 1.5351
a=1.68855 1.5351 1.5351 1.5351 1.5351 1.6886 1.6886 1.6886
a=1.85185 1.6886 1.6886 1.6886 1.6886 1.8519 1.8519 1.8519
a=2.02536 1.8519 1.8519 1.8519 1.8519 2.0254 2.0254 2.0254
a=2.20937 2.0254 2.0254 2.0254 2.0254 2.2094 2.2094 2.2094
a=2.40421 2.2094 2.2094 2.2094 2.2094 2.4042 2.4042 2.4042
a=2.61017 2.4042 2.4042 2.4042 2.4042 2.6102 2.6102 2.6102
a=2.82758 2.6102 2.6102 2.6102 2.6102 2.8276 2.8276 2.8276
a=3.05674 2.8276 2.8276 2.8276 2.8276 3.0567 3.0567 3.0567
a=3.29795 3.0567 3.0567 3.0567 3.0567 3.298 3.298 3.298
a=3.55153 3.298 3.298 3.298 3.298 3.5515 3.5515 3.5515
a=3.81779 3.5515 3.5515 3.5515 3.5515 3.8178 3.8178 3.8178
a=4.09704 3.8178 3.8178 3.8178 3.8178 4.097 4.097 4.097
a=4.38957 4.097 4.097 4.097 4.097 4.3896 4.3896 4.3896
a=4.69572 4.3896 4.3896 4.3896 4.3896 4.6957 4.6957 4.6957
a=5.01577 4.6957 4.6957 4.6957 4.6957 5.0158 5.0158 5.0158
a=5.35005 5.0158 5.0158 5.0158 5.0158 5.0158 5.3501 5.3501
a=5.69886 5.3501 5.3501 5.3501 5.3501 5.3501 5.6989 5.6989
a=6.06251 5.6989 5.6989 5.6989 5.6989 5.6989 6.0625 6.0625
a=6.44131 6.0625 6.0625 6.0625 6.0625 6.0625 6.4413 6.4413
a=6.83557 6.4413 6.4413 6.4413 6.4413 6.4413 6.8356 6.8356
a=7.2456 6.8356 6.8356 6.8356 6.8356 6.8356 7.2456 7.2456
a=7.67171 7.2456 7.2456 7.2456 7.2456 7.2456 7.6717 7.6717
a=8.1142 7.6717 7.6717 7.6717 7.6717 7.6717 8.1142 8.1142
a=8.57339 8.1142 8.1142 8.1142 8.1142 8.1142 8.5734 8.5734
a=9.04958 8.5734 8.5734 8.5734 8.5734 8.5734 9.0496 9.0496
a=9.54309 9.0496 9.0496 9.0496 9.0496 9.0496 9.0496 9.5431
a=10.0542 9.5431 9.5431 9.5431 9.5431 9.5431 9.5431 10.054
a=10.5833 10.054 10.054 10.054 10.054 10.054 10.054 10.583
a=11.1306 10.583 10.583 10.583 10.583 10.583 10.583 11.131
a=11.6964 10.583 11.131 11.131 11.131 11.131 11.131 11.696
a=12.2812 11.131 11.696 11.696 11.696 11.696 11.696 12.281
a=12.885 11.696 12.281 12.281 12.281 12.281 12.281 12.885
a=13.5084 12.281 12.885 12.885 12.885 12.885 12.885 13.508
a=14.1516 12.885 13.508 13.508 13.508 13.508 13.508 13.508
a=14.8148 13.508 13.508 14.152 14.152 14.152 14.152 14.152
a=15.4985 14.152 14.152 14.815 14.815 14.815 14.815 14.815
a=16.2028 14.815 14.815 15.498 15.498 15.498 15.498 15.498
a=16.9282 15.498 15.498 16.203 16.203 16.203 16.203 16.203
a=17.675 16.203 16.203 16.928 16.928 16.928 16.928 16.928
a=18.4433 16.928 16.928 16.928 17.675 17.675 17.675 17.675
a=19.2337 17.675 17.675 17.675 18.443 18.443 18.443 18.443
a=20.0462 18.443 18.443 18.443 19.234 19.234 19.234 19.234
a=20.8814 19.234 19.234 19.234 20.046 20.046 20.046 20.046
a=21.7394 20.046 20.046 20.046 20.881 20.881 20.881 20.881
a=22.6207 20.881 20.881 20.881 21.739 21.739 21.739 21.739
a=23.5254 21.739 21.739 21.739 21.739 22.621 22.621 22.621
a=24.4539 22.621 22.621 22.621 22.621 23.525 23.525 23.525
a=25.4065 23.525 23.525 23.525 23.525 24.454 24.454 24.454
a=26.3836 24.454 24.454 24.454 24.454 25.407 25.407 25.407
a=27.3854 25.407 25.407 25.407 25.407 25.407 26.384 26.384
a=28.4123 26.384 26.384 26.384 26.384 26.384 27.385 27.385
a=29.4645 27.385 27.385 27.385 27.385 27.385 28.412 28.412
a=30.5423 28.412 28.412 28.412 28.412 28.412 29.464 29.464
a=31.6462 29.464 29.464 29.464 29.464 29.464 30.542 30.542
a=32.7763 30.542 30.542 30.542 30.542 30.542 30.542 31.646
a=33.933 31.646 31.646 31.646 31.646 31.646 31.646 32.776
a=35.1166 32.776 32.776 32.776 32.776 32.776 32.776 33.933
a=36.3274 33.933 33.933 33.933 33.933 33.933 33.933 35.117
a=37.5657 35.117 35.117 35.117 35.117 35.117 35.117 36.327
a=38.8319 36.327 36.327 36.327 36.327 36.327 36.327 37.566
a=40.1262 37.566 37.566 37.566 37.566 37.566 37.566 37.566
a=41.4489 38.832 38.832 38.832 38.832 38.832 38.832 38.832
a=42.8004 40.126 40.126 40.126 40.126 40.126 40.126 40.126
a=44.181 41.449 41.449 41.449 41.449 41.449 41.449 41.449
a=45.5909 42.8 42.8 42.8 42.8 42.8 42.8 42.8
a=47.0305 44.181 44.181 44.181 44.181 44.181 44.181 44.181
a=48.5001 45.591 45.591 45.591 45.591 45.591 45.591 45.591
a=50 47.031 47.031 47.031 47.031 47.031 47.031 47.031
Table for Policy function at l=1
z1 z2 z3 z4 z5 z6 z7
_______ _______ _______ _______ _______ _______ _______
a=0 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236
a=5.15305e-05 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236
a=0.000412244 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236
a=0.00139132 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236
a=0.00329795 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236
a=0.00644131 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236
a=0.0111306 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236
a=0.017675 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236
a=0.0263836 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236
a=0.0375657 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236
a=0.0515305 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236
a=0.0685871 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236
a=0.0890447 0.80516 0.80516 0.71236 0.71236 0.71236 0.71236 0.71236
a=0.113213 0.80516 0.80516 0.80516 0.80516 0.80516 0.80516 0.80516
a=0.1414 0.80516 0.80516 0.80516 0.80516 0.80516 0.80516 0.80516
a=0.173915 0.80516 0.80516 0.80516 0.80516 0.80516 0.80516 0.80516
a=0.211069 0.9057 0.9057 0.80516 0.80516 0.80516 0.80516 0.80516
a=0.253169 0.9057 0.9057 0.9057 0.9057 0.9057 0.9057 0.9057
a=0.300526 0.9057 0.9057 0.9057 0.9057 0.9057 0.9057 0.9057
a=0.353448 1.0143 1.0143 1.0143 1.0143 1.0143 1.0143 0.9057
a=0.412244 1.0143 1.0143 1.0143 1.0143 1.0143 1.0143 1.0143
a=0.477224 1.1312 1.1312 1.1312 1.0143 1.0143 1.0143 1.0143
a=0.548697 1.1312 1.1312 1.1312 1.1312 1.1312 1.1312 1.1312
a=0.626972 1.2568 1.2568 1.2568 1.2568 1.2568 1.2568 1.1312
a=0.712358 1.2568 1.2568 1.2568 1.2568 1.2568 1.2568 1.2568
a=0.805164 1.3913 1.3913 1.3913 1.3913 1.3913 1.3913 1.3913
a=0.9057 1.3913 1.3913 1.3913 1.3913 1.3913 1.3913 1.3913
a=1.01427 1.5351 1.5351 1.5351 1.5351 1.5351 1.5351 1.5351
a=1.1312 1.6886 1.6886 1.6886 1.6886 1.6886 1.6886 1.5351
a=1.25678 1.6886 1.6886 1.6886 1.6886 1.6886 1.6886 1.6886
a=1.39132 1.8519 1.8519 1.8519 1.8519 1.8519 1.8519 1.8519
a=1.53515 2.0254 2.0254 2.0254 2.0254 2.0254 2.0254 2.0254
a=1.68855 2.2094 2.2094 2.0254 2.0254 2.0254 2.0254 2.0254
a=1.85185 2.2094 2.2094 2.2094 2.2094 2.2094 2.2094 2.2094
a=2.02536 2.4042 2.4042 2.4042 2.4042 2.4042 2.4042 2.4042
a=2.20937 2.6102 2.6102 2.6102 2.6102 2.6102 2.6102 2.6102
a=2.40421 2.8276 2.8276 2.8276 2.8276 2.8276 2.8276 2.8276
a=2.61017 2.8276 2.8276 2.8276 2.8276 2.8276 2.8276 2.8276
a=2.82758 3.0567 3.0567 3.0567 3.0567 3.0567 3.0567 3.0567
a=3.05674 3.298 3.298 3.298 3.298 3.298 3.298 3.298
a=3.29795 3.5515 3.5515 3.5515 3.5515 3.5515 3.5515 3.5515
a=3.55153 3.8178 3.8178 3.8178 3.8178 3.8178 3.8178 3.8178
a=3.81779 4.097 4.097 4.097 4.097 4.097 4.097 4.097
a=4.09704 4.3896 4.3896 4.3896 4.3896 4.3896 4.3896 4.3896
a=4.38957 4.3896 4.3896 4.3896 4.3896 4.3896 4.3896 4.3896
a=4.69572 4.6957 4.6957 4.6957 4.6957 4.6957 4.6957 4.6957
a=5.01577 5.0158 5.0158 5.0158 5.0158 5.0158 5.0158 5.0158
a=5.35005 5.3501 5.3501 5.3501 5.3501 5.3501 5.3501 5.3501
a=5.69886 5.6989 5.6989 5.6989 5.6989 5.6989 5.6989 5.6989
a=6.06251 6.0625 6.0625 6.0625 6.0625 6.0625 6.0625 6.0625
a=6.44131 6.4413 6.4413 6.4413 6.4413 6.4413 6.4413 6.4413
a=6.83557 6.8356 6.8356 6.8356 6.8356 6.8356 6.8356 6.8356
a=7.2456 7.2456 7.2456 7.2456 7.2456 7.2456 7.2456 7.2456
a=7.67171 7.6717 7.6717 7.6717 7.6717 7.6717 7.6717 7.6717
a=8.1142 8.1142 8.1142 8.1142 8.1142 7.6717 8.1142 8.1142
a=8.57339 8.1142 8.1142 8.1142 8.1142 8.1142 8.5734 8.5734
a=9.04958 8.5734 8.5734 8.5734 8.5734 8.5734 9.0496 9.0496
a=9.54309 9.0496 9.0496 9.0496 9.0496 9.0496 9.0496 9.5431
a=10.0542 9.5431 9.5431 9.5431 9.5431 9.5431 9.5431 10.054
a=10.5833 10.054 10.054 10.054 10.054 10.054 10.054 10.583
a=11.1306 10.583 10.583 10.583 10.583 10.583 10.583 11.131
a=11.6964 11.131 11.131 11.131 11.131 11.131 11.131 11.696
a=12.2812 11.696 11.696 11.696 11.696 11.696 11.696 12.281
a=12.885 12.281 12.281 12.281 12.281 12.281 12.281 12.885
a=13.5084 12.885 12.885 12.885 12.885 12.885 12.885 13.508
a=14.1516 13.508 13.508 13.508 13.508 13.508 13.508 13.508
a=14.8148 14.152 14.152 14.152 14.152 14.152 14.152 14.152
a=15.4985 14.815 14.815 14.815 14.815 14.815 14.815 14.815
a=16.2028 15.498 15.498 15.498 15.498 15.498 15.498 15.498
a=16.9282 16.203 16.203 16.203 16.203 16.203 16.203 16.203
a=17.675 16.928 16.928 16.928 16.928 16.928 16.928 16.928
a=18.4433 17.675 17.675 17.675 17.675 17.675 17.675 17.675
a=19.2337 18.443 18.443 18.443 18.443 18.443 18.443 18.443
a=20.0462 19.234 19.234 19.234 19.234 19.234 19.234 19.234
a=20.8814 20.046 20.046 20.046 20.046 20.046 20.046 20.046
a=21.7394 20.881 20.881 20.881 20.881 20.881 20.881 20.881
a=22.6207 21.739 21.739 21.739 21.739 21.739 21.739 21.739
a=23.5254 22.621 22.621 22.621 21.739 22.621 22.621 22.621
a=24.4539 23.525 23.525 22.621 22.621 23.525 23.525 23.525
a=25.4065 24.454 24.454 23.525 23.525 24.454 24.454 24.454
a=26.3836 25.407 24.454 24.454 24.454 25.407 25.407 25.407
a=27.3854 25.407 25.407 25.407 25.407 25.407 26.384 26.384
a=28.4123 26.384 26.384 26.384 26.384 26.384 27.385 27.385
a=29.4645 27.385 27.385 27.385 27.385 27.385 28.412 28.412
a=30.5423 28.412 28.412 28.412 28.412 28.412 29.464 29.464
a=31.6462 29.464 29.464 29.464 29.464 29.464 30.542 30.542
a=32.7763 30.542 30.542 30.542 30.542 30.542 30.542 31.646
a=33.933 31.646 31.646 31.646 31.646 31.646 31.646 32.776
a=35.1166 32.776 32.776 32.776 32.776 32.776 32.776 33.933
a=36.3274 33.933 33.933 33.933 33.933 33.933 33.933 35.117
a=37.5657 35.117 35.117 35.117 35.117 35.117 35.117 36.327
a=38.8319 36.327 36.327 36.327 36.327 36.327 36.327 37.566
a=40.1262 37.566 37.566 37.566 37.566 37.566 37.566 37.566
a=41.4489 38.832 38.832 38.832 38.832 38.832 38.832 38.832
a=42.8004 40.126 40.126 40.126 40.126 40.126 40.126 40.126
a=44.181 41.449 41.449 41.449 41.449 41.449 41.449 41.449
a=45.5909 42.8 42.8 42.8 42.8 42.8 42.8 42.8
a=47.0305 44.181 44.181 44.181 44.181 44.181 44.181 44.181
a=48.5001 45.591 45.591 45.591 45.591 45.591 45.591 45.591
a=50 47.031 47.031 47.031 47.031 47.031 47.031 47.031
Table for Consumption at l=0
z1 z2 z3 z4 z5 z6 z7
_________ _________ _________ _________ _________ _________ _________
a=0 0.0011724 0.0014819 0.0018518 0.002294 0.0028225 0.0024751 0.0032303
a=5.15305e-05 0.0016336 0.0019526 0.0023312 0.0027812 0.0023377 0.0029758 0.0037364
a=0.000412244 0.002688 0.0031071 0.0035878 0.0041413 0.0047809 0.0055226 0.0063857
a=0.00139132 0.0044279 0.005037 0.0057197 0.0064872 0.0073528 0.0083324 0.0094448
a=0.00329795 0.0067154 0.0075809 0.0085382 0.0095997 0.01078 0.012095 0.013566
a=0.00644131 0.0094984 0.010679 0.011975 0.0134 0.014969 0.016702 0.018619
a=0.0111306 0.019287 0.014296 0.015991 0.017844 0.019872 0.022097 0.024543
a=0.017675 0.02513 0.027111 0.029264 0.022899 0.025454 0.028245 0.031296
a=0.0263836 0.031689 0.034152 0.036822 0.03972 0.031686 0.035112 0.038846
a=0.0375657 0.038941 0.041941 0.045185 0.048698 0.052505 0.042672 0.047162
a=0.0515305 0.046861 0.050452 0.054328 0.058517 0.063048 0.067956 0.05622
a=0.0685871 0.055426 0.059662 0.064227 0.069153 0.074474 0.080225 0.065992
a=0.0890447 0.064614 0.069547 0.074859 0.080583 0.086757 0.093421 0.10062
a=0.113213 0.0744 0.080085 0.0862 0.092783 0.099874 0.10752 0.11577
a=0.1414 0.084763 0.091253 0.098228 0.10573 0.1138 0.1225 0.13187
a=0.173915 0.1282 0.10303 0.11092 0.1194 0.12852 0.13833 0.1489
a=0.211069 0.14428 0.15254 0.12425 0.13377 0.144 0.155 0.16683
a=0.253169 0.16119 0.17041 0.13821 0.14882 0.16022 0.17247 0.18565
a=0.300526 0.17889 0.18913 0.20011 0.16453 0.17717 0.19074 0.20532
a=0.353448 0.19737 0.20868 0.2208 0.18087 0.19481 0.20977 0.22583
a=0.412244 0.21661 0.22904 0.24235 0.19783 0.21313 0.22954 0.24716
a=0.477224 0.23657 0.25017 0.26474 0.28035 0.2321 0.25003 0.26928
a=0.548697 0.25724 0.27208 0.28795 0.23348 0.2517 0.27123 0.29217
a=0.626972 0.2786 0.29472 0.31196 0.33042 0.27191 0.29309 0.3158
a=0.712358 0.30063 0.31808 0.33674 0.27133 0.29271 0.31562 0.34016
a=0.805164 0.32331 0.34214 0.36227 0.38382 0.31408 0.33877 0.36523
a=0.9057 0.34661 0.36687 0.38854 0.41172 0.33599 0.36254 0.39097
a=1.01427 0.37051 0.39226 0.41551 0.44038 0.35842 0.38689 0.41738
a=1.1312 0.39499 0.41828 0.44318 0.4698 0.38135 0.41181 0.44442
a=1.25678 0.42003 0.44491 0.47151 0.49994 0.40476 0.43728 0.47208
a=1.39132 0.44561 0.47214 0.50048 0.53078 0.42863 0.46327 0.50034
a=1.53515 0.4717 0.49993 0.53008 0.56231 0.45293 0.48976 0.52916
a=1.68855 0.49829 0.52827 0.56029 0.59449 0.47765 0.51674 0.55854
a=1.85185 0.52536 0.55713 0.59107 0.62732 0.50276 0.54417 0.58845
a=2.02536 0.55288 0.5865 0.62241 0.66077 0.52824 0.57204 0.61886
a=2.20937 0.58083 0.61636 0.65429 0.69481 0.55408 0.60033 0.64976
a=2.40421 0.60919 0.64668 0.68669 0.72943 0.58024 0.62901 0.68113
a=2.61017 0.63794 0.67743 0.71959 0.7646 0.6067 0.65806 0.71294
a=2.82758 0.66706 0.70861 0.75296 0.80031 0.63346 0.68747 0.74517
a=3.05674 0.69653 0.74019 0.78678 0.83652 0.66048 0.7172 0.7778
a=3.29795 0.72632 0.77214 0.82104 0.87323 0.68774 0.74725 0.81081
a=3.55153 0.75642 0.80446 0.85571 0.91041 0.71523 0.77758 0.84418
a=3.81779 0.7868 0.8371 0.89077 0.94804 0.74291 0.80818 0.87789
a=4.09704 0.81745 0.87006 0.9262 0.98609 0.77078 0.83903 0.91191
a=4.38957 0.84834 0.90332 0.96197 1.0246 0.7988 0.8701 0.94622
a=4.69572 0.87945 0.93685 0.99807 1.0634 0.82696 0.90137 0.98081
a=5.01577 0.91076 0.97063 1.0345 1.1026 0.85524 0.93282 1.0156
a=5.35005 0.94225 1.0046 1.0712 1.1422 1.2179 0.96444 1.0507
a=5.69886 0.9739 1.0389 1.1081 1.182 1.2609 0.99619 1.086
a=6.06251 1.0057 1.0733 1.1453 1.2222 1.3042 1.0281 1.1215
a=6.44131 1.0376 1.1078 1.1828 1.2627 1.3479 1.06 1.1571
a=6.83557 1.0696 1.1426 1.2204 1.3034 1.3919 1.0921 1.1929
a=7.2456 1.1016 1.1774 1.2582 1.3443 1.4362 1.1242 1.2288
a=7.67171 1.1338 1.2124 1.2961 1.3855 1.4808 1.1563 1.2648
a=8.1142 1.1659 1.2474 1.3342 1.4268 1.5256 1.1885 1.3009
a=8.57339 1.1981 1.2825 1.3725 1.4684 1.5707 1.2206 1.337
a=9.04958 1.2303 1.3176 1.4108 1.5101 1.616 1.2527 1.3732
a=9.54309 1.2624 1.3528 1.4492 1.5519 1.6614 1.7783 1.4094
a=10.0542 1.2945 1.3879 1.4876 1.5938 1.7071 1.8279 1.4456
a=10.5833 1.3265 1.4231 1.5261 1.6358 1.7529 1.8777 1.4818
a=11.1306 1.3584 1.4582 1.5646 1.6779 1.7988 1.9277 1.5179
a=11.6964 1.9375 1.4932 1.603 1.7201 1.8449 1.9779 1.554
a=12.2812 1.9877 1.5282 1.6415 1.7622 1.891 2.0283 1.59
a=12.885 2.0381 1.563 1.6799 1.8044 1.9372 2.0788 1.6258
a=13.5084 2.0892 1.5977 1.7182 1.8466 1.9834 2.1294 1.6616
a=14.1516 2.1413 1.6323 1.7564 1.8887 2.0297 2.18 2.3403
a=14.8148 2.1945 2.3099 1.7945 1.9308 2.076 2.2308 2.3959
a=15.4985 2.2487 2.3642 1.8325 1.9728 2.1223 2.2816 2.4515
a=16.2028 2.3039 2.4195 1.8704 2.0147 2.1685 2.3324 2.5072
a=16.9282 2.3601 2.4757 1.908 2.0565 2.2147 2.3833 2.5631
a=17.675 2.4174 2.533 1.9455 2.0981 2.2607 2.4341 2.6189
a=18.4433 2.4758 2.5913 2.7295 2.1396 2.3067 2.4849 2.6748
a=19.2337 2.5352 2.6507 2.7889 2.1809 2.3526 2.5356 2.7307
a=20.0462 2.5956 2.7112 2.8494 2.222 2.3983 2.5863 2.7866
a=20.8814 2.6572 2.7727 2.9109 2.2628 2.4439 2.6368 2.8425
a=21.7394 2.7198 2.8353 2.9735 2.3035 2.4893 2.6873 2.8983
a=22.6207 2.7835 2.899 3.0372 2.3443 2.5345 2.7376 2.9541
a=23.5254 2.8482 2.9638 3.102 3.2671 2.5794 2.7877 3.0098
a=24.4539 2.9141 3.0297 3.1678 3.333 2.6241 2.8377 3.0653
a=25.4065 2.9811 3.0967 3.2348 3.4 2.6686 2.8875 3.1207
a=26.3836 3.0492 3.1648 3.3029 3.468 2.7128 2.937 3.176
a=27.3854 3.1184 3.234 3.3721 3.5372 3.7346 2.9863 3.2311
a=28.4123 3.1887 3.3043 3.4424 3.6075 3.8049 3.0353 3.286
a=29.4645 3.2601 3.3757 3.5139 3.679 3.8764 3.0841 3.3407
a=30.5423 3.3327 3.4483 3.5864 3.7516 3.9489 3.1325 3.3951
a=31.6462 3.4064 3.522 3.6601 3.8253 4.0227 3.1807 3.4493
a=32.7763 3.4813 3.5968 3.735 3.9001 4.0975 4.3334 3.5032
a=33.933 3.5573 3.6729 3.811 3.9761 4.1735 4.4094 3.5568
a=35.1166 3.6344 3.75 3.8882 4.0533 4.2507 4.4866 3.6101
a=36.3274 3.7128 3.8283 3.9665 4.1316 4.329 4.5649 3.663
a=37.5657 3.7922 3.9078 4.046 4.2111 4.4085 4.6444 3.7156
a=38.8319 3.8729 3.9885 4.1266 4.2918 4.4892 4.7251 3.7688
a=40.1262 3.9548 4.0703 4.2085 4.3736 4.571 4.8069 5.0889
a=41.4489 4.0378 4.1534 4.2915 4.4567 4.654 4.89 5.172
a=42.8004 4.122 4.2376 4.3758 4.5409 4.7383 4.9742 5.2562
a=44.181 4.2075 4.323 4.4612 4.6263 4.8237 5.0596 5.3416
a=45.5909 4.2941 4.4097 4.5478 4.7129 4.9103 5.1463 5.4283
a=47.0305 4.3819 4.4975 4.6357 4.8008 4.9982 5.2341 5.5161
a=48.5001 4.471 4.5866 4.7247 4.8899 5.0872 5.3232 5.6052
a=50 4.5613 4.6769 4.815 4.9801 5.1775 5.4135 5.6955
Table for Consumption at l=1
z1 z2 z3 z4 z5 z6 z7
_______ _______ _______ _______ _______ _______ _______
a=0 0.28764 0.28764 0.28764 0.28764 0.28764 0.28764 0.28764
a=5.15305e-05 0.28769 0.28769 0.28769 0.28769 0.28769 0.28769 0.28769
a=0.000412244 0.28806 0.28806 0.28806 0.28806 0.28806 0.28806 0.28806
a=0.00139132 0.28906 0.28906 0.28906 0.28906 0.28906 0.28906 0.28906
a=0.00329795 0.29101 0.29101 0.29101 0.29101 0.29101 0.29101 0.29101
a=0.00644131 0.29421 0.29421 0.29421 0.29421 0.29421 0.29421 0.29421
a=0.0111306 0.299 0.299 0.299 0.299 0.299 0.299 0.299
a=0.017675 0.30567 0.30567 0.30567 0.30567 0.30567 0.30567 0.30567
a=0.0263836 0.31455 0.31455 0.31455 0.31455 0.31455 0.31455 0.31455
a=0.0375657 0.32596 0.32596 0.32596 0.32596 0.32596 0.32596 0.32596
a=0.0515305 0.3402 0.3402 0.3402 0.3402 0.3402 0.3402 0.3402
a=0.0685871 0.3576 0.3576 0.3576 0.3576 0.3576 0.3576 0.3576
a=0.0890447 0.28566 0.28566 0.37847 0.37847 0.37847 0.37847 0.37847
a=0.113213 0.31031 0.31031 0.31031 0.31031 0.31031 0.31031 0.31031
a=0.1414 0.33906 0.33906 0.33906 0.33906 0.33906 0.33906 0.33906
a=0.173915 0.37223 0.37223 0.37223 0.37223 0.37223 0.37223 0.37223
a=0.211069 0.30959 0.30959 0.41013 0.41013 0.41013 0.41013 0.41013
a=0.253169 0.35253 0.35253 0.35253 0.35253 0.35253 0.35253 0.35253
a=0.300526 0.40084 0.40084 0.40084 0.40084 0.40084 0.40084 0.40084
a=0.353448 0.34624 0.34624 0.34624 0.34624 0.34624 0.34624 0.45482
a=0.412244 0.40621 0.40621 0.40621 0.40621 0.40621 0.40621 0.40621
a=0.477224 0.35557 0.35557 0.35557 0.47249 0.47249 0.47249 0.47249
a=0.548697 0.42847 0.42847 0.42847 0.42847 0.42847 0.42847 0.42847
a=0.626972 0.38273 0.38273 0.38273 0.38273 0.38273 0.38273 0.50831
a=0.712358 0.46983 0.46983 0.46983 0.46983 0.46983 0.46983 0.46983
a=0.805164 0.42994 0.42994 0.42994 0.42994 0.42994 0.42994 0.42994
a=0.9057 0.53249 0.53249 0.53249 0.53249 0.53249 0.53249 0.53249
a=1.01427 0.49942 0.49942 0.49942 0.49942 0.49942 0.49942 0.49942
a=1.1312 0.46527 0.46527 0.46527 0.46527 0.46527 0.46527 0.61868
a=1.25678 0.59336 0.59336 0.59336 0.59336 0.59336 0.59336 0.59336
a=1.39132 0.5673 0.5673 0.5673 0.5673 0.5673 0.5673 0.5673
a=1.53515 0.54049 0.54049 0.54049 0.54049 0.54049 0.54049 0.54049
a=1.68855 0.51295 0.51295 0.69697 0.69697 0.69697 0.69697 0.69697
a=1.85185 0.67952 0.67952 0.67952 0.67952 0.67952 0.67952 0.67952
a=2.02536 0.66165 0.66165 0.66165 0.66165 0.66165 0.66165 0.66165
a=2.20937 0.64338 0.64338 0.64338 0.64338 0.64338 0.64338 0.64338
a=2.40421 0.62471 0.62471 0.62471 0.62471 0.62471 0.62471 0.62471
a=2.61017 0.8348 0.8348 0.8348 0.8348 0.8348 0.8348 0.8348
a=2.82758 0.8274 0.8274 0.8274 0.8274 0.8274 0.8274 0.8274
a=3.05674 0.81992 0.81992 0.81992 0.81992 0.81992 0.81992 0.81992
a=3.29795 0.81238 0.81238 0.81238 0.81238 0.81238 0.81238 0.81238
a=3.55153 0.80477 0.80477 0.80477 0.80477 0.80477 0.80477 0.80477
a=3.81779 0.79711 0.79711 0.79711 0.79711 0.79711 0.79711 0.79711
a=4.09704 0.7894 0.7894 0.7894 0.7894 0.7894 0.7894 0.7894
a=4.38957 1.0878 1.0878 1.0878 1.0878 1.0878 1.0878 1.0878
a=4.69572 1.0939 1.0939 1.0939 1.0939 1.0939 1.0939 1.0939
a=5.01577 1.1003 1.1003 1.1003 1.1003 1.1003 1.1003 1.1003
a=5.35005 1.107 1.107 1.107 1.107 1.107 1.107 1.107
a=5.69886 1.114 1.114 1.114 1.114 1.114 1.114 1.114
a=6.06251 1.1213 1.1213 1.1213 1.1213 1.1213 1.1213 1.1215
a=6.44131 1.1288 1.1288 1.1288 1.1288 1.1288 1.1288 1.1571
a=6.83557 1.1367 1.1367 1.1367 1.1367 1.1367 1.1367 1.1929
a=7.2456 1.1449 1.1449 1.1449 1.1449 1.1449 1.1449 1.2288
a=7.67171 1.1534 1.1534 1.1534 1.1534 1.1534 1.1563 1.2648
a=8.1142 1.1623 1.1623 1.1623 1.1623 1.6048 1.1885 1.3009
a=8.57339 1.6307 1.6307 1.6307 1.6307 1.6307 1.2206 1.337
a=9.04958 1.6572 1.6572 1.6572 1.6572 1.6572 1.2527 1.3732
a=9.54309 1.6844 1.6844 1.6844 1.6844 1.6844 1.7783 1.4094
a=10.0542 1.7122 1.7122 1.7122 1.7122 1.7122 1.8279 1.4456
a=10.5833 1.7407 1.7407 1.7407 1.7407 1.7529 1.8777 1.4818
a=11.1306 1.7699 1.7699 1.7699 1.7699 1.7988 1.9277 1.5179
a=11.6964 1.7998 1.7998 1.7998 1.7998 1.8449 1.9779 1.554
a=12.2812 1.8303 1.8303 1.8303 1.8303 1.891 2.0283 1.59
a=12.885 1.8616 1.8616 1.8616 1.8616 1.9372 2.0788 1.6258
a=13.5084 1.8935 1.8935 1.8935 1.8935 1.9834 2.1294 1.6616
a=14.1516 1.9262 1.9262 1.9262 1.9262 2.0297 2.18 2.3403
a=14.8148 1.9595 1.9595 1.9595 1.9595 2.076 2.2308 2.3959
a=15.4985 1.9936 1.9936 1.9936 1.9936 2.1223 2.2816 2.4515
a=16.2028 2.0284 2.0284 2.0284 2.0284 2.1685 2.3324 2.5072
a=16.9282 2.064 2.064 2.064 2.064 2.2147 2.3833 2.5631
a=17.675 2.1002 2.1002 2.1002 2.1002 2.2607 2.4341 2.6189
a=18.4433 2.1372 2.1372 2.1372 2.1396 2.3067 2.4849 2.6748
a=19.2337 2.175 2.175 2.175 2.1809 2.3526 2.5356 2.7307
a=20.0462 2.2135 2.2135 2.2135 2.222 2.3983 2.5863 2.7866
a=20.8814 2.2528 2.2528 2.2528 2.2628 2.4439 2.6368 2.8425
a=21.7394 2.2928 2.2928 2.2928 2.3035 2.4893 2.6873 2.8983
a=22.6207 2.3336 2.3336 2.3336 2.3443 2.5345 2.7376 2.9541
a=23.5254 2.3752 2.3752 2.3752 3.2671 2.5794 2.7877 3.0098
a=24.4539 2.4176 2.4176 3.3223 3.333 2.6241 2.8377 3.0653
a=25.4065 2.4608 2.4608 3.3893 3.4 2.6686 2.8875 3.1207
a=26.3836 2.5047 3.4574 3.4574 3.468 2.7128 2.937 3.176
a=27.3854 3.5266 3.5266 3.5266 3.5372 3.7346 2.9863 3.2311
a=28.4123 3.5969 3.5969 3.5969 3.6075 3.8049 3.0353 3.286
a=29.4645 3.6683 3.6683 3.6683 3.679 3.8764 3.0841 3.3407
a=30.5423 3.7409 3.7409 3.7409 3.7516 3.9489 3.1325 3.3951
a=31.6462 3.8146 3.8146 3.8146 3.8253 4.0227 3.1807 3.4493
a=32.7763 3.8895 3.8895 3.8895 3.9001 4.0975 4.3334 3.5032
a=33.933 3.9655 3.9655 3.9655 3.9761 4.1735 4.4094 3.5568
a=35.1166 4.0426 4.0426 4.0426 4.0533 4.2507 4.4866 3.6101
a=36.3274 4.121 4.121 4.121 4.1316 4.329 4.5649 3.663
a=37.5657 4.2005 4.2005 4.2005 4.2111 4.4085 4.6444 3.7156
a=38.8319 4.2811 4.2811 4.2811 4.2918 4.4892 4.7251 3.7688
a=40.1262 4.363 4.363 4.363 4.3736 4.571 4.8069 5.0889
a=41.4489 4.446 4.446 4.446 4.4567 4.654 4.89 5.172
a=42.8004 4.5302 4.5302 4.5302 4.5409 4.7383 4.9742 5.2562
a=44.181 4.6157 4.6157 4.6157 4.6263 4.8237 5.0596 5.3416
a=45.5909 4.7023 4.7023 4.7023 4.7129 4.9103 5.1463 5.4283
a=47.0305 4.7901 4.7901 4.7901 4.8008 4.9982 5.2341 5.5161
a=48.5001 4.8792 4.8792 4.8792 4.8899 5.0872 5.3232 5.6052
a=50 4.9695 4.9695 4.9695 4.9801 5.1775 5.4135 5.6955
Table for Savings/COH at l=0
z1 z2 z3 z4 z5 z6 z7
_______ _______ _______ _______ _______ _______ _______
a=0 0.26015 0.21764 0.18208 0.15233 0.12744 0.35984 0.30105
a=5.15305e-05 0.2015 0.17432 0.15027 0.12909 0.37311 0.31859 0.27133
a=0.000412244 0.34107 0.30929 0.27943 0.25147 0.22542 0.20124 0.1789
a=0.00139132 0.42687 0.39568 0.36572 0.33704 0.30964 0.28356 0.25881
a=0.00329795 0.48958 0.45937 0.43001 0.40155 0.37404 0.34749 0.32195
a=0.00644131 0.53956 0.51034 0.48172 0.45375 0.42647 0.39992 0.37414
a=0.0111306 0.36592 0.55284 0.52501 0.49763 0.47074 0.4444 0.41866
a=0.017675 0.41292 0.39466 0.37655 0.53535 0.50896 0.48297 0.45741
a=0.0263836 0.45432 0.43583 0.41742 0.39913 0.54245 0.51688 0.49162
a=0.0375657 0.49101 0.47249 0.45396 0.43548 0.41707 0.54702 0.52213
a=0.0515305 0.52373 0.50529 0.48679 0.46826 0.44974 0.43127 0.54955
a=0.0685871 0.55306 0.5348 0.51641 0.49794 0.47943 0.4609 0.57435
a=0.0890447 0.5795 0.56147 0.54327 0.52494 0.50651 0.48801 0.46948
a=0.113213 0.60344 0.58569 0.56773 0.54959 0.5313 0.5129 0.49441
a=0.1414 0.62521 0.60777 0.59008 0.57217 0.55407 0.53581 0.51744
a=0.173915 0.52449 0.62798 0.61058 0.59293 0.57505 0.55698 0.53875
a=0.211069 0.54656 0.53273 0.62945 0.61208 0.59445 0.57659 0.55853
a=0.253169 0.56699 0.55329 0.64687 0.62979 0.61242 0.59479 0.57694
a=0.300526 0.58595 0.57239 0.55853 0.64621 0.62912 0.61174 0.5941
a=0.353448 0.60359 0.59018 0.57646 0.66149 0.64467 0.62755 0.61015
a=0.412244 0.62003 0.60679 0.59323 0.67573 0.6592 0.64234 0.62518
a=0.477224 0.63538 0.62233 0.60894 0.59522 0.67279 0.6562 0.63928
a=0.548697 0.64976 0.63689 0.62368 0.7015 0.68553 0.66921 0.65254
a=0.626972 0.66324 0.65057 0.63753 0.62415 0.6975 0.68144 0.66503
a=0.712358 0.6759 0.66343 0.65058 0.72417 0.70876 0.69297 0.67681
a=0.805164 0.68783 0.67554 0.66289 0.64986 0.71938 0.70385 0.68794
a=0.9057 0.69907 0.68698 0.67451 0.66166 0.72941 0.71414 0.69848
a=1.01427 0.70968 0.69779 0.68551 0.67284 0.7389 0.72388 0.70846
a=1.1312 0.71972 0.70802 0.69592 0.68344 0.74788 0.73311 0.71794
a=1.25678 0.72923 0.71771 0.7058 0.6935 0.75639 0.74187 0.72694
a=1.39132 0.73825 0.72692 0.71519 0.70307 0.76448 0.7502 0.7355
a=1.53515 0.74681 0.73566 0.72412 0.71217 0.77218 0.75813 0.74366
a=1.68855 0.75495 0.74398 0.73262 0.72085 0.7795 0.76568 0.75144
a=1.85185 0.7627 0.75191 0.74072 0.72912 0.78648 0.77289 0.75886
a=2.02536 0.77009 0.75947 0.74845 0.73702 0.79314 0.77976 0.76596
a=2.20937 0.77713 0.76668 0.75583 0.74457 0.7995 0.78634 0.77274
a=2.40421 0.78386 0.77358 0.76289 0.75179 0.80558 0.79263 0.77924
a=2.61017 0.7903 0.78017 0.76964 0.75871 0.8114 0.79865 0.78546
a=2.82758 0.79646 0.78648 0.77611 0.76534 0.81697 0.80442 0.79143
a=3.05674 0.80235 0.79253 0.78232 0.7717 0.82232 0.80996 0.79716
a=3.29795 0.80801 0.79834 0.78827 0.7778 0.82745 0.81528 0.80266
a=3.55153 0.81343 0.80391 0.79399 0.78367 0.83237 0.82038 0.80795
a=3.81779 0.81864 0.80926 0.79948 0.7893 0.83711 0.82529 0.81304
a=4.09704 0.82364 0.8144 0.80476 0.79473 0.84166 0.83002 0.81794
a=4.38957 0.82846 0.81935 0.80985 0.79995 0.84604 0.83457 0.82267
a=4.69572 0.83309 0.82411 0.81475 0.80499 0.85026 0.83896 0.82722
a=5.01577 0.83755 0.8287 0.81947 0.80984 0.85433 0.84319 0.83161
a=5.35005 0.84185 0.83313 0.82402 0.81452 0.80463 0.84727 0.83585
a=5.69886 0.846 0.8374 0.82841 0.81904 0.80927 0.8512 0.83994
a=6.06251 0.85 0.84152 0.83266 0.82341 0.81377 0.85501 0.84389
a=6.44131 0.85386 0.8455 0.83675 0.82763 0.81811 0.85869 0.84772
a=6.83557 0.8576 0.84934 0.84072 0.83171 0.82231 0.86224 0.85142
a=7.2456 0.8612 0.85306 0.84455 0.83566 0.82637 0.86569 0.855
a=7.67171 0.86469 0.85666 0.84826 0.83948 0.83031 0.86902 0.85847
a=8.1142 0.86807 0.86014 0.85185 0.84318 0.83412 0.87224 0.86183
a=8.57339 0.87134 0.86352 0.85533 0.84677 0.83782 0.87537 0.86509
a=9.04958 0.87451 0.86679 0.8587 0.85024 0.84141 0.8784 0.86825
a=9.54309 0.87758 0.86995 0.86197 0.85362 0.84489 0.83577 0.87132
a=10.0542 0.88056 0.87303 0.86514 0.85689 0.84826 0.83925 0.87429
a=10.5833 0.88345 0.87601 0.86822 0.86007 0.85154 0.84263 0.87718
a=11.1306 0.88625 0.8789 0.87121 0.86315 0.85472 0.84592 0.87999
a=11.6964 0.84526 0.88171 0.87411 0.86615 0.85782 0.84911 0.88272
a=12.2812 0.84848 0.88445 0.87693 0.86906 0.86083 0.85222 0.88538
a=12.885 0.85161 0.8871 0.87967 0.8719 0.86375 0.85524 0.88796
a=13.5084 0.85462 0.88968 0.88234 0.87465 0.8666 0.85818 0.89047
a=14.1516 0.85749 0.89219 0.88494 0.87733 0.86937 0.86104 0.85233
a=14.8148 0.86025 0.85397 0.88746 0.87994 0.87207 0.86383 0.85521
a=15.4985 0.86289 0.85685 0.88992 0.88249 0.8747 0.86654 0.85802
a=16.2028 0.86542 0.85961 0.89231 0.88496 0.87726 0.86919 0.86075
a=16.9282 0.86784 0.86226 0.89465 0.88737 0.87975 0.87177 0.86342
a=17.675 0.87017 0.8648 0.89692 0.88973 0.88219 0.87429 0.86602
a=18.4433 0.87241 0.86724 0.86115 0.89202 0.88456 0.87674 0.86856
a=19.2337 0.87456 0.86959 0.86372 0.89426 0.88687 0.87913 0.87103
a=20.0462 0.87663 0.87184 0.86618 0.89644 0.88913 0.88147 0.87345
a=20.8814 0.87862 0.874 0.86855 0.89857 0.89133 0.88375 0.87581
a=21.7394 0.88053 0.87609 0.87083 0.90065 0.89349 0.88598 0.87812
a=22.6207 0.88238 0.87809 0.87302 0.90266 0.89559 0.88816 0.88037
a=23.5254 0.88416 0.88002 0.87513 0.86935 0.89764 0.89028 0.88257
a=24.4539 0.88588 0.88188 0.87716 0.87158 0.89965 0.89236 0.88472
a=25.4065 0.88753 0.88368 0.87912 0.87373 0.90161 0.89439 0.88683
a=26.3836 0.88913 0.88541 0.88101 0.8758 0.90353 0.89638 0.88888
a=27.3854 0.89068 0.88708 0.88283 0.87779 0.87184 0.89832 0.8909
a=28.4123 0.89217 0.8887 0.88458 0.87971 0.87396 0.90022 0.89286
a=29.4645 0.89362 0.89026 0.88628 0.88157 0.876 0.90208 0.89479
a=30.5423 0.89502 0.89177 0.88792 0.88336 0.87797 0.9039 0.89668
a=31.6462 0.89637 0.89323 0.8895 0.88509 0.87987 0.90568 0.89853
a=32.7763 0.89768 0.89464 0.89104 0.88676 0.88171 0.87575 0.90033
a=33.933 0.89895 0.89601 0.89252 0.88838 0.88349 0.8777 0.90211
a=35.1166 0.90018 0.89733 0.89395 0.88994 0.8852 0.8796 0.90384
a=36.3274 0.90138 0.89862 0.89534 0.89146 0.88686 0.88142 0.90554
a=37.5657 0.90254 0.89986 0.89669 0.89292 0.88846 0.88319 0.90721
a=38.8319 0.90366 0.90107 0.89799 0.89434 0.89002 0.8849 0.90882
a=40.1262 0.90475 0.90224 0.89926 0.89572 0.89152 0.88656 0.88069
a=41.4489 0.90581 0.90338 0.90048 0.89705 0.89298 0.88816 0.88247
a=42.8004 0.90684 0.90448 0.90167 0.89834 0.89439 0.88971 0.88418
a=44.181 0.90785 0.90555 0.90283 0.89959 0.89575 0.89121 0.88584
a=45.5909 0.90882 0.90659 0.90395 0.90081 0.89708 0.89267 0.88745
a=47.0305 0.90977 0.90761 0.90504 0.90199 0.89837 0.89408 0.88901
a=48.5001 0.91069 0.90859 0.9061 0.90313 0.89962 0.89545 0.89052
a=50 0.91159 0.90955 0.90713 0.90425 0.90083 0.89678 0.89198
Table for Savings/COH at l=1
z1 z2 z3 z4 z5 z6 z7
_______ _______ _______ _______ _______ _______ _______
a=0 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236 0.71236
a=5.15305e-05 0.71232 0.71232 0.71232 0.71232 0.71232 0.71232 0.71232
a=0.000412244 0.71206 0.71206 0.71206 0.71206 0.71206 0.71206 0.71206
a=0.00139132 0.71135 0.71135 0.71135 0.71135 0.71135 0.71135 0.71135
a=0.00329795 0.70997 0.70997 0.70997 0.70997 0.70997 0.70997 0.70997
a=0.00644131 0.70771 0.70771 0.70771 0.70771 0.70771 0.70771 0.70771
a=0.0111306 0.70436 0.70436 0.70436 0.70436 0.70436 0.70436 0.70436
a=0.017675 0.69974 0.69974 0.69974 0.69974 0.69974 0.69974 0.69974
a=0.0263836 0.69369 0.69369 0.69369 0.69369 0.69369 0.69369 0.69369
a=0.0375657 0.68607 0.68607 0.68607 0.68607 0.68607 0.68607 0.68607
a=0.0515305 0.67679 0.67679 0.67679 0.67679 0.67679 0.67679 0.67679
a=0.0685871 0.66578 0.66578 0.66578 0.66578 0.66578 0.66578 0.66578
a=0.0890447 0.73812 0.73812 0.65304 0.65304 0.65304 0.65304 0.65304
a=0.113213 0.72181 0.72181 0.72181 0.72181 0.72181 0.72181 0.72181
a=0.1414 0.70367 0.70367 0.70367 0.70367 0.70367 0.70367 0.70367
a=0.173915 0.68385 0.68385 0.68385 0.68385 0.68385 0.68385 0.68385
a=0.211069 0.74525 0.74525 0.66253 0.66253 0.66253 0.66253 0.66253
a=0.253169 0.71982 0.71982 0.71982 0.71982 0.71982 0.71982 0.71982
a=0.300526 0.69321 0.69321 0.69321 0.69321 0.69321 0.69321 0.69321
a=0.353448 0.74551 0.74551 0.74551 0.74551 0.74551 0.74551 0.6657
a=0.412244 0.71403 0.71403 0.71403 0.71403 0.71403 0.71403 0.71403
a=0.477224 0.76084 0.76084 0.76084 0.6822 0.6822 0.6822 0.6822
a=0.548697 0.72528 0.72528 0.72528 0.72528 0.72528 0.72528 0.72528
a=0.626972 0.76656 0.76656 0.76656 0.76656 0.76656 0.76656 0.68996
a=0.712358 0.72789 0.72789 0.72789 0.72789 0.72789 0.72789 0.72789
a=0.805164 0.76393 0.76393 0.76393 0.76393 0.76393 0.76393 0.76393
a=0.9057 0.72321 0.72321 0.72321 0.72321 0.72321 0.72321 0.72321
a=1.01427 0.75453 0.75453 0.75453 0.75453 0.75453 0.75453 0.75453
a=1.1312 0.78398 0.78398 0.78398 0.78398 0.78398 0.78398 0.71275
a=1.25678 0.73997 0.73997 0.73997 0.73997 0.73997 0.73997 0.73997
a=1.39132 0.7655 0.7655 0.7655 0.7655 0.7655 0.7655 0.7655
a=1.53515 0.78935 0.78935 0.78935 0.78935 0.78935 0.78935 0.78935
a=1.68855 0.81158 0.81158 0.74398 0.74398 0.74398 0.74398 0.74398
a=1.85185 0.76478 0.76478 0.76478 0.76478 0.76478 0.76478 0.76478
a=2.02536 0.78419 0.78419 0.78419 0.78419 0.78419 0.78419 0.78419
a=2.20937 0.80225 0.80225 0.80225 0.80225 0.80225 0.80225 0.80225
a=2.40421 0.81904 0.81904 0.81904 0.81904 0.81904 0.81904 0.81904
a=2.61017 0.77206 0.77206 0.77206 0.77206 0.77206 0.77206 0.77206
a=2.82758 0.78698 0.78698 0.78698 0.78698 0.78698 0.78698 0.78698
a=3.05674 0.80089 0.80089 0.80089 0.80089 0.80089 0.80089 0.80089
a=3.29795 0.81384 0.81384 0.81384 0.81384 0.81384 0.81384 0.81384
a=3.55153 0.8259 0.8259 0.8259 0.8259 0.8259 0.8259 0.8259
a=3.81779 0.83713 0.83713 0.83713 0.83713 0.83713 0.83713 0.83713
a=4.09704 0.84758 0.84758 0.84758 0.84758 0.84758 0.84758 0.84758
a=4.38957 0.8014 0.8014 0.8014 0.8014 0.8014 0.8014 0.8014
a=4.69572 0.81106 0.81106 0.81106 0.81106 0.81106 0.81106 0.81106
a=5.01577 0.82009 0.82009 0.82009 0.82009 0.82009 0.82009 0.82009
a=5.35005 0.82856 0.82856 0.82856 0.82856 0.82856 0.82856 0.82856
a=5.69886 0.83649 0.83649 0.83649 0.83649 0.83649 0.83649 0.83649
a=6.06251 0.84392 0.84392 0.84392 0.84392 0.84392 0.84392 0.84389
a=6.44131 0.85088 0.85088 0.85088 0.85088 0.85088 0.85088 0.84772
a=6.83557 0.85742 0.85742 0.85742 0.85742 0.85742 0.85742 0.85142
a=7.2456 0.86355 0.86355 0.86355 0.86355 0.86355 0.86355 0.855
a=7.67171 0.8693 0.8693 0.8693 0.8693 0.8693 0.86902 0.85847
a=8.1142 0.87471 0.87471 0.87471 0.87471 0.82701 0.87224 0.86183
a=8.57339 0.83266 0.83266 0.83266 0.83266 0.83266 0.87537 0.86509
a=9.04958 0.83802 0.83802 0.83802 0.83802 0.83802 0.8784 0.86825
a=9.54309 0.84308 0.84308 0.84308 0.84308 0.84308 0.83577 0.87132
a=10.0542 0.84787 0.84787 0.84787 0.84787 0.84787 0.83925 0.87429
a=10.5833 0.85242 0.85242 0.85242 0.85242 0.85154 0.84263 0.87718
a=11.1306 0.85672 0.85672 0.85672 0.85672 0.85472 0.84592 0.87999
a=11.6964 0.86081 0.86081 0.86081 0.86081 0.85782 0.84911 0.88272
a=12.2812 0.86469 0.86469 0.86469 0.86469 0.86083 0.85222 0.88538
a=12.885 0.86837 0.86837 0.86837 0.86837 0.86375 0.85524 0.88796
a=13.5084 0.87187 0.87187 0.87187 0.87187 0.8666 0.85818 0.89047
a=14.1516 0.8752 0.8752 0.8752 0.8752 0.86937 0.86104 0.85233
a=14.8148 0.87837 0.87837 0.87837 0.87837 0.87207 0.86383 0.85521
a=15.4985 0.88139 0.88139 0.88139 0.88139 0.8747 0.86654 0.85802
a=16.2028 0.88427 0.88427 0.88427 0.88427 0.87726 0.86919 0.86075
a=16.9282 0.88701 0.88701 0.88701 0.88701 0.87975 0.87177 0.86342
a=17.675 0.88963 0.88963 0.88963 0.88963 0.88219 0.87429 0.86602
a=18.4433 0.89213 0.89213 0.89213 0.89202 0.88456 0.87674 0.86856
a=19.2337 0.89451 0.89451 0.89451 0.89426 0.88687 0.87913 0.87103
a=20.0462 0.89679 0.89679 0.89679 0.89644 0.88913 0.88147 0.87345
a=20.8814 0.89897 0.89897 0.89897 0.89857 0.89133 0.88375 0.87581
a=21.7394 0.90106 0.90106 0.90106 0.90065 0.89349 0.88598 0.87812
a=22.6207 0.90306 0.90306 0.90306 0.90266 0.89559 0.88816 0.88037
a=23.5254 0.90498 0.90498 0.90498 0.86935 0.89764 0.89028 0.88257
a=24.4539 0.90681 0.90681 0.87194 0.87158 0.89965 0.89236 0.88472
a=25.4065 0.90857 0.90857 0.87407 0.87373 0.90161 0.89439 0.88683
a=26.3836 0.91026 0.87613 0.87613 0.8758 0.90353 0.89638 0.88888
a=27.3854 0.87811 0.87811 0.87811 0.87779 0.87184 0.89832 0.8909
a=28.4123 0.88003 0.88003 0.88003 0.87971 0.87396 0.90022 0.89286
a=29.4645 0.88187 0.88187 0.88187 0.88157 0.876 0.90208 0.89479
a=30.5423 0.88365 0.88365 0.88365 0.88336 0.87797 0.9039 0.89668
a=31.6462 0.88537 0.88537 0.88537 0.88509 0.87987 0.90568 0.89853
a=32.7763 0.88704 0.88704 0.88704 0.88676 0.88171 0.87575 0.90033
a=33.933 0.88865 0.88865 0.88865 0.88838 0.88349 0.8777 0.90211
a=35.1166 0.8902 0.8902 0.8902 0.88994 0.8852 0.8796 0.90384
a=36.3274 0.89171 0.89171 0.89171 0.89146 0.88686 0.88142 0.90554
a=37.5657 0.89316 0.89316 0.89316 0.89292 0.88846 0.88319 0.90721
a=38.8319 0.89458 0.89458 0.89458 0.89434 0.89002 0.8849 0.90882
a=40.1262 0.89594 0.89594 0.89594 0.89572 0.89152 0.88656 0.88069
a=41.4489 0.89727 0.89727 0.89727 0.89705 0.89298 0.88816 0.88247
a=42.8004 0.89855 0.89855 0.89855 0.89834 0.89439 0.88971 0.88418
a=44.181 0.8998 0.8998 0.8998 0.89959 0.89575 0.89121 0.88584
a=45.5909 0.90101 0.90101 0.90101 0.90081 0.89708 0.89267 0.88745
a=47.0305 0.90218 0.90218 0.90218 0.90199 0.89837 0.89408 0.88901
a=48.5001 0.90332 0.90332 0.90332 0.90313 0.89962 0.89545 0.89052
a=50 0.90443 0.90443 0.90443 0.90425 0.90083 0.89678 0.89198
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'); 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'); 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_2OC'; profsave(profile('info'), st_file_name); end
ans(:,:,1) =
-112.6174 -103.3361 -95.1882 -87.9347 -81.4030 -75.2255 -68.5180
-103.6891 -96.6424 -90.1347 -84.1009 -78.1601 -71.6884 -66.0478
-83.2627 -78.4528 -73.9510 -69.7187 -65.7201 -61.9379 -58.3568
-68.1110 -64.5537 -61.1999 -58.0291 -55.0170 -52.1542 -49.4311
-57.5919 -54.8005 -52.1318 -49.5895 -47.1621 -44.8474 -42.6405
-49.8681 -47.6630 -45.4765 -43.3598 -41.3212 -39.3674 -37.4990
-43.7478 -42.1362 -40.3593 -38.5733 -36.8244 -35.1325 -33.5057
-38.8730 -37.5559 -36.2348 -34.7462 -33.2398 -31.7557 -30.3150
-35.0351 -33.8655 -32.7147 -31.5647 -30.2894 -28.9854 -27.6981
-31.9535 -30.8887 -29.8444 -28.8137 -27.7820 -26.6544 -25.5012
-29.4230 -28.4395 -27.4755 -26.5250 -25.5805 -24.6447 -23.6184
-27.3037 -26.3859 -25.4861 -24.5989 -23.7168 -22.8460 -21.9776
-25.4947 -24.6340 -23.7882 -22.9531 -22.1221 -21.3020 -20.4971
-23.8980 -23.0844 -22.2873 -21.5012 -20.7190 -19.9465 -19.1882
-22.4981 -21.7229 -20.9634 -20.2143 -19.4684 -18.7313 -18.0079
-21.2145 -20.5158 -19.7906 -19.0735 -18.3588 -17.6523 -16.9588
-20.0378 -19.4057 -18.7404 -18.0540 -17.3669 -16.6871 -16.0196
-18.9634 -18.3779 -17.7457 -17.0876 -16.4282 -15.7751 -15.1336
-17.9649 -17.4116 -16.8368 -16.2102 -15.5732 -14.9410 -14.3198
-17.0440 -16.5184 -15.9796 -15.3987 -14.7881 -14.1799 -13.5804
-16.1875 -15.6863 -15.1761 -14.6167 -14.0232 -13.4317 -12.8499
-15.3834 -14.9016 -14.4120 -13.8973 -13.3313 -12.7570 -12.1904
-14.6343 -14.1703 -13.7007 -13.2051 -12.6519 -12.0960 -11.5472
-13.9187 -13.4693 -13.0155 -12.5454 -12.0232 -11.4838 -10.9487
-13.2450 -12.8092 -12.3705 -11.9179 -11.3913 -10.8603 -10.3360
-12.5966 -12.1722 -11.7455 -11.3066 -10.7995 -10.2821 -9.7711
-11.9755 -11.5622 -11.1469 -10.7203 -10.2249 -9.7169 -9.2149
-11.3809 -10.9767 -10.5710 -10.1542 -9.6685 -9.1732 -8.6840
-10.8016 -10.4068 -10.0105 -9.6034 -9.1460 -8.6645 -8.1858
-10.2469 -9.8610 -9.4738 -9.0775 -8.6195 -8.1438 -7.6722
-9.7088 -9.3303 -8.9507 -8.5618 -8.1041 -7.6356 -7.1726
-9.1810 -8.8099 -8.4377 -8.0557 -7.6144 -7.1569 -6.7041
-8.6685 -8.3045 -7.9397 -7.5656 -7.1399 -6.6898 -6.2431
-8.1750 -7.8179 -7.4595 -7.0921 -6.6635 -6.2198 -5.7804
-7.6879 -7.3365 -6.9842 -6.6231 -6.2043 -5.7693 -5.3382
-7.2106 -6.8650 -6.5190 -6.1645 -5.7583 -5.3304 -4.9056
-6.7463 -6.4065 -6.0664 -5.7185 -5.3249 -4.9049 -4.4877
-6.2946 -5.9603 -5.6260 -5.2845 -4.8963 -4.4805 -4.0672
-5.8553 -5.5254 -5.1957 -4.8590 -4.4714 -4.0620 -3.6558
-5.4189 -5.0937 -4.7689 -4.4372 -4.0600 -3.6592 -3.2606
-4.9893 -4.6688 -4.3491 -4.0233 -3.6623 -3.2707 -2.8794
-4.5689 -4.2533 -3.9388 -3.6198 -3.2769 -2.8940 -2.5093
-4.1595 -3.8487 -3.5396 -3.2278 -2.9013 -2.5254 -2.1460
-3.7619 -3.4557 -3.1519 -2.8470 -2.5312 -2.1598 -1.7840
-3.3756 -3.0736 -2.7745 -2.4752 -2.1596 -1.7901 -1.4167
-2.9968 -2.6981 -2.4026 -2.1069 -1.7927 -1.4262 -1.0562
-2.6220 -2.3266 -2.0345 -1.7421 -1.4342 -1.0723 -0.7067
-2.2524 -1.9603 -1.6715 -1.3827 -1.0838 -0.7282 -0.3676
-1.8888 -1.5998 -1.3145 -1.0297 -0.7376 -0.3913 -0.0367
-1.5320 -1.2461 -0.9641 -0.6833 -0.3970 -0.0623 0.2856
-1.1825 -0.8996 -0.6209 -0.3440 -0.0629 0.2612 0.6152
-0.8407 -0.5607 -0.2851 -0.0120 0.2643 0.5763 0.9353
-0.5069 -0.2296 0.0429 0.3125 0.5842 0.8830 1.2458
-0.1812 0.0934 0.3630 0.6293 0.8965 1.1830 1.5487
0.1363 0.4084 0.6752 0.9383 1.2013 1.4870 1.8538
0.4455 0.7153 0.9795 1.2397 1.5001 1.7851 2.1517
0.7480 1.0160 1.2783 1.5365 1.7951 2.0771 2.4424
1.0457 1.3120 1.5726 1.8289 2.0857 2.3642 2.7268
1.3385 1.6032 1.8620 2.1165 2.3717 2.6490 3.0067
1.6261 1.8893 2.1464 2.3992 2.6530 2.9303 3.2815
1.9084 2.1701 2.4257 2.6771 2.9305 3.2078 3.5512
2.1876 2.4456 2.6998 2.9501 3.2043 3.4812 3.8156
2.4629 2.7158 2.9687 3.2183 3.4738 3.7502 4.0746
2.7340 2.9806 3.2324 3.4816 3.7387 4.0145 4.3279
3.0010 3.2402 3.4909 3.7398 3.9987 4.2739 4.5757
3.2635 3.4944 3.7441 3.9930 4.2536 4.5282 4.8199
3.5217 3.7439 3.9921 4.2411 4.5033 4.7774 5.0623
3.7755 3.9919 4.2349 4.4840 4.7478 5.0214 5.3017
4.0250 4.2351 4.4728 4.7218 4.9870 5.2603 5.5374
4.2705 4.4753 4.7057 4.9545 5.2210 5.4941 5.7688
4.5119 4.7113 4.9338 5.1821 5.4497 5.7226 5.9957
4.7496 4.9437 5.1571 5.4048 5.6731 5.9460 6.2178
4.9833 5.1721 5.3799 5.6226 5.8914 6.1642 6.4350
5.2133 5.3966 5.5982 5.8361 6.1047 6.3773 6.6473
5.4396 5.6177 5.8142 6.0451 6.3131 6.5854 6.8547
5.6623 5.8352 6.0264 6.2498 6.5167 6.7884 7.0571
5.8814 6.0494 6.2356 6.4503 6.7156 6.9866 7.2547
6.0970 6.2602 6.4412 6.6500 6.9098 7.1800 7.4474
6.3092 6.4677 6.6436 6.8461 7.1000 7.3688 7.6354
6.5180 6.6719 6.8431 7.0407 7.2862 7.5532 7.8188
6.7236 6.8731 7.0396 7.2318 7.4687 7.7332 7.9978
6.9260 7.0713 7.2337 7.4209 7.6481 7.9092 8.1724
7.1253 7.2669 7.4250 7.6068 7.8268 8.0813 8.3428
7.3219 7.4598 7.6138 7.7905 8.0026 8.2501 8.5092
7.5158 7.6501 7.7998 7.9714 8.1771 8.4157 8.6720
7.7070 7.8378 7.9834 8.1501 8.3489 8.5785 8.8313
7.8957 8.0228 8.1644 8.3261 8.5190 8.7392 8.9873
8.0816 8.2053 8.3430 8.5000 8.6866 8.8991 9.1403
8.2650 8.3853 8.5190 8.6714 8.8523 9.0569 9.2909
8.4458 8.5627 8.6927 8.8405 9.0157 9.2136 9.4390
8.6241 8.7377 8.8639 9.0074 9.1772 9.3682 9.5851
8.7998 8.9102 9.0328 9.1721 9.3365 9.5215 9.7291
8.9731 9.0803 9.1993 9.3345 9.4940 9.6728 9.8727
9.1439 9.2481 9.3636 9.4948 9.6494 9.8226 10.0148
9.3123 9.4134 9.5256 9.6529 9.8028 9.9705 10.1561
9.4782 9.5765 9.6854 9.8089 9.9543 10.1168 10.2959
9.6418 9.7372 9.8430 9.9629 10.1040 10.2614 10.4346
9.8031 9.8958 9.9984 10.1148 10.2517 10.4044 10.5718
9.9621 10.0520 10.1517 10.2647 10.3975 10.5456 10.7078
10.1188 10.2061 10.3029 10.4126 10.5415 10.6852 10.8423
ans(:,:,2) =
-12.8083 -12.4821 -12.1518 -11.8074 -11.4237 -11.0294 -10.6360
-12.8080 -12.4817 -12.1514 -11.8071 -11.4234 -11.0290 -10.6356
-12.8056 -12.4794 -12.1491 -11.8047 -11.4210 -11.0266 -10.6333
-12.7991 -12.4729 -12.1426 -11.7983 -11.4146 -11.0202 -10.6268
-12.7867 -12.4605 -12.1302 -11.7858 -11.4021 -11.0077 -10.6144
-12.7664 -12.4402 -12.1099 -11.7656 -11.3819 -10.9875 -10.5941
-12.7368 -12.4106 -12.0803 -11.7360 -11.3523 -10.9579 -10.5645
-12.6967 -12.3704 -12.0401 -11.6958 -11.3121 -10.9177 -10.5243
-12.6452 -12.3190 -11.9887 -11.6444 -11.2607 -10.8663 -10.4729
-12.5823 -12.2560 -11.9257 -11.5814 -11.1977 -10.8033 -10.4100
-12.5082 -12.1819 -11.8516 -11.5073 -11.1236 -10.7292 -10.3358
-12.4237 -12.0975 -11.7672 -11.4228 -11.0391 -10.6448 -10.2514
-12.3204 -12.0022 -11.6737 -11.3293 -10.9456 -10.5513 -10.1579
-12.1687 -11.8505 -11.5287 -11.1937 -10.8211 -10.4362 -10.0520
-12.0131 -11.6949 -11.3731 -11.0381 -10.6655 -10.2806 -9.8964
-11.8565 -11.5383 -11.2165 -10.8815 -10.5089 -10.1240 -9.7398
-11.6896 -11.3779 -11.0614 -10.7264 -10.3538 -9.9689 -9.5847
-11.4636 -11.1519 -10.8368 -10.5088 -10.1429 -9.7637 -9.3848
-11.2541 -10.9424 -10.6273 -10.2993 -9.9335 -9.5542 -9.1754
-11.0183 -10.7133 -10.4051 -10.0844 -9.7267 -9.3567 -8.9820
-10.7574 -10.4524 -10.1442 -9.8235 -9.4658 -9.0958 -8.7262
-10.5149 -10.2166 -9.9152 -9.5951 -9.2374 -8.8674 -8.4978
-10.2163 -9.9180 -9.6166 -9.3043 -8.9627 -8.6033 -8.2412
-9.9455 -9.6524 -9.3565 -9.0495 -8.7087 -8.3521 -7.9910
-9.6306 -9.3374 -9.0415 -8.7345 -8.3937 -8.0371 -7.6791
-9.3192 -9.0314 -8.7409 -8.4387 -8.0995 -7.7474 -7.3953
-9.0098 -8.7220 -8.4315 -8.1293 -7.7901 -7.4380 -7.0859
-8.6653 -8.3829 -8.0979 -7.8017 -7.4727 -7.1292 -6.7848
-8.3478 -8.0705 -7.7901 -7.4991 -7.1782 -6.8396 -6.4975
-8.0121 -7.7348 -7.4544 -7.1634 -6.8425 -6.5039 -6.1634
-7.6535 -7.3803 -7.1045 -6.8180 -6.4983 -6.1639 -5.8285
-7.3067 -7.0380 -6.7669 -6.4856 -6.1730 -5.8454 -5.5163
-6.9764 -6.7119 -6.4421 -6.1608 -5.8482 -5.5206 -5.1915
-6.6101 -6.3456 -6.0791 -5.8030 -5.4976 -5.1753 -4.8507
-6.2493 -5.9891 -5.7270 -5.4560 -5.1581 -4.8419 -4.5231
-5.8999 -5.6429 -5.3842 -5.1168 -4.8221 -4.5085 -4.1920
-5.5529 -5.2995 -5.0445 -4.7809 -4.4889 -4.1803 -3.8693
-5.2115 -4.9580 -4.7031 -4.4394 -4.1475 -3.8388 -3.5279
-4.8413 -4.5915 -4.3406 -4.0817 -3.7978 -3.4963 -3.1915
-4.4780 -4.2321 -3.9854 -3.7324 -3.4593 -3.1658 -2.8670
-4.1237 -3.8817 -3.6395 -3.3927 -3.1310 -2.8449 -2.5515
-3.7793 -3.5410 -3.3032 -3.0624 -2.8104 -2.5301 -2.2408
-3.4443 -3.2094 -2.9755 -2.7394 -2.4933 -2.2164 -1.9298
-3.1151 -2.8828 -2.6517 -2.4184 -2.1729 -1.8971 -1.6117
-2.7817 -2.5494 -2.3183 -2.0850 -1.8395 -1.5636 -1.2782
-2.4386 -2.2088 -1.9804 -1.7497 -1.5056 -1.2318 -0.9490
-2.0987 -1.8715 -1.6458 -1.4179 -1.1782 -0.9082 -0.6288
-1.7631 -1.5385 -1.3154 -1.0910 -0.8573 -0.5925 -0.3173
-1.4326 -1.2105 -0.9902 -0.7692 -0.5406 -0.2824 -0.0125
-1.1079 -0.8882 -0.6706 -0.4529 -0.2290 0.0215 0.2856
-0.7893 -0.5720 -0.3571 -0.1426 0.0771 0.3213 0.6152
-0.4773 -0.2623 -0.0500 0.1616 0.3772 0.6143 0.9353
-0.1720 0.0407 0.2506 0.4593 0.6709 0.9002 1.2458
0.1263 0.3369 0.5445 0.7504 0.9583 1.1830 1.5487
0.4176 0.6262 0.8316 1.0351 1.2417 1.4870 1.8538
0.7065 0.9151 1.1205 1.3240 1.5297 1.7851 2.1517
0.9975 1.2047 1.4086 1.6104 1.8148 2.0771 2.4424
1.2848 1.4906 1.6929 1.8933 2.0963 2.3642 2.7268
1.5679 1.7724 1.9733 2.1722 2.3740 2.6490 3.0067
1.8468 2.0499 2.2495 2.4470 2.6530 2.9303 3.2815
2.1211 2.3230 2.5213 2.7177 2.9305 3.2078 3.5512
2.3907 2.5915 2.7887 2.9843 3.2043 3.4812 3.8156
2.6567 2.8554 3.0514 3.2466 3.4738 3.7502 4.0746
2.9188 3.1145 3.3096 3.5046 3.7387 4.0145 4.3279
3.1768 3.3690 3.5632 3.7582 3.9987 4.2739 4.5757
3.4308 3.6188 3.8121 4.0072 4.2536 4.5282 4.8199
3.6805 3.8638 4.0563 4.2517 4.5033 4.7774 5.0623
3.9261 4.1044 4.2959 4.4915 4.7478 5.0214 5.3017
4.1681 4.3421 4.5309 4.7266 4.9870 5.2603 5.5374
4.4061 4.5758 4.7615 4.9570 5.2210 5.4941 5.7688
4.6406 4.8061 4.9877 5.1828 5.4497 5.7226 5.9957
4.8714 5.0324 5.2094 5.4048 5.6731 5.9460 6.2178
5.0986 5.2551 5.4269 5.6226 5.8914 6.1642 6.4350
5.3223 5.4745 5.6421 5.8361 6.1047 6.3773 6.6473
5.5425 5.6905 5.8537 6.0451 6.3131 6.5854 6.8547
5.7593 5.9032 6.0624 6.2498 6.5167 6.7884 7.0571
5.9728 6.1127 6.2676 6.4503 6.7156 6.9866 7.2547
6.1831 6.3190 6.4696 6.6500 6.9098 7.1800 7.4474
6.3901 6.5222 6.6701 6.8461 7.1000 7.3688 7.6354
6.5940 6.7224 6.8687 7.0407 7.2862 7.5532 7.8188
6.7948 6.9218 7.0645 7.2318 7.4687 7.7332 7.9978
6.9935 7.1185 7.2579 7.4209 7.6481 7.9092 8.1724
7.1908 7.3127 7.4484 7.6068 7.8268 8.0813 8.3428
7.3853 7.5041 7.6365 7.7905 8.0026 8.2501 8.5092
7.5773 7.6931 7.8218 7.9714 8.1771 8.4157 8.6720
7.7667 7.8795 8.0048 8.1501 8.3489 8.5785 8.8313
7.9535 8.0633 8.1851 8.3261 8.5190 8.7392 8.9873
8.1377 8.2446 8.3631 8.5000 8.6866 8.8991 9.1403
8.3194 8.4233 8.5385 8.6714 8.8523 9.0569 9.2909
8.4986 8.5997 8.7117 8.8405 9.0157 9.2136 9.4390
8.6753 8.7736 8.8823 9.0074 9.1772 9.3682 9.5851
8.8495 8.9451 9.0507 9.1721 9.3365 9.5215 9.7291
9.0213 9.1141 9.2167 9.3345 9.4940 9.6728 9.8727
9.1907 9.2809 9.3806 9.4948 9.6494 9.8226 10.0148
9.3577 9.4453 9.5421 9.6529 9.8028 9.9705 10.1561
9.5224 9.6075 9.7014 9.8089 9.9543 10.1168 10.2959
9.6847 9.7673 9.8585 9.9629 10.1040 10.2614 10.4346
9.8447 9.9250 10.0135 10.1148 10.2517 10.4044 10.5718
10.0025 10.0805 10.1664 10.2647 10.3975 10.5456 10.7078
10.1581 10.2338 10.3172 10.4126 10.5415 10.6852 10.8423