score_matrix function

Calculate matrix of log Likelihood derivatives