function xdot = ex1rhs(t,x)
    m = 2;
    k = 4;
    k1 = -6;
    k2 = -4;
    
    xdot = [0 1; -k/m 0]*x + [0; (k1*x(1)+k2*x(2))/m];

end