function [Y] = findExact( t )

Y = t^2 + (1.0/3.0)*exp(-20.0*t);