function f = fsonderfall(t, y)
k1 = 3; 
k2 = 0.5;
f = [-k1*y(1) k1*y(1)-k2*y(2) k2*y(2)]'; 