function xdot = varcoef3rhs(t,x)
    xdot = zeros(2,1);
    xdot(1) = x(2);
    xdot(2) = -10/(t^2+9)*x(1);
end