% DGL-Funktion für ode45: 
function  dy_dt = Zerfall_Fun( t, y, zk )   
   dy_dt(1,1) = - zk * y(1);
end   