function xdot = pendulum0(x,t)
  % k = g / l
  k=1;
  xdot(1) = x(2);
  xdot(2) = -k * sin(x(1));
endfunction
