function [U] = u1Exact( t )

U = (7.0/4.0)*t + (1.0/2.0)*(t^3)*log(t) - (3.0/4.0)*(t^3);