function xdot = varcoef1rhs(t,x)
    xdot = zeros(2,1);
    xdot(1) = x(2);
    xdot(2) = -t*x(1);
end
