%% Task 5.0
clear all
load deciles


%% Task 5.1
N = 10;
count=0;
for i=1:N
    count=count+i;
end
count

%% Task 5.2 

r=.05;
cr=100*ones(10,1);
cr(1)=cr(1)*(1+r);
for i=2:10
    cr(i)=cr(i-1)*(1+r);
end
cr

plot(y)


%% Task 5.3
clear all
load deciles
deciles=[dec01 dec02 dec03 dec04 dec05 dec06 dec07 dec08 dec09 dec10];
T=size(deciles,1);
K=size(deciles,2);
gr=deciles+1;
cr=zeros(size(deciles));
cr(1,:)=gr(1,:);
for i=2:T
    for j=1:K
        cr(i,j)=cr(i-1,j)*gr(i,j);
    end
end
plot(cr)

%% Task 5.4
T=1000;
K=10;
e=randn(T,K);
rw1=zeros(T,K);
rw1(1,:)=e(1,:);
for i=2:T
    for j=1:K
        rw1(i,j)=rw1(i-1,j)+e(i,j);
    end
end
    
rw2=zeros(1000,10);
rw2(1,:)=e(1,:)
for i=2:T
    rw2(i,:)=rw2(i-1,:)+e(i,:);
end


rw3=cumsum(e);

plot(rw1);
max(max(abs(rw1-rw2)))
max(max(abs(rw1-rw3)))