function [F] = f( t, y )

%F = (2.0 - 2.0*t*y)/(t*t + 1.0);
F = t^2 - 2.0*exp(-2.0*t) ;