%
% Exact solution for problem
%
% y'''(t) = -t - y(t)
% y(0)   = y'(0) = 0
% y''(0) = 1
%
function res = exact( t )
  res = -t + (2/sqrt(3))*exp(t/2).*sin((sqrt(3)/2)*t) ;
end
