function dx = deriv(t,x)

xi = 0.5;
A = [0 1; -1 -2.0*xi];
B = [0; 1];

dx = A*x + B;