added dRlm, d2Rlm, dIlm, d2Ilm.  Still need deriv of interaction matrix