function [F] = f( t, y )

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