function dy = ODEPend(t,y)
  k  = 1; alpha = 0.1;
  dy = [y(2);-k*y(1)-alpha*y(2)];
