function [F] = f( t, y )

F = 1 + (t-y)^2;