function xdot = prog1013(t,x)
    xdot = zeros(2,1);
    xdot(1) = x(2);
    xdot(2) = cos(t) - 8*x(1) - 3*x(2);
end