function F = ocfcn(t,y)
F = 1-2*t*y(1,:) ;
% This function supplies the 
%   the solution satisfying y(0)=0 is
%   
%       y(t) = exp(-t^2) \int_0^1 exp(s^2) ds