function ukk = rk4_inline(tk,uk,omega,gamma,Dt)

finline = inline('[0,1;-omega^2,-2*gamma]*uk','tk','uk','omega','gamma');

k1 = finline(tk,uk,omega,gamma);
k2 = finline(tk,uk+0.5*Dt*k1,omega,gamma);
k3 = finline(tk,uk+0.5*Dt*k2,omega,gamma);
k4 = finline(tk,uk+Dt*k3,omega,gamma);

ukk = uk + Dt*(k1/6+k2/3+k3/3+k4/6);