% Dieses Programm erzeugt zufaellig gestoerte Daten fuer den 
% radioaktiven Zerfallsprozess, der in der Aufgabe beschrieben ist.
%
% Die "exakte" Isotopenverteilung zum Zeitpunkt t=0 ist 7.0, 2.0, 30.0

t = [0:100];
y1 = 7*exp(-t/10);
y2 = 2*exp(-t/20);
y3 = 30*exp(-t/5);
r1 = randn(1,101);
r2 = randn(1,101);
r3 = randn(1,101);
r4 = randn(1,101);
z1 = y1 + 0.11*sqrt(y1).*r1;
z2 = y2 + 0.11*sqrt(y2).*r2;
z3 = y3 + 0.11*sqrt(y3).*r1+0.05*r4;


y = z1+z2+z3;
yy = y1+y2+y3;

plot(t,yy,t,y,t,y1,t,y2,t,y3)

save daten t y