function [F] = f( t, y )

F = ( 3.0*t*t + 4.0*t + 2.0 ) / ( 2.0 * (y-1.0) );