function [U2] = u1Exact( t )

U2 = 0.5*exp(2*t) + (t^2) - 0.5;