function u = funz(tk,uk,omega,gamma)

u = [0, 1; -omega^2, -2*gamma]*uk;