LinMatrixH function

Model linearization with respect to epsilon.