function xdot = series2rhs(t,x)
    xdot = zeros(2,1);
    xdot(1) = x(2);
    xdot(2) = -x(1)/(t+1);
end