McCullagh_second_order_log_l_wrt_psi_delta function

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