McCullagh_second_order_log_l_wrt_psi_delta_vec function

Second derivative of log(likelihood) wrt psi[i1, j1] and delta_vec[k].