function xdot = series3rhs(t,x)
    xdot = zeros(2,1);
    xdot(1) = x(2);
    xdot(2) = -2*x(1)/(t^2+3);
end