McCullagh_second_order_log_l_wrt_psi_alpha function

Second derivative of log(likelihoood) wrt ps[i1, j1] and alpha[index].