function dx = edsODE(t,x)
% EDSODE  Computes the RHS of a system of 3 first-order ODES which came
% from a scalar 3rd order ODE.  x(1) = y, x(2) = v, x(3) = w

dx = zeros(3,1);
dx(1) = x(2);
dx(2) = x(3);
dx(3) = - x(1) + 5 * x(2) - 5 * x(3);