function [F] = f( t, y )

F = -20.0*(y-t^2)+2*t;