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

F = -u1 + u2 + 4*t ;