function dy = pendulumode(t, x)
%PENDULUMODE ODE for simple pendulum
%   DY = PENDULUMODE(X) defines the two coupled ODEs describing
%   the simple pendulum problem.
%
%   T is unused
%   X(1) is theta
%   X(2) is omega
%   DY(1) is d(theta)/dt
%   DY(2) is d(omega)/dt

L = 1;
dy = [x(2); -9.81 / L * sin(x(1))];