function [U1] = u1Exact( t )

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