function f = fsonderfallstyv( y)
global k1 k2
f = [-k1*y(1) k1*y(1)-k2*y(2)]'; 