function [F] = f1( t, u1, u2 )

F = t^2 - 2.0*exp(-2.0*t) ;