function xdot = pend(t,x)

  alpha = 0.1;

  Theta  = x(1);
  Lambda = x(2);

  Thetadot = Lambda;
  Lambdadot = -sin(Theta) - alpha * Lambda;

  xdot = [Thetadot; Lambdadot];
