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

F = u1 - u2 + 2;
