McCullagh_second_order_pi_wrt_alpha_c function

Second derivaitve of pi[i, j] wrt alpha[index] and c.