% data for optimal investment to fund an excess stream problem
n  = 6;
T  = 12;
rp = 0.05/12;
rn = 0.09/12;

E  = [1 1 8 5 1 2 2 8 6 0 8 1]';                   % expense
C  = [0.0060 0.0050 0.0030 0.0070 0.0050 0.0040]'; % coupon
P  = [0.9870 0.9805 0.9761 0.9946 0.9783 0.9680]'; % bond price
M  = [     3      4      6      6      9     10]'; % time to maturity


A  = zeros(T,n); % A: (time x coupon_type) A=[a_{t,i}]
for i = 1:n
    A(1:M(i)-1,i) = C(i);
    A(M(i)    ,i) = C(i)+1;
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% computing bond prices
%  use this formula to generate a bond price vector P
%
% Y = [0.125 0.120 0.085 0.095 0.090 0.088]'/12; % bond yield
% P = C.*(1-(1+Y).^-M)./Y+(1+Y).^-M;             % bond price