function [Y] = findExact( t )

%Y = (2.0*t + 1) / (t*t + 2.0);

Y = (t^3)/3.0 + exp(-2.0*t);