function [F] = f( t, y )

F = 1 + y/t;