function f = gv(t,u)

% stiff ode   du/dt = (1 - 4*t/3)*u

f = (1 - 4*t/3)*u;
