function [req]=mc8(sn)


tm1=cputime;


% Anzahl der Simulationen
sn;

% betrachteter Zeitraum in Jahren jn
jn=71;


% Zinssatz r
r=0.02;


% Zufallsgenerator Aktien
epseq=randn(jn,sn);
dt=1;
mueeq=0.06;
sigeq=0.18;
for s1=1:jn
req(s1,:)=(mueeq-1/2*sigeq^2)*dt+sigeq*epseq(s1,:)*dt^0.5;
end;

% Selektion
r8(1,:)=req(1,:);
r8(2,:)=req(11,:);
r8(3,:)=req(21,:);
r8(4,:)=req(31,:);
r8(5,:)=req(41,:);
r8(6,:)=req(51,:);
r8(7,:)=req(61,:);
r8(8,:)=req(71,:);
req=r8;


% Zufallsgenerator Arbeitseinkommen
alter=(20:99)';

% Arbeitseinkommen stochastischer Anteil 1a: permanente Schocks xsikorr
xsi=randn(jn,sn);
rhoeql=0.0058;
xsikorr=(rhoeql*epseq+(1-rhoeql^2)^0.5*xsi);
% Arbeitseinkommen stochastischer Anteil 1b: permanente Schocks omega
omega=randn(jn,sn);
% stochastischer Anteil 1: permanente Schocks uit, sigma* 0.0106
uit=(xsikorr+omega)*0.0106;
%uit=(xsikorr+omega)*0.00106*2;
%uit=(xsikorr+omega)*0.00106*4.5;
% stochastischer Anteil 1: permanente Schocks vit
vit=zeros(jn,sn);
vit(1,:)=uit(1,:);
for s4=2:jn
vit(s4,:)=vit(s4-1,:)+uit(s4,:);
end;

% Arbeitseinkommen stochastischer Anteil 2: temporaere Schocks eps, sigma=0.0438
epsl=randn(jn,sn)*0.0438;
%epsl=randn(jn,sn)*0.00438*2;
%epsl=randn(jn,sn)*0.00438*4.5;

% deterministischer Anteil 1: mit großer Familie -0.0236 und Ehestatus +0.4437
eit=ones(jn,sn)*(2.7004-0.0236+0.4437);
% deterministischer Anteil 2: mit mittlerer Ausbildung und abhaengig vom Alter
fit=zeros(jn,sn);
jm=min(jn,46);
for s5=1:jm
fit(s5,:)=-2.17+0.1682*alter(s5,1)-0.0323*(alter(s5,1)^2)/10+0.0020*(alter(s5,1)^3)/100;
%fit(s5,:)=-1.87+0.1522*alter(s5,1)-0.0293*(alter(s5,1)^2)/10+0.0018*(alter(s5,1)^3)/100;
end;
for s5=47:jn
fit(s5,:)=fit(jm,:);
end;

% das Arbeitseinkommen 
jm=min(jn,46);
L=zeros(jm,sn);
for s6=1:jm
L(s6,:)=1000*exp(eit(s6,:)+fit(s6,:)+vit(s6,:)+epsl(s6,:));
end;
for s6=47:jn
L(s6,:)=0.68212*L(jm,:);
end;

% Selektion
L8(1,:)=L(1,:);
L8(2,:)=L(11,:);
L8(3,:)=L(21,:);
L8(4,:)=L(31,:);
L8(5,:)=L(41,:);
L8(6,:)=L(51,:);
L8(7,:)=L(61,:);
L8(8,:)=L(71,:);
L=L8;

Lm=mean(L,2);

jn=8;


% Aktienindex, beginnend in t-1 mit XEX=100
XEX=zeros(jn,sn);
XEX(1,:)=100;
for s7=2:(jn+1)
XEX(s7,:)=exp(req(s7-1,:)).*XEX(s7-1,:);
end;

Xm=mean(XEX,2);

% sichere Anlage, beginnend in t-1 mit RSX=100
RSX=zeros(jn,1);
RSX(1)=100;
for s9=2:(jn+1)
RSX(s9)=exp(r).*RSX(s9-1);
end;


tm2=cputime-tm1


save RLdat sn jn r RSX mueeq sigeq req XEX Xm L Lm


% Funktionsplot
close all

high1=max(max(XEX,[],2));
low1=min(min(XEX,[],2));
oben1=high1+0.175*abs(high1-low1);

figure(1)
semilogy((1:(jn+1)),XEX)
axis([0 (jn+1) 25 oben1]);
grid("on");
title('log Aktienindex', "fontsize", 18);
xlabel('Perioden', "fontsize", 18);
ylabel('log Indexstand', "fontsize", 18);
print -mono -deps -FHelvetica:18 bildmc1.eps

high2=max(Xm);
low2=min(Xm);
oben2=high2+0.175*abs(high2-low2);

figure(2)
semilogy((1:(jn+1)),Xm, "linewidth", 3, (1:(jn+1)),RSX, "linewidth", 3)
axis([0 (jn+1) 25 oben2]);
grid("on");
title('log durchschnittlicher Aktienindex und Fixverzinsung', "fontsize", 18);
xlabel('Perioden', "fontsize", 18);
ylabel('log Indexstand', "fontsize", 18);
print -mono -deps -FHelvetica:18 bildmc2.eps

high3=max(max(L,[],2));
low3=min(min(L,[],2));
oben3=high3+0.175*abs(high3-low3);
unten3=low3-0.175*abs(high3-low3);

figure(3)
plot((1:jn),L)
axis([0 (jn+1) unten3 oben3]);
grid("on");
title('Arbeitseinkommen', "fontsize", 18);
xlabel('Perioden', "fontsize", 18);
ylabel('Arbeitseinkommen', "fontsize", 18);
print -mono -deps -FHelvetica:18 bildmc3.eps

high4=max(Lm);
low4=min(Lm);
oben4=high4+0.175*abs(high4-low4);
unten4=low4-0.175*abs(high4-low4);

figure(4)
plot((1:jn),Lm, "linewidth", 3)
axis([0 (jn+1) unten4 oben4]);
grid("on");
title('durchschnittliches Arbeitseinkommen', "fontsize", 18);
xlabel('Perioden', "fontsize", 18);
ylabel('Arbeitseinkommen', "fontsize", 18);
print -mono -deps -FHelvetica:18 bildmc4.eps

end;
