% exact solution y(t) = cos(t) + (y0 - 1)*exp(-100t)
%       with ic y(0)=y0

function yp = scalarStiffTestOde(t,y)

      yp = -100*(y - cos(t)) - sin(t);