gradf_eps function

Model linearization with respect to epsilon.