LinMatrixLH function

Model linearization with respect to epsilon and eta.