function [FP] = fPrime( t, y )

FP = -( 3.0*t*t + 4.0*t + 2.0 ) / ( 2.0 * (y-1.0)^2 );